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Abstract 


The  purpose  of  this  project  was  to  develop  a  high  repetition 
rate  CO^  TEA  laser.  The  laser  was  pumped  by  a  Marx  generator  and  employed 
spark  gaps  for  the  dual  purpose  of  preionizing  the  gas  discharge  and 
switching  the  Marx  generator.  Typical  operating  characteristics  of  the 
laser  were  studied.  Discharge  instabilities  proved  to  be  a  major  problem 
and  limited  repetition  rates  to  lower  values  than  were  initially  hoped  for. 

The  discharge  instabilities  were  then  studied  in  more  detail  for 
the  purpose  of  identifying  their  causes.  Optical  measurements  of  discharge 
generated  shock  waves  in  the  laser  showed  that  they  were  not  the  cause  as 
they  had  sufficient  time  to  dissipate  between  pulses  and  so  would  not  affect 
following  pulses.  Subsequent  tests  in  a  smaller  discharge  system  showed 
that  the  major  cause  of  instability  in  the  high  repetition  rate  system 
was  due  to  not  allowing  sufficient  convection  time  for  discharge  products 
to  clear  out  of  the  discharge  region  between  pulses.  The  repetition  rate 
drops  off  with  increasing  energy  density  and  Dzakowic  and  Wutzke  have 
proposed  a  theory  to  explain  this  effect.  However  in  this  study  the 
repetition  rate  fell  off  more  quickly  than  their  theory  predicts,  largely 
due  to  the  system  geometry.  That  is,  when  the  energy  density  gets  large 
enough  the  discharge  products  must  clear  not  only  the  discharge  electrodes 
but  the  spark  gaps  as  well. 

Gas  poisoning  was  also  observed  in  this  study  and  caused  arcing 
in  the  closed  cycle  high-repetition-rate  CO^  TEA  discharge  after  only  a  few 
minutes . 

A  comparison  was  done  between  a  spark  pin  preionized  discharge  con¬ 
figuration  with  solid  electrodes  and  a  dielectric  discharge  preionized 
configuration  with  a  screen  electrode.  Higher  repetition  rates  could  be 
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attained  with  the  solid  electrode  system  than  with  the  screen  electrode 
system.  This  was  probably  due  to  discharge  products  getting  trapped  behind 
the  screen  and  causing  discharge  deterioration.  Both  systems  showed 
poisoning  after  about  the  same  length  of  time. 

Some  preliminary  work  was  done  using  a  streak  camera  and  streak 
photographs  of  repetition  rate  induced  instabilities  were  obtained. 
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CHAPTER  1 


INTRODUCTION 

HISTORICAL  SURVEY 

Laser  action  in  a  CO^  gas  discharge  was  first  reported  by 
C.K.N.  Patel  in  1964  (1,2).  Average  power  output  by  the  laser  was  about 
lmW  at  10.4ym  and  9.4ym,  both  continuous  wave  and  pulsed. 

The  next  major  advances  were  the  addition  of  N^  and  He  to  the  CO^ 
resulting  in  large  increases  in  power  output  (3, 4, 5, 6).  Collisions 
between  vibrationally  excited  N^  molecules  and  ground  state  CO^  molecules 
proved  to  be  a  very  effective  means  of  exciting  the  CO^  molecules  to  the 
upper  laser  level.  Helium  aids  the  laser  transition  by  rapidly  depopulating 
the  lower  level  of  the  transition,  thereby  preventing  a  build-up  of 
molecules  in  the  lower  level  which  would  act  as  a  bottleneck  to  transitions 
from  the  upper  level. 

The  use  of  fast  flowing  gas  in  CO^  lasers  to  provide  convective 
cooling  was  reported  in  1969  (7,8,9).  This  was  much  more  effective  than 
the  diffusion  cooling  relied  upon  in  sealed  and  slow  flow  systems  and 
resulted  in  an  increase  in  output  power  per  unit  length  of  several  orders 
of  magnitude. 

Early  CC>2  lasers  ran  at  low  pressures.  High  pressure  lasers  were 
developed  in  order  to  provide  a  higher  density  of  molecules  capable  of 
lasing  and  hence  higher  output  power.  In  addition  large  volume  atmospheric 
pressure  lasers  are  easier  to  construct  than  large  volume  low  pressure 
lasers  because  they  do  not  have  to  be  as  well  sealed  or  capable  of  with¬ 
standing  as  great  pressure  differences. 

A  CC>2  laser  works  best  with  an  E/P  ratio  of  about  10  to  50  V/ (cm  torr)  ? 

where  E  is  electric  field  and  P  is  gas  pressure.  So  as  gas 
pressures  rose  this  meant  working  with  very  high  voltages  for 
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longitudinally  excited  plasmas.  To  avoid  this  problem  Transversely 
Excited  Atmospheric  (TEA)  lasers  were  developed  (10,11)  which  use  large 
surface  area  electrodes  placed  close  together  to  excite  a  large  volume 
of  plasma  with  much  lower  voltages  than  are  needed  for  longitudinal 
excitation.  Since  the  TEA  laser  electrodes  are  closer  together,  lower 
voltages  will  produce  just  as  large  fields. 

Another  major  advancement  in  the  construction  of  TEA  lasers  was 
the  idea  of  conditioning  the  gas  with  some  form  of  preionization  in  order 
to  obtain  a  stable,  uniform,  large  volume  plasma  more  easily  (12).  Pre¬ 
ionization  has  been  accomplished  by  a  number  of  techniques  including 
the  use  of  UV  (13),  plasma  jets  (14),  electron  beams  (15),  neutrons  (16), 
and  organic  gas  additives  (17) . 

Most  TEA  laser  development  to  date  has  gone  into  high  energy  single 
pulse  CC^  lasers.  However,  a  few  high  repetition  rate  systems  have  been 
built  (18-22) .  The  laser  described  in  reference  21  is  not  a  CO^  laser 
but  the  authors  have  used  high-repetition-rate  CC^  TEA  laser  technology, 
only  substituting  rare  gases  for  the  CO^  laser  mixture.  These  systems 
all  use  gas  convection  to  clear  undesirable  discharge  products  from  the 
discharge  region  between  pulses.  The  limitation  on  maximum  repetition 
rate  has  been  due  to  the  onset  of  discharge  arcing,  and  this  arcing 
usually  occurs  at  lower  repetition  rates  than  would  be  expected  from  the 
velocity  of  the  gas  flowing  between  the  electrodes  (18,19).  The  best 
explanation  of  this  effect  offered  so  far  is  that  of  Dzakowic  and  Wutzke 
(18) .  Another  problem  which  can  occur  in  a  closed-cycle  system  is 
degradation  of  discharge  stability  due  to  dissociation  of  the  CO^  into 
CO  and  0^  (22)  .  A  solution  to  the  problem  of  C02  dissociation  is  the 
use  of  some  catalyst  to  recombine  the  CO  and  0^  back  into  CO^  (23) . 
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1.1  SUMMARY  OF  EXPERIMENTS 

In  this  thesis  the  design  and  testing  of  a  high  repetition  rate 
closed  cycle  CO^  TEA  laser  is  described.  Although  the  laser  did  operate 
at  high  repetition  rates,  discharge  arcing  limited  it  to  a  lower 
repetition  rate  than  was  initially  estimated.  A  poisoning  effect  due  to 
the  dissociation  of  the  laser  gases  was  also  a  suspected  source  of 
problems.  In  order  to  understand . the  causes  of  the  arcing,  the  discharge 
was  studied  in  further  experiments. 

To  facilitate  the  discharge  studies  a  smaller  more  flexible 
discharge  system  was  designed  which  modeled  the  essential  features  of 
the  laser  discharge  system.  Using  this  system  a  number  of  key  parameters 
that  could  contribute  to  discharge  arcing  were  studied.  These  parameters 
were  gas  flow  velocity,  pulse  length,  discharge  energy  density,  system 
geometry  and  preionization  technique.  Discharge  afterglows  were 
monitored  using  a  photomultiplier  to  observe  convective  clearing  of  the 
discharge  products  from  the  discharge  region  and  to  observe  discharge 
expansion,  an  effect  which  was  predicted  by  Dzakowic  and  Wutzke  (18) . 

A  streak  camera  was  also  used  in  a  attempt  to  monitor  discharge  expansion. 
Although  the  streak  camera  failed  to  show  the  discharge  expansion  it 
did  produce  some  streak  photographs  which  illustrate  the  various  types 
of  arcing  that  occur. 

In  studying  the  effects  of  system  geometry  and  preionization  on 
arcing  three  different  techniques  were  compared.  They  were  solid  electrodes 
with  spark  pin  preionization  both  upstream  and  downstream,  solid  electrodes 
with  just  downstream  preionization,  and  a  system  where  preionization  was 
provided  by  an  auxilliary  discharge  behind  a  screen  electrode. 
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Besides  the  parameters  studied  on  the  small  discharge  system, 
the  presence  of  shock  waves  was  investigated  on  the  laser  system  to 
see  if  they  contributed  in  any  way  to  discharge  arcing. 

1.2  LASER  SYSTEM 

This  project  started  with  the  design  and  testing  of  a  closed  cycle 
high  repetition  rate  CC^  TEA  laser.  The  laser  had  three  basic  subsystems: 
a  closed  cycle  recirculating  gas  system,  the  electrical  discharge 
system,  and  the  optical  system. 

1.3  RECIRCULATING  GAS  SYSTEM 

The  recirculating  gas  system  is  a  small  scale  wind  tunnel  which 
circulates  the  laser  gas  mixture  through  the  electrical  discharge  region. 
It  consists  of  a  centrifugal  blower,  flow  ducting,  and  a  heat  exchanger. 
The  layout  of  the  components  is  illustrated  in  Fig.  1.1.  The  blower 
was  capable  of  delivering  a  volumetric  gas  flow  of  1500  Ji/sec,  which 
corresponded  to  a  velocity  of  100  m/sec  between  the  electrodes. 

1.4  ELECTRICAL  DISCHARGE  SYSTEM 

The  electrical  discharge  system  used  to  pump  the  laser  is  illustrated 
in  Figs.  1.2  and  1.3.  It  is  a  Marx  generator  with  its  switching  gaps 
located  downstream  from  the  discharge  electrodes  in  the  flow  duct. 

Locating  the  switching  gaps  in  the  flow  duct  has  a  dual  purpose  -  to 
extinguish  the  switching  gaps  between  pulses  using  gas  flow  and  to  pre¬ 
ionize  the  discharge  using  ultraviolet  light  produced  by  the  spark  gaps. 

A  time  delay  is  needed  between  the  spark  gaps  firing  and  the  formation  of 
the  discharge  in  order  for  the  preionization  to  take  effect.  This  time 
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Fig.  1.1  Physical  Layout  of  C02  TEA  Laser  System. 

The  system  is  7  feet  high,  4  feet  wide  and  3  feet  deep. 
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Fig.  1.2  Schematic  of  Marx  Generator.  Typical  component  values  are 
Cp  =  . 03yf ,  C  =  . 3yf ,  and  L  =  30mh. 
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to  additional 
spark  pins 


Fig.  1.3  Schematic  of  Marx  Triggering  Circuit.  Pins  1,  3,  5,  7 
and  9  of  the  Marx  are  triggered. 
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delay,  and  a  certain  amount  of  pulse  shaping  as  well,  is  supplied  by 
the  peaking  capacitor.  In  order  for  the  preionization  to  be  effective 
the  Marx  spark  gaps  must  fire  simultaneously  and  to  accomplish  this  it 
was  necessary  to  trigger  more  than  one  of  the  spark  gaps.  A  minimum  of 
3  gaps  have  to  be  triggered  but  typically  5  are.  A  schematic  of  the 
trigger  circuit  used  is  shown  in  Fig.  1.3. 

The  discharge  electrodes  are  61  cm  long  and  10  cm  wide.  They  are 
separated  by  2.5  cm  to  give  a  discharge  volume  of  approximately  0.37£. 

1.5  OPTICAL  SYSTEM 

Both  stable  and  unstable  resonators  were  tested  as  optical  systems 
for  the  laser.  Internal  windows  on  the  laser  cavity  were  also  tested 
but  introduced  too  much  loss  into  the  low  gain  cavity,  severely  reducing 
power  output.  As  a  result  it  was  necessary  to  attach  the  optics  directly 
to  the  laser  cavity  with  bellows.  This  resulted  in  problems  with 
vibration  being  transmitted  to  the  optics,  causing  unstable  mode  output. 
The  unstable  resonator  provided  a  greater  power  output  than  the  stable 
resonator . 


CHAPTER  2 


SYSTEM  DESIGN 
2.1  FLOW  SYSTEM  DESIGN 

The  laser  gas  was  contained  and  circulated  in  a  large  sealed  system. 
Although  not  a  high  vacuum  system,  it  was  reasonably  airtight  and  could 
be  pumped  down  to  a  pressure  of  1  or  2  torr.  After  being  pumped  out  the 
system  was  backfilled  to  the  appropriate  partial  pressure  for  each  gas 
depending  on  the  desired  laser  mixture. 

Gas  flow  through  the  system  was  driven  by  a  centrifugal  fan  which 
was  powered  by  a  15  HP  electric  motor.  In  order  to  avoid  building  a 
complicated  shaft  seal  the  motor  was  installed  in  the  fan  box  along  with 
the  fan.  A  heat  exchanger  was  also  provided  to  remove  heat  introduced 
into  the  system  by  the  discharge  and  fan. 

A  sealed  recirculating  system  such  as  this  is  very  similar  to  a 
closed  cycle  subsonic  wind  tunnel  and  so  was  designed  according  to 
principles  outlined  in  Pope  (24)  and  Pankhurst  (25)  for  subsonic  wind 
tunnels.  Two  other  sources  which  proved  useful  were  the  Trane  Air 
Conditioning  Manual  (26)  which  contained  useful  information  on  pressure 
measurement  and  fan  characteristics;  and  the  ASHRAE  Handbook  (27)  which 
provided  information  on  ducting  design.  Essentially  the  design  criterion 
consisted  of  trying  to  offer  as  little  resistance  to  gas  flow  in  the 
ducting  as  possible,  in  order  to  achieve  a  high  flow  velocity. 

In  a  recirculating  gas  system  the  gas  flows  under  the  influence  of 
a  pressure  head  developed  by  the  system  fan.  Opposing  the  gas  flow  is 
skin  friction  resulting  from  the  movement  of  gas  through  the  ducting. 
Besides  the  head  loss  due  to  skin  friction  considerable  head  loss  occurs 
due  to  changes  in  ducting  cross-sectional  area  and  changes  of  direction 
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of  the  duct.  These  losses  are  proportional  to  velocity  squared  and  are 
described  by  eqn.  2.1.1 


AU  1  2 
AH  =  —  cpu 


2.1.1 


AH  is  the  pressure  head  loss,  p  is  the  gas  density,  u  is  the  gas  velocity, 
and  c  is  a  friction  coefficient. 

Initially  the  gas  flow  is  converged  from  37  cm.  wide  to  2.5  cm. 
wide  between  the  electrodes.  To  reduce  friction  losses  wide  return 
passages  were  used.  Converging  the  gas  flow  is  very  efficient  and 
produces  little  head  loss.  However  expanding  the  gas  flow  for  the  return 
trip  is  much  more  lossy.  In  order  to  minimize  this  loss  diffusers  are 
designed,  where  space  permits,  with  small  angles  between  the  walls  because 
larger  angles  cause  flow  separation  which  is  very  lossy.  Such  diffusers 
are  called  narrow  angled  diffusers  and  have,  typically,  7°to  10°  angles 
between  the  walls.  Wider  angled  diffusers  usually  make  use  of  vanes  to 
help  reduce  losses  and  the  angle  between  the  vanes  is  kept  around  7° to 
10°.  A  10°  angle  diffuser  that  expanded  from  2.5  cm.  to  20  cm.,  the 
width  of  the  return  passage,  would  be  very  long.  In  order  to  conserve 
space  a  wide-angle  vaned  diffuser  was  designed  according  to  the  design 
criterion  outlined  in  a  paper  by  Feil  (28) . 

Using  eqn.  2.1.1  an  estimate  of  system  performance  was  obtained. 

In  this  estimate  only  the  main  components  contributing  to  system  head 
loss  were  considered.  They  are:  the  throttle,  or  region  between  the 
electrodes,  the  wide  angled  diffuser,  the  two  elbows,  and  the  long  narrow 
angled  diffuser  that  acts  as  a  return  duct.  All  the  estimates  can  be 
related  to  the  velocity  between  the  electrodes  by  use  of  the  continuity 
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equation.  This  just  means  that  the  volumetric  flow  of  gas  is  the  same 
through  all  parts  of  the  system,  VA  =  constant.  So,  for  instance,  the 

cross  sectional  area  in  the  elbows  is  8  times  that  between  the  electrodes 

♦ 

so  the  velocity  in  the  elbows  will  be  1/8  that  between  the  electrodes. 

The  head  loss  for  the  various  parts  of  the  system  are  summarized  in 
table  2.1. 


c 

V 

!  ^ 
hpv 

high  velocity  region 

.1 

V 

.050pV2 

wide  angled  diffuser 

.35 

v/2 

. 044pv2 

elbow 

.3 

v/8 

.002pv2 

elbow 

.3 

v/8 

.002pv2 

long  diffuser 

.3 

v/8 

.002pv2 

total 

O.lOpv2 

Table  2.1 

Pressure  Head  Loss  Estimates  Throughout  the 
Flow  System 


The  values  of  c  are  typical  values  and  are  obtained  from  references 
25,  27  and  28.  The  values  of  v  for  the  diffusers  are  the  velocities 
at  their  inlets.  It  can  be  seen  that  the  main  sources  of  pressure  head 
loss  are  the  high  velocity  region  between  the  electrodes  and  the  wide 
angled  diffuser  because  of  the  high  gas  velocities  in  them.  Losses  in 
the  elbows  and  long  diffuser  are  almost  negligible.  The  high  velocity 
region  accounts  for  half  the  head  loss  by  itself. 
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The  total  head  loss  then,  as  a  function  of  velocity,  is  given  by 

2 

AH  =  O.lpv  2.1.2 

The  cross  sectional  area  between  the  electrodes  is  l"x24",  and  using  this, 

the  head  loss  can  be  converted  to  a  function  of  flow  rate  in  CFM. 

3 

Standard  air  density  is  0.075  lb (mass) /ft  .  Using  these  values,  and  making 
suitable  transformations  to  get  the  correct  units  this  becomes 

AH  =  4.33x10_7Q2  2.1.3 

where  AH  is  pressure  loss  in  inches  of  H^O  and  Q  is  flow  rate  in  CFM. 

This  is  plotted,  in  Fig.  2.1,  on  the  manufacturer  supplied  characteristics 
for  the  fan  used  in  this  system.  The  fan  is  a  13^  inch  diameter  centrifugal 
blower  running  at  4628  rpm.  It  can  be  seen  that  the  system  losses  and 
fan  static  pressure  match  at  a  gas  flow  rate  of  550  CFM  which  corresponds  to 
a  flow  velocity  of  160  m/sec.  To  achieve  this  velocity  in  air  would  require 
20  HP  but  fortunately  considerably  less  power  is  needed  in  the  laser 
mixture  because  the  power  is  proportional  to  the  gas  density.  Since  a 
typical  laser  mixture  is  mostly  He  (typically  80%  He,  10%  N  ,  10%  CO^)  with 
an  atomic  weight  of  2  and  air  is  mostly  with  an  atomic  weight  of  14,  the 
laser  mixture  will  have  roughly  1/7  the  density  of  air  so  will  only  require 
about  3  HP  to  achieve  this  flow  velocity. 

2.2  MARX  GENERATOR  DESIGN 

This  laser  is  essentially  a  double  discharge  design  switched  by 
spark  gaps.  Since  the  spark  gaps  must  be  blown  in  order  to  recover  quickly 


Static  Pressure  inches  of  water 
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fan  speed  4628  rpm 
gas  density  .075  lbs/ft 


Fig.  2.1  Fan  Characteristic  and  System  Pressure  Head  Loss.  SP  is 

the  fan  static  pressure  curve,  HL  is  the  system  head  loss, 
and  HP  is  the  power  required  to  drive  the  fan.  The  static 
pressure  and  system  loss  curves  intersect  at  approximately 
5500  CFM  and  13  inches  of  water. 
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for  high  repetition  rate  operation  it  was  decided  to  include  them  in  the 
same  flow  system  as  the  discharge  electrodes  which  also  rely  on  convection 
to  remove  discharge  products  from  between  them.  By  using  a  Marx  generator, 
which  is  a  voltage  multiplying  circuit,  any  desired  pulse  voltage  can 
be  achieved  for  optimum  performance  and  in  addition  the  Marx  switching  gaps 
serve  as  a  source  of  ultraviolet  light  to  preionize  the  discharge  region. 
This  then  combines  the  two  essential  features  of  discharge  preionization  and 
Marx  switching  into  one  set  of  spark  gaps,  arranged  along  the  length  of 
the  electrodes. 

Design  of  a  Marx  generator  is  a  straight  forward  procedure  and  is 
outlined  in  references  29  and  30.  A  schematic  of  the  Marx  used  for  this 
laser  is  shown  in  Figs.  1.2  and  1.3.  A  Marx  generator  consists  of  a  number 
of  capacitors  which  are  charged  up  in  parallel  through  hold-off  inductors, 
then  discharged  in  series  through  the  spark  gaps  to  provide  voltage  multi¬ 
plication.  The  charging  network  can  be  either  resistors  or  inductors 
but  for  a  high  repetition  rate  system  the  power  loss  in  resistor  charging 
would  be  prohibitive  so  inductors  are  used  instead.  For  a  N  stage  Marx 
the  output  voltage  will  be  approximately  N  times  the  Marx  charging  voltage. 
If  the  hold-off  inductors  are  too  small  or  the  time  lag  in  the  spark  gaps 
firing  is  too  long, enough  charge  can  be  lost  through  circulating  currents 
in  the  Marx  to  reduce  the  output  voltage  appreciably. 

The  parameters  of  this  Marx  generator  were  determined  largely  by 
experimentation.  However  an  estimate  of  the  discharge  time  through  the 
hold-off  inductors  was  needed  to  determine  that  they  would  not  shunt  out 
the  discharge. 

The  current  in  an  LC  circuit  is  given  by 


•  MW  « 


•  PP 
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and  so  has  period 

x  =  2tt/lC 

The  smallest  discharge  loop  in  this  circuit  would  consist  of  two  capacitors 
and  one  inductor  in  series.  Initially C  was  chosen  as  . lpf  and  L  as  30mh. 
This  gave  a  discharge  time  of 

i  =  2i\\f  30x10  ^  x  xlO  sec 
-4 

=  2.4x10  sec 

Since  the  discharge  pulse  length  is  on  the  order  of  1  usee  this  is 
sufficiently  long  that  no  appreciable  charge  is  lost  through  the  hold- 
off  inductors  during  the  discharge  pulse. 


■ 


CHAPTER  3 


SYSTEM  TESTING 


3.1  FLOW  TESTING 

Flow  velocity  was  measured  using  a  pitot  tube  and  a  capacitance 
manometer.  Connecting  the  static  pressure  from  the  tube  to  one  port  of 
the  manometer  head  and  the  total  pressure  to  the  other  port  resulted  in 
a  differential  measurement  that  gave  the  velocity  pressure.  The  velocity 
can  then  be  calculated  from  the  velocity  pressure  using 


2 


AP  =  ifpv 


or 


3.1.1 


2  3 

where  AP  is  pressure  in  N/m  ,  p  is  gas  density  in  kg/m  ,  and  v  is 


in  m/sec. 


Results  of  the  velocity  measurements  are  given  in  table  3.1  and 
are  reasonably  close  to  the  original  design  estimates.  The  average  gas 
flow  velocity  between  the  electrodes  is  approximately  100  m/sec.  The 
design  estimate  was  160  m/sec. 

Flow  uniformity  was  examined  between  the  electrodes,  at  the  exit 
from  the  wide  angled  diffuser,  and  in  the  narrow  angled  diffuser  connecting 
the  two.  Since  these  measurements  were  made  in  air,  and  the  system 
drive  motor  is  only  rated  for  15  HP  it  was  necessary  to  slow  the  fan  down 
to  avoid  overloading  the  motor.  This  was  accomplished  by  exchanging 
the  fan  and  motor  pulleys.  Although  the  velocity  will  be  slower  the 
velocity  distribution  should  not  be  affected.  Fig.  3.1  shows  velocity 
distributions  between  the  electrodes  and  at  the  exit  from  the  wide  angled 
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diffuser.  The  flow  velocity  between  the  electrodes  is  uniform  but  the 


Air  730  torr 


distance  from  endwall 
inches 

0 

2 

4 

6 

8 

10 

velocity  pressure 
torr 

42 

49 

50 

52 

52 

50 

AP  =  50  torr 
av 

velocity  =  105  m/sec 


laser  mix  70:70:590  torr 
C02 :N2 : He 


distance  from  endwall 
inches 

0 

2 

4 

6 

8 

10 

velocity  pressure 
torr 

17 

19 

19 

19.5 

20 

20 

AP  =19.5  torr 
av 

velocity  =  107  m/sec 

The  average  values  of  velocity  pressure  were  calculated  without 
using  the  values  at  the  endwall  (distance  =  0)  because  those 
values  would  be  affected  by  the  boundary  layer. 


Table  3.1  Velocity  Measurements 
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inches  along  diffuser 
(b) 


Fig.  3.1  Gas  Flow  Velocity  Distributions  (a)  Between  the  Electrodes 


and  (b)  Leaving  the  Wide  Angled  Diffuser. 
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flow  velocity  leaving  the  wide  angled  diffuser  was  very  nonuniform.  The 
dead  spaces  near  the  endwalls  indicate  that  flow  separation  has  occurred, 
and  this  reduces  the  efficiency  of  the  diffuser. 

Referring  back  to  Fig.  1.1,  the  wide  angled  diffuser  is  connected 
to  the  electrodes  by  a  narrow  angled  diffuser  that  increases  the  channel 
width  from  2.5  cm  to  5  cm.  Further  tests  showed  that  the  separation  was 
already  occurring  back  in  the  narrow  angled  diffuser,  and  the  gas  entering 
the  wide  angled  diffuser  was  already  separated  from  the  end  walls.  The 
wide  angled  diffuser  only  accented  the  effect  as  it  slowed  down  the  gas 
still  further.  Because  the  narrow  angled  diffuser  was  constructed  from 
perspex  it  was  transparent  and  a  simple  technique  could  be  used  to  visualize 
the  flow  separation  occurring  in  it.  The  method  was  to  tape  short  tufts 
of  yarn,  about  5  cm  long,  to  the  walls  of  the  diffuser.  They  readily 
showed  the  direction  of  flow  and  degree  of  turbulence  in  the  diffuser.  A 
resulting  flow  map  is  shown  in  Fig.  3.2.  The  electrodes  were  to  a  large 
extent  responsible  for  this  separation.  This  was  demonstrated  by  replacing 
the  electrodes  with  flat  plates  and  taping  shut  the  spark  pin  holes  in 
the  diffuser.  The  resulting  flow  map  is  shown  in  Fig.  3.3  and  shows  a 
considerable  improvement  over  the  flow  with  the  electrodes  in  place.  Since 
the  electrodes  are  a  necessary  part  of  the  system,  flow  separation  caused 
by  them  must  be  tolerated.  An  attempt  was  made  to  reduce  the  separation 
by  adding  cross  vanes  to  the  narrow  angled  diffuser  to  act  as  flow 
straighteners .  But  since  no  improvement  in  flow  velocity  was  observed 
with  the  vanes  they  were  removed. 

Although  better  diffuser  design  should  result  in  reduced  flow 
separation  losses  and  hence  a  higher  velocity  the  system  is  adequate  for 
the  purposes  of  this  experiment.  The  flow  between  the  electrodes  meets 
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The  dotted  lines  are  flow  lines  and  the  shaded  areas  are  turbulent  regions. 
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The  dotted  lines  are  flow  lines  and  the  shaded  areas  are  turbulent 
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the  requirements  of  being  uniform  and  fast. 

Another  potential  problem  considered  to  be  a  possibility  during 
laser  operation  was  the  slowing  down  of  gas  flow,  and  hence  a  reduction 
in  attainable  repetition  rate,  due  to  shock  waves  from  the  discharge. 

To  check  this  a  run  was  made  with  a  pitot  tube  inserted  into  the  return 
duct  to  record  velocity  pressure  while  the  laser  was  operating  and  the 
velocity  pressure  was  observed  on  an  inclined  U-tube  manometer.  The  U 
tube  manometer  was  used  because  the  capacitance  manometer  could  not  operate 
in  the  noisy  electrical  environment  of  the  running  laser.  No  change  in 
velocity  pressure  was  observed  with  a  change  in  discharge  repetition  rate 
from  0  to  200  Hz,  so  flow  choking  due  to  shock  waves  does  not  occur  in 
this  system. 

3.2  MARX  GENERATOR  TESTING 

Many  of  the  parameters  of  this  Marx  generator  were  determined 
by  experimentation.  They  were  the  number  of  stages,  triggering  method, 
spark  pin  location,  spark  pin  material,  and  the  need  for  a  peaking 
capacitor . 

Originally  a  6  stage  Marx  was  designed  but  5  and  10  stages  were 
tested  as  well,  10  stages  being  chosen  for  the  final  design.  The  10 
stage  Marx  ran  more  stably  than  the  5  stage  Marx  of  the  same  capacitance. 
This  was  because  the  spark  gaps  could  more  easily  hold  off  the  lower 
voltage  of  4  to  6  KV  needed  for  the  10  stage  Marx  than  the  lh  to  8  KV 
needed  for  the  5  stage  Marx.  The  discharge  was  also  more  uniform  for 
the  10  stage  Marx  due  to  the  increased  preionization. 

Three  important  results  were  observed  with  respect  to  the  Marx 
operation  and  the  quality  of  the  discharge.  Firstly  the  Marx  spark 
pins  must  be  downstream  from  the  discharge.  Secondly  adding  a  peaking 
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capacitor  in  parallel  with  the  discharge  helped  stabilize  it.  Finally 
the  method  of  triggering  the  Marx  affected  the  discharge  stability 
because  it  affected  the  Marx  commutation  time. 

Originally  the  Marx  was  designed  with  its  spark  pins  upstream  from 
the  discharge.  However  moving  the  pins  downstream  greatly  reduced  the 
amount  of  discharge  arcing  that  occurred.  This  implies  that  discharge 
products  from  the  spark  gaps  being  swept  into  the  main  discharge  region 
was  a  cause  of  arcing. 

The  addition  of  a  peaking  capacitor  (C  in  Fig.  1.2)  in  parallel 
with  the  main  discharge  helped  stabilize  it  and  reduce  arcing.  Small 
peaking  capacitors,  typically  1000  pf  to  5000  pf,  were  used.  Originally 
it  was  assumed  that  the  peaking  capacitor  introduced  a  time  delay  needed 
for  the  preionization  to  take  effect.  This  explaination  is  not  likely 
however,  because  the  UV  only  needs  on  the  order  of  10  ^  sec.  to  reach 
the  discharge  from  the  spark  pins.  A  more  likely  explanation  is  that  the 
peaking  capacitor  acts  as  a  low  impedance  source  which  allows  the  discharge 
to  form  rapidly. 

The  small  peaking  capacitor  results  in  a  double  humped  current 
pulse  (Fig.  3.4a)  and  in  order  to  understand  this  waveform  a  simple  circuit 
model,  described  in  section  3.3,  was  devised.  Essentially  this  model 
showed  that  the  initial  peak  was  due  to  the  rapid  discharge  of  the  peaking 
capacitor,  followed  by  a  slower  peak  due  to  the  discharge  of  the  residual 
charge  left  in  the  Marx.  This  suggested  that  if  all  the  charge  could  be 
transferred  from  the  Marx  to  the  peaking  capacitor  before  breakdown 
occurred  a  much  cleaner  pulse  would  result.  This  was  most  easily  done 
by  making  the  trigger  capacitor,  C  ,  equal  in  value  to  the  total  capacitance 
of  the  Marx  with  all  stages  in  series.  This  resulted  in  a  cleaner  and 
shorter  current  pulse,  shown  in  Fig.  3.4b. 
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Fig.  3.4  Output  Pulses  for  10  Stage  Marx,  O.lyf  Per  Stage,  Using  Two 


In  (a)  C  =  1750pf  and  in  (b)  C  =  O.Olyf. 
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Initially  low  inductance  capacitors  were  constructed  for  this 
purpose  using  sheets  of  brass  separated  by  glass  dielectric.  Although 
they  worked  well  enough  to  demonstrate  the  technique  they  failed  quickly 
due  to  the  glass  shattering,  either  from  insufficient  dielectric 
strength,  or  from  mechanical  shock  caused  by  rapid  charging  and  dis¬ 
charging.  Eventually  capacitors  were  ordered  for  this  purpose  from 
Condenser  Products. 

Finally,  the  method  of  triggering  the  Marx  was  found  to  be  important 
in  forming  a  good  discharge.  Initially  only  one  gap  of  the  Marx  was 
triggered.  The  literature  (29,30)  indicated  that  this  plus  the  UV  coupling 
between  the  spark  pins  would  be  sufficient  for  fast  switching  of  the 
Marx.  However  the  commutation  time  of  the  Marx  proved  to  be  about  15  ysec. 
(Fig.  3.5).  This  long  time  delay  was  detrimental  to  discharge  formation, 
possibly  because  it  gave  the  preionization  electron  density  time  to 
recombine.  Triggering  several  Marx  gaps  resulted  in  improved  discharges 
and  a  Marx  switching  time  of  less  than  100  nsec.  For  the  10  stage  Marx 
as  many  as  all  10  gaps  were  triggered  but  a  minimum  of  3  triggered 
gaps  was  necessary  for  successful  operation.  In  the  final  design  5  gaps 
were  triggered.  A  disadvantage  with  triggering  a  number  of  gaps  is 
that  they  cannot  all  be  switched  to  ground.  As  a  result  high  voltage 
pulses  feed  back  into  them  when  the  Marx  switches  and  it  is  difficult 
to  build  a  reliable  high  repetition  rate  trigger  pulser  capable  of 
surviving  this.  As  a  result  pulser  failures  were  a  continual  problem. 

Fig.  3.6  shows  the  construction  of  a  trigger  pin  used  in  the  Marx.  A 
pulser  unit  that  used  an  HY60  thyratron  was  fairly  resistant  to  damage 
(Fig.  3.7). 
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5  |j  s  e  c/d  i  v 

Fig.  3.5  Commutation  Delay  in  the  Marx  Generator  with  a  Single  Trigger  Pin. 
The  upper  trace  is  the  voltage  at  the  first  Marx  pin  and  the 
lower  trace  is  the  Marx  output  voltage.  These  waveforms  are 
for  a  six  stage  Marx  generator.  The  first  pin  switches  at  the 
beginning  of  the  square  pulse  and  the  output  pulse  occurs  at 
the  end  of  the  square  pulse  approximately  15  ysec.  later. 
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Figure  3.7  HY60  Pulse  Generator  Used  as  a  Trigger  Unit  for  the  Marx  Generator 
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For  the  Marx  pins  carbon  and  brass  rods  were  both  tried.  The 
brass  rods  were  chosen  for  the  final  design  because  they  were  more 
durable  and  caused  less  discharge  instability  than  the  carbon  rods. 


3.3  CIRCUIT  MODEL 

In  order  to  understand  the  waveform  resulting  from  the  Marx 
generator  and  peaking  capacitor  combination  a  simple  circuit  model,  shown 
in  Fig.  3.8,  was  devised.  The  discharge  load  was  represented  by  a  switch 
in  series  with  a  resistor.  The  Marx  generator  was  modelled  by  a  capacitor, 
inductor,  and  switch  which  represented  the  Marx  spark  gaps.  C^  represents 
the  peaking  capacitor.  Initially  C ^  is  at  0  volts  and  C^  is  at  voltage 
V  .  Closing  at  t=0  is  equivalent  to  the  Marx  gaps  firing  and  C^  then 
charges  from  C^. 

Using  Laplace  transforms  voltage  and  current  waveforms  can  be  solved 
for  this  model.  The  initial  conditions  at  t=0  are  shown  in  Fig.  3.9a, 
and  the  circuit  equation  is 
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The  inverse  transform,  from  tables,  is 
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Fig.  3.8  Circuit  Model. 
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Fig.  3.9  Initial  Conditions  for  the  Circuit  Model  of  Fig.  3.8 


at  (a)  t=0  and  (b)  t=tB* 
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The  negative  sign  comes  from  the  direction  around  the  loop  for  which 
voltage  drops  are  taken  to  be  positive. 
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At  the  time  of  breakdown,  tfi ,  switch  closes.  This  gives  a 
new  circuit,  shown  in  Fig.  3.9b,  with  the  following  intial  conditions 
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Taking  t  to  be  t=0  for  the  new  circuit  gives  two  loop  equations. 
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Finally,  combining  3.3.5  and  3.3.6  gives 
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The  solution  i2(t)  is  the  inverse  transform  of  this  expression. 

Since  the  denominator  is  a  cubic  equation  there  are  four  possible  solutions, 
depending  on  Its  roots.  If  the  roots  are  real  and  unequal  the  solution 
will  have  the  form  i(t)  =  Ae  +  Be  +  Ce  1  where  -a,  -3,  -y  are  the 
roots  of  the  equation.  If  the  roots  are  real  and  two  are  equal  the 
solution  will  have  the  form  i(t)  =  Ae  at  4-  (Bt  +  C)  e  ^ .  If  the  roots 
are  real  and  all  three  are  equal  the  solution  will  have  the  form 

2  —  ct  t 

i(t)  =  (At  +  Bt  +  C)e  .  Finally  if  one  root  is  real  and  two  complex 


conjugates,  the  solution  has  the  form  i(t)  =  Ae  Be  cos(3t-^)  where  the 

roots  are  -y,  -a-Sj ,  -a+Sj . 

Of  these  possibilities  only  the  last  one  can  give  a  shape  similar 
to  the  observed  waveform.  Initially  then,  the  voltage  across  the 
electrodes  rises  according  to  eqn.  3.3.2 
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until  the  breakdown  voltage  is  reached.  Note  that  can  be  higher  than 

V  ,  the  initial  Marx  voltage.  If  C  <<C..  breakdown  will  still  occur  for 

V  almost  as  high  as  2V  .  When  breakdown  occurs  the  discharge  voltage 
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has  the  form 
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v(t)  =  Ri(t)  =  Ae  +  Be  at  cos(Bt-iJj) 


3.3.8 


So  the  initial  spike  is  due  to  the  rapid  discharge  of  the  peaking 
capacitor,  followed  by  the  slower  discharge  of  the  residual  charge  left 
in  the  Marx  generator. 

Eqn.  3.3.7  has  the  form 
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where  N^Cs)  =  a^s  +  a^s  +  a 


and  the  denominator  is  a  cubic  equation  with  roots  ~Y ,  -a-Bj ,  -a+Bj .  The 
inverse  transform  of  this  equation  is  (ref.  31). 
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A  small  computer  program  was  written  to  evaluate  the  roots  of  the 
denominator  of  3.3.7,  calculate  the  coefficients  A,  B,  y,  a,  3  and  ip,  and 
then  to  calculate  the  voltage  across  the  load  resistor.  Typical  waveforms 
are  shown  in  Fig.  3.10. 

Although  qualitatively  this  model  offers  insight  into  the  discharge 
waveform  the  quantitative  agreement  is  poor.  The  main  points  of  dis¬ 
agreement  are  that  the  model  predicts  an  intial  current  spike  that  is 
much  too  short,  and  current  values  that  are  too  high. 

A  better  model  is  suggested  in  Fig.  3.11.  The  Laplace  transform 
for  i2(t)  is 
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The  derivation  is  identical  to  that  for  the  simpler  model. 

Since  3.3.9  has  a  4th  order  eqn.  in  the  denominator  it  can  have 
4  roots  with  imaginary  components  -a  +  a j ,  -a  -  aj ,  -b  +  3 j ,  -b  -  3 j . 
This  would  give  a  solution  of  the  form 
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Fig.  3.11  Improved  Circuit  Model. 
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This  model  is  better  because  it  would  give  the  initial  sinusoidal  current 

rise  which  the  other  model  ignores,  and  the  inductance  would  widen  the 
initial  current  spike. 

3. A  PULSE  STRETCHING 

For  atmospheric  pressure  gas  discharges  the  Marx  generator  typically 
produced  current  pulse  lengths  of  between  0.5  and  1  ysec  duration. 

Attempts  were  made  to  increase  this  pulse  length  by  adding  a  series 
inductor  to  the  circuit  (Fig.  3.12a).  A  small  capacitor  shunts  the  inductor 
to  help  start  the  discharge.  For  pulse  lengths  longer  than  6  to  7  ysec 
in  the  1:1:8  CO^iN^tHe  laser  mix  arcing  occurred.  Increasing  the  He  to 
90%  of  the  mixture  permitted  pulse  lengths  of  up  to  12  ysec  before  arcing 
(Fig.  3.13)  but  for  12  ysec  pulses  the  discharge  was  nonuniform. 

At  reduced  gas  pressures  the  discharge  is  more  tolerant  of  long 
pulses.  At  200  torr  with  90%  He  18  ysec  pulses  produce  a  good  discharge 
and  up  to  30  ysec  pulses  could  be  produced. 

For  high  repetition  rate  operation  at  atmospheric  pressure  short 
pulse  lengths  are  required  for  a  stable  discharge  in  this  system. 

Another  serious  problem  occurred  which  made  pulse  stretching  in 
this  system  somewhat  impractical.  With  the  circuit  configuration  in 
Fig.  3.12a,  arcing  would  occur  from  the  last  Marx  spark  pin  to  the  cathode, 
effectively  bypassing  the  pulse  stretching  inductor.  This  problem  was 
somewhat  alleviated  by  distributing  the  inductance  throughout  the  Marx 
circuit  as  shown  in  3.12b. 

However  when  the  inductance  got  large  the  problem  returned  with 
the  arcing  being  to  some  other  Marx  pin.  This  problem  is  serious  because 
arcing  from  a  Marx  pin  to  the  cathode  produces  a  carbon  track  in  the 
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Fig.  3.13  Voltage  and  Current  Waveforms  of  Stretched  Pulse. 
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perspex  insulation  that  acts  as  a  short  circuit  and  the  system  is 
disabled  until  the  carbon  track  is  cleaned  off.  The  only  solution 
to  this  problem  is  to  replace  the  perspex  with  a  more  damage  resistant 
material,  such  as  ceramic,  and  to  move  the  spark  pins  further  away  from 
the  electrodes,  which  has  the  undesirable  effect  of  reducing  the  intensity 
of  the  preionizing  ultraviolet  radiation. 


3.5  ELECTRODES 

Also  considered  with  the  discharge  circuitry  are  the  electrodes. 
Three  electrode  materials  were  tried:  carbon,  aluminum  and  copper.  The 
carbon  electrodes  are  very  resistant  to  arcing  damage  but  have  two  draw¬ 
backs.  One  is  that  a  fine  carbon  dust  is  formed  which  coats  the  inside 
of  the  laser  system  (including  optics) .  The  other  problem  is  the  grainy 
structure  of  the  carbon  promotes  the  formation  of  many  small  cathode 
filaments  which  are  potential  arc  sites.  Aluminum  can  be  polished  to  a 
very  smooth  surface  which  is  less  prone  to  arc  but  it  is  easily  damaged 
when  arcing  does  occur  and  a  large  arc  will  burn  a  pit  into  the  surface. 
Then  to  eliminate  the  arc  it  is  necessary  to  remove  and  repolish  the 
electrode.  Copper  seems  to  be  the  best  choice  of  the  three.  It  can  be 
polished  to  a  smooth  surface,  is  less  prone  to  arc  than  carbon,  and  is 
more  resistant  to  arcing  damage  than  aluminum. 

Copper  was  finally  chosen  as  the  electrode  material  for  this  system. 
The  electrodes  for  this  experiment  were  formed  by  plating  copper  over  an 
aluminum  electrode. 

3.6  REPETITION  RATE 

One  goal  of  this  project  was  to  see  how  high  a  repetition  rate 
could  be  attained  with  this  discharge  configuration.  Although  repetition 
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rates  achieved  were  high  by  ordinary  TEA  laser  standards  they  were 
lower  than  expected  from  the  clearing  time  of  gas  between  the  electrodes. 
The  discharge  dimension  in  the  direction  of  flow  is  5  cm.  and  the  flow 
velocity  is  100  m/sec.  So  this  region  should  be  swept  clear  in  0.5  msec, 
implying  a  maximum  repetition  rate  of  2000  Hz.  Although  this  was 
achieved  in  He,  observed  repetition  rates  were  considerably  lower  in 
the  laser  mixtures.  Typical  maximum  repetition  rates  are  given  in  table 
3.2.  In  the  laser  gas  mixtures  the  maximum  repetition  rate  was  limited 
by  discharge  arcing  but  in  He  the  maximum  repetition  rate  was  limited 
by  the  Marx  generator  spark  gaps  going  into  continuous  conduction.  Since 
the  spark  gaps  are  5  cm.  downstream  from  the  discharge  the  2000  Hz 
repetition  rate  corresponds  to  the  time  it  takes  ionization  to  blow  from 
the  discharge  region  to  the  spark  pins,  reducing  their  breakdown  voltage, 
and  causing  them  to  go  into  continuous  conduction. 

A  related  problem  which  occurred  at  high  repetition  rates  and  for 
pulse  stretching  experiments  was  arcing  from  one  of  the  electrodes  to 
the  Marx  pins.  This  arcing  took  the  form  of  surface  tracking  along  the 
perspex.  When  it  occurred  the  system  had  to  be  dismantled  and  the  arc 
cleaned  off  because  it  would  short  out  the  circuit. 


Marx  Stages 

Energy  per  Pulse 

J 

Gas  Mixture 

Maximum  Repetition 
Rate 

Hz 

10 

20 

Helium 

700  torr 

2000 

7 

3 

1:1:8 

COg :^2 : He 

700  torr 

500 

10 

54 

70:60:600 

CO2 :^2 :He 

730  torr 

200 

Table  3.2 
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It  was  also  observed  that  after  several  minutes  of  running,  a 
stable  high  repetition  rate  discharge  would  begin  to  arc  sporadically. 
This  was  a  good  indication  that  gas  poisoning  was  occurring. 

In  no  case  was  it  possible  to  reliably  maintain  a  high  repetition 
rate  for  long  periods.  For  instance,  a  typical  run  in  a  700  torr  laser 
mix  for  a  10  stage  Marx  with  .3  yf  per  stage  was  as  follows.  286  pps 
was  maintained  for*  40  seconds  arc  free  after  the  initial  fill.  Then 
sporadic  arcing  began  in  the  discharge.  After  3  more  minutes  the  run  had 
to  be  terminated  because  an  electrical  feedthrough  in  the  perspex  got 
too  hot  causing  the  perspex  to  melt. 

One  of  the  factors  limiting  attainable  repetition  rate  seemed  to 
be  Marx  charging  time,  which  was  limited  by  the  LC  time  constant  of  the 
circuit  to  typically  3.5  msec,  for  a  total  Marx  capacitance  of  3  yf.  L 
in  this  case  is  the  power  supply  inductance  and  from  the  time  constant 
is  calculated  to  be  approximately  1.7  H.  An  attempt  to  pulse  faster 
than  this  moves  the  operating  region  of  the  Marx  from  the  plateau  A 
to  the  rising  voltage  B  on  the  charging  waveform  (Fig.  3.14).  On  the 
rising  voltage  portion  of  the  waveform  it  is  difficult  to  get  the  Marx  to 
operate  reliably. 

3.7  OPTICS 

Three  different  optical  resonators  were  tested  for  this  system: 
an  unstable  resonator,  a  stable  resonator  with  a  germanium  output  mirror, 
and  a  stable  resonator  with  a  hole  coupled  output  mirror. 

Originally  an  unstable  resonator,  shown  in  Fig.  3.15,  was  designed 
using  the  procedure  outlined  by  Siegman  (32)  .  It  had  a  cavity  length 
of  1.2  m  and  an  output  coupling  of  10%.  The  mirrors  were  gold  coated 
beryllium  copper,  one  being  concave  with  a  14.4  m  radius  of  curvature 
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Fig.  3.14  Typical  Marx  Generator  Charging  Waveform. 
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Figure  3.15  Optical  System  -  (a)  Unstable  resonator  and 

(b)  Stable  resonator. 
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and  the  other  convex  with  a  12  m  radius.  The  output  was  coupled  out 
through  a  salt  window  by  means  of  a  skimmer  plate.  A  typical  discharge 
current  pulse  and  laser  output  pulse  shape  are  shown  in  Fig.  3.16. 

The  laser  output  pulse  was  observed  using  a  gold  doped  germanium  detector. 

Two  skimmer  plates  were  tried,  one  with  a  round  hole  and  one  with  a 
rectangular  hole.  The  one  with  the  rectangular  output  worked  much  better 
because  it  used  more  of  the  cavity  volume.  It  had  an  efficiency  of  6% 
compared  to  2%  for  the  round  skimmer.  Efficiency  is  a  convenient  means 
of  comparing  the  performance  of  the  various  optical  systems.  It  is 
simply  optical  energy  out  divided  by  electrical  energy  in  per  pulse, 
converted  to  percent. 

Bellows  were  used  to  attach  the  unstable  resonator  directly  to 
the  discharge  cavity.  However  vibration  transmitted  by  the  bellows  made 
it  difficult  to  maintain  alignment  of  the  mirrors.  To  isolate  the  optics 
from  the  flow  system  the  bellows  were  removed  and  NaCl  windows  were 
attached  to  the  flow  system.  This  resulted  in  increased  output  stability 
but  the  losses  introduced  by  the  windows  in  the  laser  cavity  reduced  output 
significantly,  dropping  the  efficiency  for  the  round  skimmer  from  2%  to  0.8%, 
so  bellows  were  reinstalled. 

Since  the  unstable  resonator  was  very  sensitive  to  vibration  a 
stable  resonator  was  tested  for  comparison.  The  stable  resonator  was 
constructed  using  the  14.4  m  concave  mirror  and  a  hole  coupled  output 
mirror.  The  hole  coupled  output  mirror  is  simply  a  flat  mirror  with  an 
array  of  small  holes  through  it.  The  fraction  of  surface  area  occupied 
by  the  holes  is  equal  to  the  fractional  output  coupling  required.  A  hole 
coupled  mirror  was  used  because  it  can  withstand  high  repetition  rates 
without  being  damaged.  Unfortunately  the  hole  coupled  mirror  worked  very 
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0.5  p  s  e  c/d  i  v 


Fig.  3.16  Optical  Output.  The  upper  trace  is  the  discharge  current 

waveform  and  the  lower  trace  is  the  optical  output  pulse  shape 
observed  with  a  gold  doped  germanium  detector.  The  approximate 
energy  of  the  optical  pulse  is  0.4  J. 
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poorly  with  an  efficiency  of  typically  1%  or  less.  For  comparison, 
replacing  the  hole  coupled  mirror  with  a  Germanium  flat,  which  has  30% 
coupling,  raised  the  efficiency  to  3%.  But  the  Germanium  flat  can  be 
run  on  a  single  shot  basis  only. 

The  various  optical  systems  are  compared  in  table  3.3.  The  energy 
per  pulse  for  the  optical  systems  in  the  table  was  measured  on  a  single 
shot  basis  using  a  Gen  Tec  Joule  meter.  The  results  show  that  the 
unstable  resonator  with  the  10%  output  coupled  square  skimmer  was  the 
best  system  of  those  tested. 

A  gas  cell  beam  detector  described  by  Offenberger,  Smy  and  Burnett  (33) 
permitted  examination  of  the  energy  output  at  high  repetition  rates.  The 
cell  is  filled  with  a  mixture  of  He  at  atmospheric  pressure  and  a  few 
torr  of  an  infrared  absorbing  gas  such  as  propylene.  When  the  laser 
beam  passes  through  the  cell  some  of  the  energy  is  absorbed  and  produces 
a  shock  wave.  The  shock  wave  strength  is  sensed  by  a  piezoelectric 
probe  and  is  proportional  to  the  energy  in  the  laser  pulse.  This  probe 
recovers  in  typically  4  msec  compared  to  .3  sec  for  the  Joulemeter. 

Using  this  probe  the  output  of  the  laser  could  be  monitored  at  high 
repetition  rates.  Fig.  3.17  shows  a  typical  output  from  the  unstable 
resonator  at  100  Hz.  The  effect  of  the  slow  scale  mechanical  vibration 
can  be  seen  quite  clearly  and  is  responsible  for  a  variation  of  100%  in 
the  output  intensity  of  the  pulses. 

The  first  set  of  bellows  used  on  this  system  consisted  of  tubes 
cut  from  rubber  gloves  and  sliding  gate  valves  to  seal  the  bellows  off 
from  the  system  during  pump-down  and  gas  changes,  because  they  could  not 
withstand  vacuum.  Although  this  system  worked  it  was  unsatisfactory 
because  the  rubber  failed  frequently.  The  main  cause  of  failure  was 
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Type  of  Resonator 


Input  Energy 
J 


Output  Energy 
J 


Efficiency 

% 


unstable 

square  skimmer  16 

10% 

unstable 

round  skimmer  19 

10% 

unstable 

round  skimmer  19 

intercavity  windows 

stable 

germanium  flat  45 

30% 

stable 

hole  coupled  25 

5% 

stable 

hole  coupled  54 

15% 

stable 

hole  coupled  54 

20% 

germanium  flat 

low  pressure  discharge  13 


.4 


2 


.16 


.8 


1.5 


3 


.3 


1 


.5 


.9 


.25 


.5 


.37 


3 


Table  3.3 

Comparison  of  performance  of  the  optical  systems  used  on  the  high 


repetition  rate  laser. 
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Fig.  3.17  Optical  Output  Pulses  From  the  Unstable  Resonator  at  100Hz 


Repetition  Rate. 
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deterioration  of  the  rubber  due  to  contact  with  ozone  produced  by  the 
discharge. But  a  misaligned  laser  beam  grazing  the  surface  of  the  rubber 
and  burning  a  hole  was  also  responsible  for  some  failures. 

Because  this  bellows  failed  so  frequently  it  was  replaced  by  a 
bellows  suggested  by  Dr.  H.J.J.  Seguin  which  proved  to  be  much  more  reliable. 
The  new  bellows  consisted  of  metal  tubes  joined  by  rubber  tube. 

This  system  was  better  because  it  eliminated  the  problem  of  the  rubber 
failing  and  also  could  withstand  vacuum  so  the  gate  valves  were  no  longer 
needed.  Both  bellows  systems  are  shown  in  Fig.  3.18. 


3.8  GAIN  MEASUREMENTS 

Small-signal  gain  of  this  laser  was  measured  using  the  apparatus 
shown  in  Fig.  3.19a.  The  output  from  a  small  continuous  wave  CO^  probe 
laser  was  chopped  to  produce  a  square  wave.  The  resulting  pulse  train 
then  made  a  double  pass  through  the  discharge  region  and  was  directed 
into  a  detector.  A  gold  doped  germanium  detector  was  used  with  a  Ge  filter 
in  front  of  it.  The  filter  eliminated  the  light  flash  from  the  discharge, 
limiting  the  light  detected  to  a  narrow  band  around  10.6  ym.  When  the 
discharge  fired  a  pulse  was  produced  on  top  of  the  square  wave  due  to 
amplification  in  the  discharge  region.  The  gain  can  be  calculated  from 
the  ratio  of  the  peak  intensity  of  the  gain  pulse  to  the  intensity  of  the 
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Fig.  3.18  Mirror  Bellows  Systems,  (a)  is  the  rubber  tube  bellows  with  a 
sliding  gate  valve  and  (b)  is  the  metal  tube  bellows. 
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Cb) 

Fig.  3.19  (a)  Shows  the  Apparatus  used  for  the  Small-signal  Gain  Measurements 
and  (b)  Shows  a  Typical  Gain  Pulse  Waveform. 
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where  Iq  is  the  intensity  of  the  square  wave,  I  is  the  intensity  of  the 
gain  pulse,  l  is  the  path  length  in  the  active  medium,  and  a  is  the  small 
signal  gain.  A  typical  small  signal  gain  pulse  is  shown  in  Fig.  3.19b. 

A  double  pass  of  the  probe  beam  through  the  discharge  served  two 
purposes.  It  increased  the  size  of  the  gain  signal  and  made  it  possible 
to  move  the  detector  away  from  the  laser  which  produced  a  lot  of 
electrical  noise. 

Measured  small  signal  gain  was  0.5%  to  0.7%  per  cm.  Measurements 
were  made  at  repetition  rates  of  up  to  200  Hz  and  gain  was  independent 
of  repetition  rate. 

3.9  SHOCK  WAVE  MEASUREMENTS 

During  the  gain  measurements  of  this  laser  an  interesting  anomaly 
accompanied  many  of  the  gain  pulses.  It  was  a  negative  swing  at  the  end 
of  the  gain  pulse.  The  negative  swing  was  due  to  discharge  induced  gas 
density  gradients  deflecting  the  probe  beam  so  that  it  partially  missed 
the  detector.  With  the  blower  running  the  density  gradients  dissipated 
in  less  than  1  msec  due  to  convection  but  if  the  blower  was  not  running 
the  gradients  had  to  dissipate  by  diffusion  and  took  about  300  msec  to 
disappear  (Fig.  3.20). 

An  experiment  was  devised  to  study  these  density  perturbations 
in  more  detail  since  if  they  persist  long  enough  they  could  degrade  the 
optical  quality,  and  possibly  even  the  discharge  uniformity,  of  subsequent 
pulses . 

A  shock  wave  in  a  gas  produces  an  index  of  refraction  gradient  that 
deflects  light  rays  passing  through  the  gas  according  to  (ref.  34) 
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Fig.  3.20  Gas  Density  Perturbations  (a)  Without  Blower  Running  and 

(b)  With  Blower  Running.  The  upper  trace  in  (a)  is  a  discharge 


current  waveform. 
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3.9.1 


where  k  is  the  Gladstone  Dale  constant,  t  is  the  thickness  (or  distance 
through  the  gas  the  beam  travels),  is  the  ambient  gas  density,  n  is 
the  ratio  of  density  at  any  instant  to  the  original  ambient  density  and 
0  is  the  angular  deflection  of  the  beam.  So,  as  the  gain  measurements 
already  indicated,  passing  a  laser  beam  through  the  discharge  region 
and  monitoring  its  deflection  provides  a  sensitive  means  of  monitoring 
discharge  produced  density  gradients.  It  was  found  that  for  this  experi¬ 
ment  an  angular  displacement  of  .3  mrad  could  easily  be  detected  for  a 

discharge  in  He  gas.  Using  3.9.1  this  indicates  that  a  density  gradient 

-2 

of  1.4x10  /mm  or  1.4%  per  mm  could  be  detected. 

Measurements  were  made  for  a  He  discharge  with  the  apparatus  shown 
in  Fig.  3.21.  A  He-Ne  laser  beam,  after  crossing  the  discharge  region 
is  directed  past  a  knife  edge,  through  a  red  filter,  and  then  to  a  photo¬ 
voltaic  detector.  The  red  filter  blocks  the  broadband  light  flash 
produced  by  the  discharge  and  only  allows  light  from  the  He-Ne  to  reach 
the  detector.  The  knife  edge  intercepts  half  the  probe  beam  and  a 
deflection  of  the  beam  toward  the  edge  will  reduce  the  amount  of  light 
striking  the  detector.  Similarly  a  deflection  away  from  the  edge  will 
increase  the  amount  of  light  striking  the  detector. 

Changing  the  orientation  of  the  knife  edge  makes  it  possible  to 
monitor  different  components  of  the  shock  wave.  Fig.  3.22  is  a  cross 
sectional  view  of  the  discharge  region  showing  the  points  at  which  shock 
wave  measurements  were  made.  Orienting  the  knife  edge  vertically  allows 
monitoring  of  shock  waves  travelling  in  the  x  direction  (horizontal) , 
and  orienting  it  horizontally  allows  monitoring  shock  waves  travelling 
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Fig.  3.21  Apparatus  for  Density  Gradient  Measurements 
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3.22 


Cross  Sectional  View  of  Discharge  Region  Showing  Points 
at  Which  Shock  Wave  Measurements  Were  Made. 
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in  the  y  direction  (vertical) .  By  moving  mirror  1  in  the  horizontal 
plane  and  the  He-Ne  probe  laser  in  the  vertical  plane  the  shock  waves 
could  be  mapped  at  various  positions  in  the  discharge  region. 

The  knife  edge  and  detector  were  28  ft  away  from  the  discharge 
chamber.  This  large  displacement  allowed  resolution  of  small  angular 
deflections  on  the  scope.  A  calibration  figure  was  obtained  by  moving 
the  knife  edge  across  the  beam  with  a  micrometer.  One  scope  division, 
at  . 2mV  per  division,  corresponded  to  a  deflection  of  .090  inches,  or 
a  sensitivity  of  0.3  mrad/div.  Since  the  beam  is  gaussian  and  has  a 
circular  cross  section  the  actual  sensitivity  becomes  non-linear  for  large 
deflections . 

Shock  waves  were  observed  in  both  static  and  flowing  gas.  A  large 
slowly  varying  density  gradient  was  observed  in  static  gas  which  was 
absent  in  the  flowing  gas  (Fig.  3.23).  The  large  scale  deflection  typically 
has  a  decay  time  constant  of  150  msec,  and  is  due  to  the  slow  diffusion  of 
the  heated  gas  remaining  after  the  discharge.  In  contrast,  with  gas  flow, 
the  heated  gas  is  quickly  removed  from  the  discharge  region  by  convection 
and  the  thermal  gradient  does  not  appear.  The  fluctuations  which  occur 
in  flowing  gas  are  due  to  shock  waves  and  mostly  damp  out  within  0.5msec. 

The  peak  values  of  thermal  deflection  are  given  in  Table  3.4.  It 
can  be  seen  that  the  largest  component  of  these  deflections  was  in  the 
vertical  direction.  A  plot  of  the  peak  values  along  the  central  axis 
(points  5,  4,  7,  9)  is  shown  in  Fig.  3.24a.  This  suggests  a  density 
gradient  and  hence  density  distribution  of  the  form  shown  in  Fig.  3.24b. 
Since  the  gas  heating  is  proportional  to  the  current  density  through  it, 
the  data  offers  some  insight  into  the  uniformity  of  the  discharge.  It 
suggests  that  the  energy  was  deposited  with  a  gaussian  profile  rather 
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Fig.  3.23  Detector  Voltage  (a)  Without  Gas  Flow  and  (b)  With  Gas  Flow. 
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Peak 

Deflection 


Position  # 

X 

inches 

y 

inches 

vertical  y 
division 

horizontal  x 
division 

1 

5/8 

1 

-1.5 

-.25 

2 

6/8 

1 

-1.0 

-.2 

3 

4/8 

1 

-1.5 

-.4 

4 

1 

1 

-1.5 

-0 

5 

1 

1/2 

-.4 

-0 

6 

1/2 

1/2 

-.1 

-.4 

7 

1 

1  1/2 

-4.0 

-.2 

8 

1/2 

1  1/2 

-4.8 

1.2 

9 

1 

1  7/8 

+2 

-.3 

Table  3.4 

Probe  Beam  Deflection  Due  to  Thermal  Density  Gradients. 
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(a) 


(b) 


Fig.  3.24  (a)  Deflection  of  Probe  Beam  Due  to  the  Thermal  Density  Gradient 
and  (b)  Model  Density  Gradient  dn/dx  with  Model  Density 


Distribution  n. 
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than  uniformly  through  the  discharge  region.  It  also  shows  that  the 
discharge  was  asymetrical,  with  most  of  the  energy  being  deposited  near 
the  top  edge.  This  was  verified  by  the  visual  observation  that  the 
discharge  was  brightest  along  the  top  edge  of  the  electrodes.  The  non¬ 
uniformity  is  most  likely  due  to  the  unsymmetrical  nature  of  the  pre¬ 
ionization,  since  the  spark  pins  are  on  one  side  only,  and  represents  a 
shortcoming  of  the  integrated  switching  and  preionization  technique. 

Kline  and  Denes  (35)  point  out  the  necessity  of  transversely  uniform  pre¬ 
ionization  in  their  study  of  glow  discharges. 

Fig.  3.25  shows  the  vertical  and  horizontal  components  of  shock  waves 
at  point  4.  In  the  vertical  direction  a  large  refractive  index  gradient 
is  followed  by  an  oscillation.  The  first  large  peak  is  due  to  gas  heating 
and  expanding  in  the  vertical  direction.  The  time  lag  between  the  first 
pulse  and  the  oscillations  which  follow  is  consistent  with  the  oscillations 
coming  from  the  spark  pins.  Typically  the  oscillations  took  104  psec  to 
arrive  in  static  gas  and  112  psec  to  arrive  in  flowing  gas. 

In  the  horizontal  direction  the  dominant  feature  of  the  shock  wave 
is  two  spikes  which  correspond  to  the  time  for  a  shock  wave  to  travel 
from  the  measuring  point  to  the  anode  and  be  reflected  back  again.  The 

time  between  spikes  for  positions  1,  2,  3,  4  is  shown  in  Table  3.5  and 

-2 

is  compared  with  transit  times  calculated  using  v=9. 65x10  “  cm/psec,  the 
velocity  of  sound  in  He.  The  data  is  consistent  with  a  shock  wave  travelling 

from  the  cathode  to  the  anode  and  then  being  reflected  but  the  velocity 

_2 

is  lower  than  expected,  having  an  average  value  of  only  7.3x10  cm/psec. 

The  likely  cause  of  this  lower  velocity  is  leaking  into  the  discharge 
system.  The  shock  wave  is  produced  at  the  cathode,  due  to  gas  heating 
in  the  cathode  fall  region.  This  result  agrees  with  cathode  shock  waves 
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0.2  m  v/d  i  v 


2  0  fJ  se  c/div 


(b) 

Fig.  3.25  Detector  Voltage  at  Postion  4  Indicating  (a)  Vertical  and 
(b)  Horizontal  Components  of  Shock  Waves. 
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;ition  # 

Distance  from  P 
to  Anode  and 

Back  to  P 

cm 

t  calc 

gsec 

t  measured 

psec 

V  measured 
cm/ psec 

4 

3.5 

36 

44 

7 . 9xl0~2 

2 

2.2 

23 

30 

7.3xl0"2 

1 

1 . 6 

16 

20 

8 . OxlO_2 

3 

1.0 

oo 

16 

6.2xl0~2 

Table  3.5 

Cathode  Shock  Wave  Propagation  Measurements. 
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observed  by  Clark  (36)  in  a  CO^  laser  discharge  using  a  holographic 
recording  technique.  Clarke  also  observed  nonuniform  gas  heating  effects 
due  to  asymmetrical  preionization. 

Most  of  the  shock  wave  measurements  are  for  static  gas. 
Measurements  made  in  flowing  gas  for  comparison  showed  that  gas  flow 
had  little  effect  on  the  shock  waves,  other  than  reducing  their 
amplitude  slightly  and  increasing  the  time  taken  for  spark  pin 
oscillations  to  reach  the  discharge  region.  In  static  gas  large 
thermal  gradients  were  present  that  took  about  0.5  sec  to  disperse  but 
in  flowing  gas  they  disappeared.  The  shock  waves  and  spark  pin 
oscillations  had  mostly  dispersed  within  0.5  msec  of  the  discharge 
in  both  flowing  and  static  gas.  So  shock  wave  induced  gas  density 
nonuniformities  should  not  cause  discharge  problems  in  this  system 
for  repetition  rates  under  2  KHz. 


3.10  DISCUSSION  AND  CONCLUSIONS 

This  system  has  shown  the  potential  of  an  integrated  switching 
and  preionization  technique  to  produce  a  simple  high  repetition  rate 
TEA  laser.  Typical  repetition  rates  of  200  Hz  were  achieved  and  the  dis¬ 
charge  in  He  suggests  that  repetition  rates  of  2  KHz  should  be  possible 
before  failure  of  the  Marx  switching  sets  in.  Pulse  stretching  was 
attempted  but  was  not  successful  and  short  pulses  of  about  1  psec 
duration  seem  to  be  best.  Maximum  output  energy  per  pulse  was  1  J  giving 
a  typical  average  output  power  of  200  W  with  an  electr ical-to-light 
conversion  efficiency  of  6%. 
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Although  the  operation  of  this  system  was  encouraging  it  still 
had  a  major  shortcoming,  the  inability  to  run  for  an  extended  period 
of  time.  Runs  were  short  and  were  usually  terminated  by  the  onset  of 
arcing  problems.  Arcing  could  occur  between  the  electrodes  or  from  an 
electrode  to  the  Marx  spark  pins.  Arcing  from  an  electrode  to  the  Marx 
pins  was  particularly  destructive  because  the  resulting  carbon  track 
on  the  perspex  insulation  had  to  be  cleaned  off  before  the  system  could 
run  again  which  meant  the  system  had  to  be  shut  down  and  dissassembled . 

Arcing  between  the  electrodes  was  less  destructive  but  still  meant  term¬ 
inating  the  run  to  avoid  damaging  the  electrodes. 

The  cause  of  the  arcing  was  unclear  but  several  possibilities 
were  considered.  One  possibility,  that  of  shock  waves  adversely  affecting 
the  discharge,  was  ruled  out  by  the  measurements  in  section  3.9  which 
showed  that  they  did  not  persist  long  enough  in  this  system  to  cause 
problems.  Other  possibilities  considered  were  inadequate  preionization, 
some  kind  of  Marx  generator  malfunction,  air  leaking  into  the  system  and 
contaminating  it,  or  gas  poisoning  due  to  discharge  products.  Some 
observations  on  this  system  and  by  other  researchers  suggested  that  gas 
poisoning  due  to  dissociation  of  the  laser  gases  was  a  possible  cause. 

But  the  vacuum  integrity  of  this  system  was  not  very  good  and  the  possibility 
of  air  leaks  degrading  the  performance  also  existed.  Also  it  was  not 
clear  whether  or  not  the  Marx  generator  was  responsible  in  some  way  for 
the  arcing.  For  instance,  the  preionization  it  provided  was  not  entirely 
uniform  transverse  to  the  discharge.  The  remainder  of  this  thesis  deals 
with  experiments  aimed  at  determining  the  cause  of  arcing. 
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CHAPTER  4 


FLOWING  GAS  ELECTRICAL  DISCHARGE  SYSTEM  DESIGN 
INTRODUCTION 

The  purpose  of  the  flowing  gas  electrical  discharge  system  is  to 
provide  a  test  facility  to  study  the  causes  of  arcing  in  a  high  repetition 
rate  glow  discharge  since  arcing  problems  had  proved  to  be  a  major 
hinderance  in  the  laser  system.  This  discharge  system  is  modelled  after 
the  one  designed  by  Olson  et  al.  (21)  and  is  essentially  a  simpler  version 
of  theirs. 

It  has  several  improvements  over  the  laser  system,  one  of  the  most 
important  being  its  smaller  size.  Because  it  is  smaller  only  one  person 
is  needed  to  dissassemble  it  and  modifications  are  easier  to  make.  The 
vacuum  integrity  of  the  discharge  system  is  much  better  than  that  of  the 
laser,  thereby  ruling  out  air  contamination  as  a  source  of  problems.  In 
order  to  model  the  laser  system,  the  discharge  was  formed  between  two 
copper  electrodes  with  preionization  provided  by  brass  spark  pins.  Rather 
than  using  a  Marx  generator,  switching  is  accomplished  with  a  thyratron 
pulser.  This  eliminated  problems  that  might  have  been  due  to  running  a 
Marx  generator  in  the  gas  flow  loop  and  also  results  in  cleaner  voltage 
waveforms  than  were  achieved  with  spark  gap  switching.  It  also  eliminated 
the  problem  of  arcing  putting  the  switching  circuit  out  of  commission 
which  happened  frequently  with  the  Marx. 

The  system  consists  of  two  main  subsections.  One  is  the  discharge 
circuitry  and  the  other  is  the  recirculating  gas  loop  which  contains 
the  electrodes,  a  heat  exchanger,  and  a  fan.  Driving  the  fan  with  a 
variable  speed  motor  allows  control  of  the  gas  flow  velocity  in  the  loop. 
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4.1  RECIRCULATING  GAS  LOOP 

The  recirculating  gas  loop  is  shown  in  Fig.  4.1.  It  is  constructed 
from  6  inch  diameter  aluminum  pipe  with  1/8"  wall  thickness.  The  fan 
housing  and  elbow  which  are  12"  in  diameter  were  constructed  by  rolling 
and  welding  aluminum  plate.  The  heat  exchanger  is  housed  in  a  box 
constructed  from  aluminum  channel.  All  sections  are  joined  together 
using  o-ring  sealed  flanges. 


4.2  PRESSURE  VESSEL  CAPABILITY 

The  pressure  that  the  system  could  be  expected  to  withstand  was 
estimated  using  a  formula  from  the  Mechanical  Engineer's  Reference  Book 
(37).  The  collapsing  pressure  for  a  thin  walled  pipe  is  given  by 


P 

cr 


2  t3E 


3  2 

d  (1-v  ) 


4.2.1 


where  t  is  wall  thickness,  d  is  diameter,  E  is  modulus  of  elasticity, 

v  is  Poisson's  ratio  which  lies  between  0  and  0.5  (between  .25  and  .36 

2 

for  metals).  For  aluminum  alloys  E  is  about  60,000  N/mm  so  P  lies  in 

cr 

2 

the  range  171  to  184  lb/in  for  6"  diameter  aluminum  pipe.  Since 

2 

atmospheric  pressure  is  14.7  lb/in  the  system  will  easily  withstand  being 
pumped  out  without  collapsing  with  a  wide  safety  margin. 


4.3  FLOW  SYSTEM  DESIGN 

In  order  to  know  what  fan  capacity  would  be  needed  to  achieve  a  good 
gas  flow  velocity  it  was  necessary  to  estimate  the  friction  loss  of  the 
piping  in  the  gas  flow  loop.  The  methods  outlined  in  the  ASHRAE  handbook 
(27)  were  used  and  are  essentially  based  on  the  formula 


f  V  (4005) 
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4.3.1 
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where  f  is  a  friction  coefficient  (-  2.5x10  ) ,  £  is  length  of  the  pipe 

in  ft.,  D  is  the  diameter  of  the  pipe  in  ft.,  v  is  gas  flow  velocity  in 
ft/min,  and  is  pressure  head  loss  in  inches  of  H^O.  Noncircular  pipe 
such  as  the  test  section  is  converted  into  an  equivalent  circular  duct 
using  the  formula 


d 

c 


1.3 


(ab) 


.625 


(a+b) 


.250 


4.3.2 


where  d  is  the  equivalent  circular  diameter  of  a  duct  with  sides  of 
c 

length  a  and  b.  Elbows  were  converted  into  an  equivalent  length  of 
straight  pipe  using  a  chart  from  the  ASHRAE  handbook, and  assuming  3 
section  elbows  with  6  inch  diameter  and  1  ft  turning  radius  led  to  an 
equivalent  length  of  8  ft  of  straight  pipe  for  each  elbow.  Fortunately 
smooth  elbows  were  obtained  which  are  better  than  the  3  section  elbows 
assumed  in  this  estimate.  Their  effective  length  is  only  5  feet  each. 
The  loss  due  to  changes  in  area  such  as  diffusers  or  converging  sections 
is  given  by 


H  =  C 


,v 


4005' 


4.3.3 


where  v  is  the  velocity  at  the  inlet  to  the  transition  section  and  C 
is  a  constant  which  is  obtained  from  the  ASHRAE  handbook,  C  is  different 
for  different  transition  sections,  depending  on  the  angle  between  the 
diffuser  or  converger  walls. 

Initial  pressure  head  loss  estimates  were  calculated  using  the 
simplified  system  model  shown  in  Fig.  4.2  The  head  loss  due  to  the  heat 
exchanger  was  ignored  at  this  stage  of  the  calculations  because  it  was 
assumed  to  be  small  and  there  was  insufficient  information  to  estimate  it. 
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Estimates  of  Flow  Velocity. 
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The  continuity  equation  requires  that  volumetric  flow  be  the  same 
everywhere  in  the  system  so  the  losses  in  each  section  were  converted  to 
a  function  of  volume  flow  instead  of  velocity  by  using 

CFM  =  V  x  cross  sectional  area 
D  2 

=  V  7T ('2')  (for  round  pipe) 

The  pressure  loss  through  the  system  in  inches  of  H^O  for  various 
gas  flow  rates  in  CFM  is  shown  in  Table  4.1. 

The  cross  sectional  area  of  the  working  section  is  l"x6"  so  the 
flow  velocity  between  the  electrodes  can  be  calculated  from  the  volume 
flow  using  the  continuity  equation.  It  can  be  seen  that  the  diffuser 
is  the  primary  source  of  system  head  loss,  accounting  for  half  the  loss 
by  itself. 

To  circulate  the  gas  a  fan  was  chosen  from  Sheldon’s  Multi-wing  series. 
This  series  of  fans  has  removable  plastic  blades  making  it  possible  to 
customize  a  wide  variety  of  sizes  and  blade  angles.  They  are  also  fairly 
cheap  and  readily  available.  The  final  choice  was  a  10  bladed  fan  with  a 
45°  angle  blade  pitch  and  12"  diameter.  From  fan  rating  tables  published 
by  Sheldon's  typical  expected  gas  flow  and  power  requirements  for  a  given 
fan  speed  can  be  estimated.  Fig.  4.3  shows  the  initial  estimates  of  gas 
flow  in  CFM  versus  fan  rpm.  This  however  was  for  a  6  bladed  12  inch 
diameter  fan  with  35°  blade  pitch.  In  addition  the  curve  for  the  10 
bladed  45°  pitch  fan  is  shown  for  comparison. 

Where  the  family  of  fan  characteristics  intersects  the  system  loss 
curve  defines  the  gas  flow  delivered  in  the  system  and  the  amount  of 
power  needed  to  deliver  that  gas  flow.  Unfortunately  it  was  necessary 
to  extrapolate  considerably  beyond  the  values  given  in  the  rating  tables 
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Item 

Pressure 

"h2o 

loss 

12"  dia.  pipe 

7.6xl0"9 

(CFM)2 

6"  dia.  pipe 

3.2xl0_6 

(CFM) 2 

l"x6"  cross  section 
working  section 

2.2xl0"6 

(CFM)2 

converger 

2.5xl0"6 

(CFM)2 

— 6 

2 

diffuser 

9.9x10 

(CFM) 

transition  6"  to 

12"  pipe 

1.6xl0"6 

(CFM)2 

transition  12"  to 

6"  pipe 

5.2xl0"7 

(CFM)2 

-5 

s  2 

TOTAL 

2.0x10 

(CFM) 

Table  4.1 

Pressure  Head  Loss  Estimates 


Throughout  the  Flowing  Gas  Discharge  System 
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Fig.  4.3  Fan  Characteristics  and  System  Loss  Curve. 

intersection  of  the  characteristics  with  the  loss  curve 

determines  the  gas  flow  in  the  system  for  a  given  fan  speed. 

The  fan  characteristics  come  from  tables  to  the  left  of  the 
vertical  bars  and  are  extrapolation  to  the  right. 
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so  the  results  will  not  be  very  accurate.  However  they  will  give  some 
idea  of  what  flow  velocity  to  expect. 

After  looking  for  a  motor  to  drive  the  fan  it  became  apparent  that 
the  size  of  the  motor  rather  than  fan  capacity  would  determine  maximum 
gas  flow  velocity.  For  instance,  to  match  the  flow  velocity  of  100  m/sec 
of  the  laser  system  in  air  would  require  about  6  HP.  Besides  the  fact 
that  this  power  input  would  require  a  very  large  heat  exchanger  (since  it 
would  have  to  remove  6  HP  of  heat  flow  from  the  gas)  such  a  large  variable 
speed  motor  is  not  readily  available.  The  simplest  way  to  vary  the  gas 
flow  was  by  varying  the  fan  speed  using  a  universal  motor  with  a  triac 
controller.  The  largest  motor  of  this  type  that  could  be  easily  obtained 
was  a  15A  industrial  disc  grinder.  This  unit  has  a  capability  of  only 
about  2  HP.  Higher  powered  motors  are  available  but  not  suitable.  For 
instance  higher  powered  induction  motors  are  easily  obtained  but  their 
speed  cannot  be  varied  except  by  going  to  complicated  and  costly  techniques 
such  as  powering  them  from  a  variable  frequency  source.  Large  D.C.  motors 
are  also  available  but  they  are  expensive  and  not  as  easy  to  obtain.  In 
addition  there  is  still  the  problem,  with  a  larger  motor,  of  being  able 
to  handle  the  heat  load  in  the  gas  flow  system. 

Using  the  2  HP  grinder  motor  an  initial  estimate  can  be  made  of 
obtainable  gas  flow.  A  2  HP  motor  should  be  able  to  deliver  approximately 
500  CFM  of  air  against  5.5"  H20  static  pressure.  This  would  result  in 
a  gas  flow  velocity  of  approximately  60  m/ sec  between  the  electrodes  in 
the  test  section. 

4.4  HEAT  EXCHANGER  DESIGN 

The  heat  exchanger  was  one  of  the  more  difficult  parts  of  the  system 
to  design,  mainly  because  of  the  lack  of  reliable  specifications  for  the 
finned  tubing.  First  an  estimate  of  the  heat  load  of  the  system  is  needed. 
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Assuming  1  KHz  operation  for  the  discharge  with  1  J/pulse  energy  leads 
to  a  heat  load  from  the  discharge  of  approximately  4000  BTU/hr.  In 
addition  there  is  of  course  the  heat  load  from  the  fan.  For  2  HP  this 
would  be  a  total  heat  load  of  approximately  10,000  BTU/hr.  The  heat  exchanger 
was  to  be  constructed  from  Aerofin  type  140  water  coils  which  is  a  finned 
tubing  with  14  fins  per  inch,  and  was  to  be  connected  in  a  crossflow 
arrangement . 


The  equation  which  describes  heat  flow  is 


q  =  MCAT 


For  gas  at  atmospheric  pr 


essure  and 


_  q _  _  10000  BTU/hr 

^ 1  Mf 

2100  lb/hr. 248  BTU/lb°F 


19 . 2°F 


So  for  a  gas  inlet  temperature  of  about  100°F  the  outlet  temperature  would 

be  81°F.  The  water  supply  for  the  heat  exchanger  is  just  tap  water.  The 

3 

tap  in  our  lab  delivers  .23  ft  /min  of  water  or  859  lb/hr  with  a  temperature 
of  59°F.  (In  fact  this  temperature  varies  with  the  time  of  year.)  For 
the  water 


=  10000  BTU/hr _ 

859  lb/hr  1  BTU/lb°F 


so  for  an  inlet  temperature  of  50 °F  the  water  will  have  an  outlet 
temperature  of  71°F.  This  gives  a  mean  temperature  difference  between  the 
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gas  and  water  of 


at  =  (100-59)  -  (81-71) 
m  „  ,100-59, 


22  °F 


The  face  of  the  heat  exchanger  is  12"  in  diameter  and  has  eight 

pipes.  The  I.D.  of  the  finned  tubes  is  9/16",  and  for  eight  pipes  in 

2 

the  face  represents  an  area  of  .014  ft  .  So  for  a  water  volume  of 
3 

.23  ft  /min  the  water  velocity  through  the  pipes  would  be  .28  ft/sec. 
From  this  information  and  from  aerofin  tables  a  heat  transfer  constant 
can  be  estimated. 


100  BTU 
hr  f t 2  °F  row 


This  estimate  is  very  crude  because  it  was  arrived  at  by  extrapolating 
from  the  charts. 


nrows  KAAT 

m 


10000  BTU/hr 

100  BTU/hr  ft2°F  x.79  ft2x22°F 

row 


5.8  rows 


where  A  is  the  face  area  of  the  heat  exchanger  and  n  is  the  number  of 
rows  required.  So  6  rows  of  pipes  are  needed  for  the  heat  exchanger 
with  8  pipes  in  parallel  in  each  row.  At  500  CFM  the  face  velocity  of 
the  gas  would  be  approximately  640  ft/min  and  the  pressure  loss  would  be 
.75"  H^O.  This  is  small  compared  to  the  other  system  losses  but  it  would 
probably  be  enough  to  drop  the  actual  gas  flow  down  to  about  400  CFM,  or 
a  flow  velocity  between  the  electrodes  of  about  50  m/sec.  Referring  back 
to  Fig.  4.3  the  system  loss  curve,  corrected  for  heat  exchanger  losses  is 


shown  on  the  graph. 


- 
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The  main  point  of  this  heat  exchanger  design  was  to  try  and  insure 
that  there  was  adequate  heat  exchanger  capacity  to  prevent  the  gas  from 
getting  unreasonably  hot.  The  upper  limit  design  goal  was  100°F.  In  the 
laser  system  insufficient  attention  was  paid  to  heat  exchanger  design  and 
as  a  result  the  flow  system  got  hot  to  the  touch  indicating  that  the 
heat  exchanger  was  inadequate  for  the  heat  load. 

4.5  FAN  SEAL 

The  fan  is  mounted  inside  the  flowing  gas  loop  but  to  help  keep  the 
system  clean  the  fan  motor  is  mounted  outside  the  loop.  This  necessitates 
bringing  the  driveshaft  out  of  the  vacuum  system  using  some  kind  of  vacuum 
seal.  The  fans  of  an  industrial  laser  being  developed  by  the  laser  group 
are  powered  by  a  Volvo  hydraulic  motor  with  the  shaft  sealed  by  a  lip 
seal.  The  seal  operated  well  and  was  capable  of  withstanding  vacuum. 
Industry  literature  indicates  that  a  lip  seal  should  be  capable  of  sealing 
against  a  vacuum  with  shaft  speeds  of  up  to  4000  rpm. 

Typical  values  quoted  for  rubber  are  as  follows.  F.ubber  lip  seals 
should  be  capable  of  withstanding  a  maximum  pressure  difference  of  10 
bar  and  a  peripheral  shaft  speed  of  about  5  m/sec.  Some  literature  quotes 
a  top  shaft  speed  for  neoprene  rubber  lip  seals  as  10  m/sec  with  an 
upper  temperature  limit  of  170°C.  For  a  1  inch  diameter  shaft  this  suggests 
a  top  speed  of  7600  rpm. 

The  initial  fan  estimates  indicated  that  roughly  2  HP  would  be 
needed  to  drive  the  fan  at  3500  rpm  in  air.  But  considerably  less  power 
would  be  needed  in  Helium  or  a  laser  mix  because  of  their  lower  density. 

So  the  maximum  shaft  speeds  quoted  here  mean  that  it  would  be  possible 
to  increase  the  fan  speed  using  pulleys  and  v-belts.  Table  4.2  gives 


. 


' 


81 


some  typical  shaft  speeds  for  a  1"  diameter  shaft. 


Fan  Speed 
rpm 


Shaft  vel. 
m/ sec 


9000 


12 


6000 


8 


4000 


5.3 


Table  4.2  Typical  shaft  speeds  for  1"  diameter  shaft. 

These  speeds  are  all  faster  than  are  usually  encountered  but  it  would 
appear  that  the  4000  rpm  and  6000  rpm  values  are  within  the  capabilities 
of  conventional  lip  seals.  It  was  decided  to  use  Volvo  seals  like  the 
ones  in  their  hydraulic  motor  because  they  have  somewhat  better  than 
normal  specifications.  They  are  supposed  to  be  good  for  a  top  speed  of 
10000  rpm  for  a  1"  diameter  shaft. 

Fig.  4.4  shows  the  fan  mounting  and  vacuum  seal  arrangement. 

Inside  the  system  the  fan  shaft  is  supported  by  a  self  aligning  bearing 
mounted  on  a  support.  To  minimize  problems  of  system  contamination  from 
the  bearing  it  was  packed  with  a  low  vapour  pressure  grease.  Outside  the 
system  the  shaft  is  supported  by  a  bearing  running  in  a  reservoir  of 
vacuum  pump  oil  contained  between  two  lip  seals.  Some  oil  leakage  from 
the  seal  is  inevitable  but  since  vacuum  pump  oil  is  used,  contamination 
would  not  be  much  worse  than  that  due  to  backstreaming  of  oil  from  the 
vacuum  pump . 


4.6  V-BELTS 


The  fan  is  driven  from  the  grinder  by  a  V-belt.  Although  speeds 


are  a  bit  higher  than  are  usually  encountered  by  V-belts,  Gates  has  a 
series  of  polyflex  belts  which  are  capable  of  running  this  fast.  A 
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Fig.  4.4  Fan  Assembly  and  Shaft  Seal. 
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7M825  belt  was  chosen.  Fig.  4.5  shows  the  grinder  mounting  and  pulley 
arrangement . 

4.7  VACUUM  PUMP 

Initially  the  volume  of  the  vacuum  system  was  estimated  to  be 
roughly  200  liters.  A  reasonable  pump  down  time  to  0.1  torr  would  be 
about  15  minutes  for  this  volume.  Using 


where  F  is  a  correction  factor  and  for  a  single  stage  pump  F=12,  a 
required  pumping  speed  of  S  =  160  Vmin,  or  better,  is  calculated.  A 
Sargent-Welch  model  1397  2  stage  rotary  pump  with  a  pumping  speed  of 
300  fi,/min  was  available  and  met  this  specification  so  it  was  used. 

4.8  MOTOR  SPEED  CONTROLLER 

The  motor  speed  controller  consists  of  a  simple  triac  controller 
as  diagrammed  in  Fig.  4.6.  With  this  controller  the  fan  speed  could  be 
varied  between  3600  rpm,  its  upper  limit  under  load,  and  about  1000  rpm. 
Although  it  was  possible  to  set  lower  speeds  they  were  not  practical  to 
use  as  they  tended  to  drift  too  much. 


4.9  FAN  SPEED  SENSING 

The  fan  speed  was  sensed  by  painting  a  black  strip  on  the  fan  pulley 
which  had  a  shiny  aluminum  finish  everywhere  else.  A  phototransistor 
was  then  used  to  generate  pulses  by  responding  to  the  light  reflected 
back  from  the  rotating  pulley.  The  pulses  were  counted  by  a  simple 


'  '  !  ;,rf  rJ  ,il 
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Fig.  4.5  Grinder  Mounting,  and  Pulley  Arrangement. 
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Fig.  4.6  Motor  Speed  Controller. 
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digital  rpm  counter  and  could  also  be  monitored  by  the  oscilloscope  for 
a  more  accurate  measurement. 

4.10  PULSE  DISCHARGE  CIRCUITRY 

A  pulse  discharge  circuit  was  needed  that  would  model  the  essential 
features  (and  hopefully  essential  problems)  of  the  laser  system.  Therefore 
a  spark  pin  preionized  discharge  configuration  was  chosen  (Fig.  4.7). 

Spark  pins  were  provided  both  upstream  and  downstream  from  the  discharge 
electrodes.  Since  the  energy  being  switched  in  this  system  was  much 
smaller  than  in  the  laser  system,  and  to  reduce  the  number  of  vacuum  feed¬ 
throughs,  a  thyratron  was  chosen  for  switching  rather  than  blown  spark 
gaps  in  the  gas  flow  loop.  Using  the  thyratron  also  resulted  in  more 
reliable  switching  and  much  less  electrical  noise  than  the  spark  gaps. 

The  pulse  discharge  circuit  consists  of  a  thyratron  discharge 
circuit,  a  thyratron  driver,  and  a  voltage  controlled  pulse  generator  to 
vary  the  repetition  rate.  The  circuits  are  shown  in  Figs.  4.8,  4.9  and 
4.10. 

Essentially  the  discharge  circuit  consists  of  a  storage  capacitor, 
the  discharge  electrodes,  and  a  thyratron.  The  storage  capacitor  is 
charged  up  from  a  D.C.  supply  through  the  charging  inductor.  The  thyratron 
acts  like  a  switch  and  when  it  fires  shorts  the  hot  side  of  the  capacitor 
to  ground.  Since  the  voltage  cannot  change  instantly  across  the  capacitor 
an  inverted  voltage  appears  across  the  electrodes.  A  discharge  forms 
between  the  electrodes,  completing  the  circuit,  and  the  storage  capacitor 
discharges  through  the  discharge. 

The  spark  pins  and  coupling  capacitors  to  either  side  of  the  discharge 
provide  preionization.  As  the  voltage  rises  across  the  electrodes  one 
of  the  spark  pins  will  break  down.  It's  coupling  capacitor  then  begins 


‘ 
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Fig.  A.  7  Discharge  System  Configuration. 

The  discharge  is  formed  between  two  copper  Rogowski  profiled 
electrodes  with  preionization  provided  by  spark  pins. 
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e  ischarge  pulse  energy  is  varied  by  changing  the  storage  capacitor  C  and 

s 

the  pulse  length  is  controlled  by  the  time  constant.  The  15  KSi  resistor 

was  later  replaced  by  a  1  mH  inductor  (see  Section  5.4) 
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to  charge.  When  the  voltage  becomes  high  enough  the  other  spark  pin 
will  then  break  down  and  the  second  triggering  capacitor  will  charge. 
Finally  the  voltage  rises  again  until  the  breakdown  voltage  of  the  main 
electrodes  is  reached.  The  discharge  then  forms  between  the  main 
electrodes  and  drains  the  charge  from  the  main  storage  capacitor  and  two 
coupling  capacitors  through  this  path. 

The  thyratron  is  fired  by  a  SCR  (Fig.  4.9),  which  in  turn  is 
controlled  by  the  voltage  controlled  pulse  generator  (Fig.  4.10).  This 
system  is  capable  of  switching  up  to  20  KV  pulses  and  the  design  goal  was 
to  switch  1  J  per  pulse  at  repetition  rates  of  up  to  1  KHz.  The  system 
is  capable  of  switching  in  excess  of  1  KHz  repetition  rates  and  will 
also  switch  more  than  1  J  per  pulse.  The  present  maximum  is  about  10  J 
per  pulse  at  250  Hz.  The  repetition  rate  limit  is  due  to  discharge  arcing 
rather  than  being  due  to  a  limitation  of  the  pulser. 


« 


CHAPTER  5 


FLOWING  GAS  DISCHARGE  SYSTEM  TESTING 
5.1  FLOW  SYSTEM  TESTING 

Initially  the  system  was  assembled  and  tested  without  turning  vanes 
in  the  big  elbow.  The  results  were  disappointing.  The  flow  through 
the  test  section  was  very  nonuniform,  most  of  it  being  along  the  right 
side  of  the  pipe  with  very  little  flow  in  the  center  and  the  left  side  of 
the  test  section  (Fig.  5.1).  The  maximum  velocity  along  the  wall  was 
13.4  m/sec,  as  measured  with  a  pitot  tube.  This  was  much  lower  than  it 
should  have  been.  In  addition  just  upstream  of  the  test  section,  before 
the  gas  flow  began  to  converge,  there  was  a  great  deal  of  swirl,  with  the 
gas  flow  travelling  in  a  helical  pattern  down  the  pipe.  This  was  easily 
observed  by  watching  the  action  of  a  tuft  of  string  held  in  the  flow  on 
the  end  of  a  rod. 

The  fan,  elbow,  and  heat  exchanger  combination  was  tested  to  see 
what  static  pressure  it  could  generate  with  the  elbow  outlet  blocked  off 
(Fig.  5.2).  Since  the  system  head  loss  is  high  and  the  volume  flow  low 
this  will  give  a  better  estimate  of  expected  system  performance  than  the 
extrapolation  done  previously  from  the  fan  tables.  Static  pressure  was 
measured  directly  in  inches  of  H^O  using  a  water-filled  U-tube  manometer. 

The  fan  only  generated  a  static  pressure  of  IV  H^O  which  is  considerably 
less  than  the  3h"  to  4"  H^O  that  had  been  anticipated.  In  addition  the 
gas  flow  out  of  the  elbow  had  a  lot  of  swirl  and  came  out  at  an  angle, 
as  shown  in  Fig.  5.3a. 

Two  changes  were  made  in  a  successful  attempt  to  remedy  these  problems. 
They  were  the  addition  of  two  turning  vanes  in  the  large  elbow  and  the 
reduction  of  clearance  between  the  fan  blades  and  the  elbow  wall.  The 
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Fig.  5.1  Initial  Gas  Flow  Velocity  Distribution  in  the 


Test  Section. 


. 


Closed 

Chamber 
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Fig.  5.2  Apparatus  for  Testing  Fan  Static  Pressure. 
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addition  of  the  turning  vanes  was  probably  the  more  significant  change 
of  the  two.  Remeasuring  the  static  pressure  generated  by  the  fan  showed 
that  it  had  risen  to  1  15/16"  H^O,  an  improvement  of  almost  3/4"  H^O. 

In  both  tests  the  fan  speed  was  2940  rpm  and  the  motor  drew  15. 8A,  which 
was  full  power.  In  addition  the  flow  out  of  the  elbow  was  now  uniform 
and  swirl  free  as  shown  in  Fig.  5.3b. 

This  eliminated  swirl  problems  downstream  and  the  flow  through  the 
test  section  became  uniform.  The  flow  velocity  through  the  test  section 
rose  to  22  m/sec  which  is  a  good  improvement  although  it  is  still  only 
half  of  what  was  estimated  in  the  design  stage. 

After  the  system  was  completed  and  tested  the  losses  were  recalculated 
using  the  final  configuration  (Fig.  4.1)  instead  of  the  crude  model 
(Fig.  4.2).  For  a  flow  rate  of  225  CFM,  which  correspondence  to  22  m/sec 
in  the  test  section,  the  calculated  system  loss  came  to  1.18"  H^O  which 
is  a  negligible  change  from  the  estimate  of  1.23"  H^O  calculated  for  the 
crude  system  model. 

These  results  show  that  there  are  two  reasons  for  the  reduction  in 
flow  velocity  from  the  initial  estimate.  The  main  reason  is  that  the 
extrapolation  of  static  pressure  from  the  fan  rating  tables  resulted  in 
an  estimate  that  was  much  too  high,  3V  H^O  instead  of  the  actual  value 
of  just  under  2"  ^0.  The  other  reason  is  that  the  actual  system  head 
losses  were  somewhat  larger  than  estimated.  This  is  not  surprising 

since  the  tables  used  to  obtain  the  friction  coefficients  were  for 
galvanized  air  ducting  and  this  system,  which  is  constructed  of  welded 

aluminum  pipe,  has  a  rougher  finish  which  would  mean  higher  friction  co¬ 
efficients.  The  test  section  also  would  have  a  higher  head  loss  than 
estimated  since  the  estimate  assumed  that  it  was  a  smooth  duct  when  in 
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fact  it  contained  recessed  wells  with  electrodes  sitting  in  them. 


5.2  GAS  FLOW  CALIBRATION 

In  designing  flowing  gas  systems  the  fan  laws  are  used  to  estimate  the 
fan  capability  required  to  achieve  a  given  gas  flow.  The  laws  are  as 
follows : 


D  RPM 

CFMa  =  CFMb  x  fcS-)  x  (^) 


D  -  RPM  „  6 

P  =  Pu  X  tV  X  (— rS-)2  X  (“T“) 


'D, 


'RPM^ 


D  RPM  6 

HP  =  HP  x  (-^T  x  (==-£)  J  x  (t^) 


■RPMb' 


where  CFM  is  gas  flow 

D  is  the  fan  diameter 
RPM  is  fan  speed 

P  is  pressure  head  generated  by  the  fan 
6  is  the  gas  density 

This  says  that  the  gas  flow  is  linearly  proportional  to  fan  speed 
and  is  independent  of  the  gas  density.  The  power  needed  to  move  the  gas 
and  the  pressure  generated  do  however  vary  with  density. 

Although  this  is  adequate  for  system  design  it  is  true  only  as  a 
first  approximation.  In  this  system  the  flow  velocity  for  a  given  fan 
speed  does  indeed  depend  on  the  gas  density  and  it  is  not  linear  with  fan 
speed.  Measurements  indicated  that  although  the  relationship  was  linear 
for  fan  speeds  above  1500  rpm,  it  was  not  for  lower  speeds.  Since  gas 
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flow  velocity  is  monitored  by  measuring  the  fan  speed  it  is  therefore 
necessary  to  get  calibration  curves  for  the  gas  flow  velocity  as  a 
function  of  fan  speed. 

The  gas  flow  velocity  between  the  electrodes  was  measured  by  using 
a  pitot-static  tube  and  the  Baratron  capacitance  manometer.  The  fan 
speed  is  observed  by  a  phototransistor  which  detects  a  timing  mark  on  the 
fan  pulley.  The  pulses  output  by  the  phototransistor  are  monitored  on 
the  scope  and  the  time  between  two  pulses  gives  the  period  of  rotation. 

The  fan  speed  in  rpm  can  then  be  determined  from  the  period,  and  calibration 
curves  plotted  which  give  gas  flow  velocity  as  a  function  of  fan  speed. 

Rather  than  expressing  the  calibration  curves  in  the  form  of  a  table 
it  is  more  convenient  to  represent  them  as  equations  giving  flow  velocity 
as  a  function  of  fan  speed.  Then  the  calibration  information  can  be 
easily  used  by  a  computer  program  for  further  data  analysis.  In  order  to 
generate  a  calibration  curve  it  was  necessary  to  try  and  fit  a  smooth 
curve  to  the  data  points  but  the  question  was  which  is  the  best  curve 
to  use.  The  IMSLIB  library  on  the  computer  has  a  routine  which  can  fit 
any  degree  of  polynomial,  up  to  degree  10,  to  a  given  set  of  data  points. 
Using  this  subroutine  polynomials  of  up  to  degree  4  were  fitted  to  the 
calibration  data.  The  program  works  by  doing  a  least  squares  fit  of  a 
series  of  orthogonal  polynomials  to  the  data  to  generate  an  equation 
fitting  the  data  points.  To  determine  how  good  the  fit  is,  the  variance 
must  be  calculated  for  the  data  points  and  the  predicted  equation.  The 
error  sum  of  square  is 

m  ~ 

=  l  (yi  -  yi  ) 

j=l  J  Jpredicted 


. 

f 

I 

. 
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and  the  variance  is 


a . 
1 


m-i-1 


where  i  is  the  degree  of  polynomial  and  m  is  the  number  of  data  points. 
Notice  that  you  are  restricted  to  i<m-2.  When  an  increase  in  the  degree 
of  the  polynomial  produces  no  significant  change  in  the  variance  then 
the  best  fit  has  been  achieved.  Using  this  technique  it  was  found  that 
the  calibration  results  were  best  described  by  a  cubic  equation.  The 
coefficients  for  the  equation  are  given  in  Table  5.1,  and  the  calibration 
curves  are  plotted  in  Figures  5.4  to  5.7.  The  curve  fitting  program  is 
listed  in  the  Appendix. 


2  3 

y  =  A  +  Bx  +  Cx  +  Dx 


gas 


B 


D 


60  torr  N, 


0.446222x10 


80  torr  CO^  -0.106198 
700  torr  He  -0.174847 


_1  0.255427xl0~2  0 . 170707xl0_5 

0.422460xl0”2  0 . 176471xl0~5 

0.356737xl0"2  0 . 168749xl0_5 


-0.192712x10 


-0.257645x10 


-9 


-0.206890x10 


-9 


400- 


700  torr  Mix  -0.649503x10 


0.456847x10  2  0.158747x10  5 


-0.231660x10 


-9 


Table  5.1  Calibration  Curve  Coefficients 

5.3  SYSTEM  VACUUM  INTEGRITY 

A  measurements  of  system  vacuum  integrity  was  important  because 
knowledge  of  the  leak  rate  gives  a  worst  case  estimate  of  impurity  levels 
that  can  be  expected  in  the  system  during  an  experiment  due  to  leaks. 


. 
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Fig.  5.4  Gas  Flow  Velocity  Calibration  Curve  for  60  torr  N2 . 


VELOCITY  M/SEC 
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Fig.  5.5  Gas  Flow  Velocity  Calibration  Curve  for  80  Torr  CO^ . 


VELOCITY  M/SEC 
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Fig.  5.6  Gas  Flow  Velocity  Calibration  Curve  for  700  Torr  He. 
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Fig.  5.7  Gas  Flow  Velocity  Calibration  Curve  for  a  1:1:8  002^2^6 
Laser  Mixture  in  the  Range  from  400  to  700  Torr. 
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First  the  system  volume  was  found.  A  one  liter  flask  of  air 
at  700  torr  pressure  was  attached  to  the  system  with  a  connecting  valve. 
Then  the  system  was  evacuated  and  the  valve  opened  to  allow  air  from  the 
flask  to  expand  into  the  system.  The  resulting  pressure  was  then  noted 
and  from  this  information  the  system  volume  could  be  calculated  using 
the  gas  law: 


P  V  = 
1  1 


P  V 
2  2 


The  system  volume  was  found  to  be  180  liters. 

The  leak  rate  is  determined  by  pumping  out  the  system  and  sealing 
it  off,  then  observing  the  pressure  rise  over  a  period  of  time.  It  is 
calculated  using  the  formula: 

i  T7  AP 

leak  rate  =  V  — — 


where  V  is  system  volume,  AP  is  rise  in  pressure,  and  AT  is  time.  The 

,  ,  _  .  .  ^  ,  ,  £  micron  _  ,  i  torr  . 

leak  rate  for  this  system  was  typically  1  -  or  3.6  — - -  •  This 

J  J  r  j  sec  hr 

was  a  rise  of  0.2  torr  in  10.2  hours.  So  the  system  has  very  good  vacuum 
integrity  and  contamination  due  to  leaks  is  not  a  problem. 


5.4  DISCHARGE  CIRCUIT 

The  thyratron  pulse  circuit  operated  fairly  reliably.  It  is  capable 
of  switching  120  nf  at  10  Kv  charging  voltage  and  repetition  rates  of 
300  pps.  The  limit  here  is  due  to  discharge  arcing  rather  than  a 
limitation  of  circuit  capability.  Smaller  capacitors  can  be  switched  at 


. 
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higher  repetition  rates  and  20  nf  at  10  Kv  charging  was  switched  at 
repetition  rates  of  up  to  1  KHz  without  problems. 

Several  modifications  were  made  to  the  circuit  during  the  course  of 
the  experiments.  Originally  there  were  two  preionization  spark  pins. 

The  experiments  demonstrated  that  discharge  products  from  the  upstream 
spark  pin  adversely  affected  the  discharge  uniformity  so  the  upstream 
spark  pin  was  removed  from  the  system. 

The  other  change  involved  replacing  what  was  originally  a  discharge 
shunt  resistor  with  a  shunt  inductor  instead.  The  shunt  element  is  needed 
in  order  for  the  main  storage  capacitor  to  charge  up.  Although  the 
shunt  resistor  worked  fine  with  He  it  was  responsible  for  limiting  capacitor 
values  to  40  nf  or  less  for  and  CO^ .  The  problem  was  that  the 
thyratron  went  into  continuous  conduction  and  was  traced  to  the  fact  that 
the  discharge  in  these  gases  would  extinguish  before  the  capacitor  had 
discharged  completely.  When  the  discharge  extinguished  the  residual 
voltage,  typically  200  V  for  CO^ ,  would  be  sufficient  to  force  the 
thyratron  into  continuous  conduction.  In  order  to  recover  the  thyratron 
needs  zero  or  reverse  voltage  across  it  for  a  minimum  of  about  50  usee. 

This  is  to  allow  the  charged  gas  in  the  tube  to  recombine  to  the  point 
of  being  non-conducting.  Replacing  the  shunt  resistor  by  the  shunt 
inductor  allowed  the  residual  charge  to  drain  off  the  capacitor  and  caused 
the  thyratron  anode  to  ring  slightly  negative  which  helped  aid  its  recovery. 

A  similar  problem  was  encountered  with  the  SCR  in  the  thyratron 
drive  circuit  locking  on  occasionally,  usually  as  a  result  of  noise  from 
high  currents  in  the  discharge  circuit  during  arcs.  This  problem  was 
traced  to  the  voltage  being  reapplied  to  the  SCR  too  quickly.  The  SCR 
remained  reverse  biased  for  only  about  40  usee  but  it  requires  typically 
75  ysec  to  recover.  Boosting  the  charging  inductor  from  7.5  mH  to  2.6  H 


' 
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increased  the  reverse  bias  period  to  over  100  ysec  which  was  sufficient 
for  the  SCR  to  recover  reliably. 

Typical  voltage  and  current  waveforms  for  the  discharge  are  shown 
in  Fig.  5.8.  These  waveforms  are  for  a  discharge  in  60  torr  but  the 
waveforms  in  other  gases  were  similar.  The  initial  current  spike  is 
due  to  the  spark  pin  coupling  capacitor  discharging  while  the  main  hump 
is  due  to  the  main  storage  capacitor  discharging. 
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5.8  Typical  Voltage  and  Current  Waveforms  for  a  60  Torr  N 
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Discharge . 


The  discharge  capacitor  was  40  nf  and  the  charging  voltage 


was  10KV. 


CHAPTER  6 


RESULTS  AND  CONCLUSIONS 

6.1  DEFINITIONS  AND  THEORY 

In  a  high  repetition  rate  pulsed  glow  discharge,  discharge  products 

are  removed  between  pulses  by  convection.  If  they  are  not  removed  arcing 

results.  The  gas  convection  time,  t  , 

conv 


Figure  6.1  Schematic  of  discharge  region  showing  terms  used  to  define  t 


conv 


is  the  time  taken  for  the  gas  to  traverse  the  discharge  region 


=  W 

conv  v 


and  is  the  inverse  of  the  ideal  repetition  rate,  f.,  ,  The  minimum 

ideal. 

time  between  arc  free  pulses  is  the  clearing  time,  x  ,  ,  and  is  the 

inverse  of  the  maximum  repetition  rate.  A  clearing  ratio  is  defined  to 
compare  the  observed  repetition  rate  with  the  gas  convection  time. 


CR  = 


minimum  time  between  arc  free  pulses 

time  for  gas  to  traverse  gap 


1 


PRF 


v 


W/v  (PRF) W 


6.1.1 
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where  PRF  is  the  maximum  observed  pulse  frequency,  or  -  .  The 

clear 

inverse  of  the  clearing  ratio  is  the  fractional  frequency 


■  i_  =  W(PRF) 
fractional  CR  v 


Usually  the  fractional  frequency  is  less  than  one,  that  is,  the 
repetition  rate  achieved  is  less  than  that  predicted  by  taking  the  time 
for  gas  to  traverse  the  electrodes  as  the  interpulse  period.  Dzakowic  and 
Wutzke  (38)  attribute  the  lower  repetition  rate  to  gas  expansion  due  to 
heating  from  the  energy  deposited  in  it  during  the  discharge.  They  predict 
that  the  resulting  clearing  ratio  is  given  by: 


CR  =  1  +  j  3[(1  +  e)1/Y  -  1] 


6.1.3 


(l-n)E. 

where:  e  =  — = - —  is  the  normalized  energy  density, 

p  C  T  vol 
v 


6  is  the  ratio  of  upstream  to  interelectrode  flow  areas, 


p 

y  =  — —  is  the  ratio  of  specific  heat  capacities, 
V 


p  is  the  density  of  the  unheated  gas, 


T  is  the  temperature  of  the  unheated  gas, 


n  is  the  electrical  to  optical  energy  conversion  efficiency 

so  (1-n)  is  the  fraction  of  input  electrical  energy  density, 

E. 

l 


vol 


,  that  appears  as  heat. 
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Figure  6.2  Geometry  for  derivation  of  gas  expansion  clearing  ratio 
(equation  6.1.3) 

A  simplified  derivation  will  show  how  this  expression  was  obtained. 
Consider  a  channel  with  cross  sectional  area,  A,  and  a  volume  of  gas  WA 
in  the  channel.  Suppose  energy  is  suddenly  deposited  by  a  discharge  in 
the  gas  along  the  width  W.  In  that  case  the  gas  will  heat  up  and  the 
pressure  in  region  WA  will  increase  but  the  density  will  stay  constant 
because  the  gas  molecules  have  not  yet  had  time  to  move. 


P  =  p  RT 
so 


The  energy  deposited  by  the  discharge  is  —  kAT  per  molecule  where  AT  is 
the  change  in  temperature.  So  the  total  energy  deposited  for  N  molecules 
is : 


AE  =  |  kNAT 


Ill 


The  new  pressure  then,  is 


2  T  +  AT  AT 

P  =  —  P  =  - - —  P  =  (1  +  — )  P 

2  T  1  T  1  ^  T  '  1 


-  &  +  r— >  pi 
2  kTN 


The  gas  will  then  expand  because  the  pressure  in  the  heated  gas  is  higher 
than  the  surrounding  pressure.  Assuming  an  adiabatic  expansion: 


P  VY  =  P  VY 

2  2  final  final 


and  let 


AE 


|  kTN 


then 


V 


Y 


P  VY 
2  2 


final  P 


final 


Y 


(1  +  OP^ 

^f inal 


Assuming  an  essentially  infinite  duct,  expansion  will  continue  until 

P  =  P  so: 

final  1 


vT. 


-  (1  +  c)  Vj  -  (1  +  e)  Vi 


Y 


f  inal 


ial 
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but 


V  =  WA 

initial 
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and 


V 


final 


(W  +  2x) A 


so 


W  +  2x  =  (1  +  e)1/Y  W 


or 


x  - 


f  [a  + 


The  clearing  time  will  be  increased  from  x  by  the  amount  of  time 

conv 


needed  to  clear  the  gas  the  additional  distance  x  past  the  electrodes 


=  x  +  W  =  x  +  W 
clear  v  v  v 


CR 


Tclear 

T 

conv 


W  +tt[(l  +  £)1/Y  -  1] 

v  _  _ 2 _ 

W/v 


-  1  +  i[(l  +  e)1/Y  -l] 


6.1.5 


The  expression 


AE _ 

|  kTN 


can  be  put  in  a  form  corresponding  to  6.1.3. 
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e 


AE _ 

fkTN 


AE 

vol 


vol 
I  kTN 


but 


,  _  R_  _  R _  _  Rn  _  Rm 

Na  N/n  N  ~  M 

Where : 

n  is  number  of  gram  moles 
N  is  number  of  molecules 
m  is  mass  of  a  molecule 
M  is  molecular  weight 
so 


but 

so 


AE  vol 


wl  1  d  !  TN 

2  R  M  TN 


vol 

mN 


c 

V 


3 

2 


R  for  a  monatomic 


gas 


e 


AE 

vol 


P 


1 


c 


so  6.1.5  is  essentially  the  same  as  Dzakowic  and  Wutzke's  expression, 
6.1.3,  except  that  the  (1  -  ri)  term  has  been  ignored  and  3  was  set  equal 
to  1. 

A  slightly  different  expression  is  applicable  for  the  gas  expansion 
in  the  flowing  gas  discharge  system  since  the  discharge  has  a  cylindrical 
geometry  (Fig.  6.3). 


■ 
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Figure  6.3  Geometry  of  cylindrical  discharge 


Going  back  to  6.1.4 


V  =  (1  +  e)^Y  V 

final  V  ;  initial 


V.  =  TT  r  l 

initial  1 


Vf inal  =  "  r2£ 


SO 


2  1/y  2 

it  =  (1  +  e)  it  r^£ 


2  ( i  x  ^/Y  2 

r2  =  (1  +  £)  ri 


r„  = 


(1  +  e)1/y  r 


1 
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r 


2 
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+  x  =  \/(l  +  r. 


x  =  ^[^/(l  +  £)1//y  -l] 


so  the  clearing  ratio  is 


CR= 


x  +  2r. 


v 

2r. 


x  +  2r. 
2r 


v 


rx[  VU  +  e)1/Y  -  l]  +  2rx 


1+|  [-/(l  +  e)  1/y  -  l] 


6.1.7 


If  the  interelectrode  and  upstream  cross-sectional  areas  are 
different  the  expression  for  clearing  time  is  modified. 


x  W 


clear  v/3  v 


where  3  is  the  ratio  of  upstream  to  interelectrode  areas.  Then  the 


clearing  ratio  becomes 

W  rC'v/q  +  e)1/Y  -l] 
v/3  v  v/3 


+ 


CR  = 


2r 


v 


W 


v 


2r 

v 


=  1  +  j  BtVd  +  e)1/Y  -  l] 


6.1.8 
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In  the  case  of  Dzakowic  and  Wutzke  equation  6.1.3  is  obtained  instead 
of  equation  6.1.5. 

For  small  energy  densities  the  clearing  ratio  should  tend  to  1  so 

the  repetition  rate  should  tend  to  f . ,  ..  according  to  the  above  formulas. 

ideal 

However  Dzakowic  and  Wutzke  argue  that  the  upper  limit  on  repetition  rate 
is 


f 


ideal 

JT 


,  or  CR  =  /J  . 


This  is  because  in  the  boundary  layer  convection  effects  are  small  and 
the  discharge  products  are  removed  by  diffusion  to  the  electrode  surfaces 
where  recombination  then  occurs.  The  diffusion  length,  6  ,  is  defined  as 
the  distance  into  the  boundary  layer  from  the  electrode  surface  at  which 
both  diffusion  time  and  convection  time  are  equal.  The  diffusion  time 
is  given  by: 


Tdiff 


6 


2 

D 


Ac 


where  X  is  the  mean  free  path  and  c  the  mean  speed  of  the  particles.  Then 
assuming  a  laminar  boundary  layer  they  calculate  a  diffusion  time  of: 


=  /311  or  CR  =  /3  . 

dif  f  v 

Notice  that  the  low  energy  limit  is  independent  of  gas  type,  gas  density, 
and  energy  density. 

Although  their  data  supported  this  low  energy  density  repetition 
rate  limit  for  energy  densities  less  than  about  30  J/i  my  own  experiments 


. 


' 

did  not.  For  the  flowing  gas  discharge  system,  if  the  energy  density 
was  low  enough,  several  pulses  per  gas  change  could  be  achieved  before 
an  instability  occurred.  As  a  result  CR<1  and  hence  repetition  rates 
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>f.,  were  observed, 
ideal 

ratios  of  <1. 


Olson  et  al  (21)  also  report  achieving  clearing 


6.2  GAS  TYPES 

The  discharge  characteristics  of  the  three  CO^  laser  gases  CO^,  , 
and  He  were  investigated  separately  and  in  combination.  This  is  because 
theoretical  calculations  are  easier  for  the  individual  gases  than  for  a 
mixture  and  also  to  see  if  the  repetition  rate  effects  could  be  attributed 
to  one  particular  gas  of  the  three.  The  repetition  rate  effects  are 
more  likely  to  be  common  to  all  three  gases  or  to  be  due  to  one  gas  than 
they  are  to  be  due  to  some  special  combination  of  the  three. 

For  the  individual  gases  60  torr  N  ,  80  torr  CO^  and  700  torr  He 
were  chosen.  The  different  pressures  were  chosen  because  they  all  had 
a  10  KV  breakdown  voltage.  The  three  gases  had  very  different  discharge 

characteristics.  Helium  was  the  most  stable  and  discharges  could  be  easily 
formed.  It  formed  the  largest  diameter  discharge  and  at  high  repetition 
rates  tended  to  form  a  nonuniform  glow  rather  than  a  bright  arc.  CC^ 
formed  the  smallest  diameter  discharge  and,  for  low  energy  densities, 
the  repetition  rate  induced  arc  tended  to  be  between  the  electrodes  rather 
than  from  the  cathode  to  the  spark  pin.  formed  a  discharge  between 

CO^  and  He  in  diameter.  A  unique  feature  of  was  the  production  of 
long  lived  metastables,  characterized  by  a  long-lived  orange  glow,  which 
were  accompanied  by  high  electron  densities  (38).  This  tended  to  promote 
tracking  problems  in  N^,  and  usually  arced  from  the  cathode  to  the 


. 
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downstream  spark  pin  rather  than  between  the  electrodes. 

6.3  EFFECT  OF  SPARK  PINS  ON  THE  SYSTEM 

In  the  laser  system  it  was  observed  that  mounting  the  preionizing 
spark  pins  upstream  from  the  discharge  appeared  to  cause  discharge 
arcing,  presumably  because  of  discharge  products  from  the  spark  pins 
blowing  downstream  into  the  discharge.  An  experiment  was  done  in  the 
flowing  gas  discharge  system  to  verify  the  effect. 

Two  spark, pins  were  mounted  in  the  system,  one  upstream  and  one 
downstream  from  the  discharge.  Discharges  were  then  observed  in  700  torr 
He,  60  torr  N^ ,  80  torr  CO^,  a  He-N^  mixture,  and  a  He-CO^  mixture. 
Discharge  uniformity  and  maximum  repetition  rate  attainable  before  bright 
arcs  occurred  were  compared  for  the  system  with  both  spark  pins  and 
with  just  the  downstream  spark  pin  operating.  The  maximum  repetition 
rates  both  with  and  without  the  upstream  spark  pin  are  shown  in  Figures 

6.4  to  6.8.  As  can  be  seen,  the  upstream  spark  pin  had  no  detrimental 
effect  on  the  He  discharge  but  in  N^ ,  CO^,  and  the  N^-He  and  CO^-He  gas 
mixtures  the  maximum  repetition  rate  was  noticeably  lower. 

For  all  the  runs  the  storage  capacitor  was  1960  pf  and  the  charging 
voltage  10  KV.  This  is  an  energy  density  of  about  0.1  J  or  an  energy 
density  of  about  1.2  J/£  for  each  pulse.  Because  the  energy  density  was 
so  low  it  was  possible  to  achieve  very  high  repetition  rates.  For 
instance,  in  700  torr  He  for  a  fan  speed  of  3500  rpm  the  gas  flow  velocity 
was  about  26  m/sec.  For  a  clearing  ratio  of  1  the  maximum  repetition 
rate  that  could  be  achieved  would  be  300  to  400  Hz  instead  of  close  to 
2000  Hz  as  was  achieved.  Since  the  energy  density  was  so  low  several 
pulses  could  be  produced  per  gas  change  before  arcing  occurred. 
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Fig.  6.4  Effect  of  Upstream  and  Downstream  Spark  Pins  on  Arc 
Free  Pvepetition  Rate  for  700  Torr  He. 
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Fig.  6.5  Effect  of  Upstream  and  Downstream  Spark  Pins  on  Arc  Free 
Repetition  Rate  for  60  Torr  1SL . 


REP  RATE  HZ 


121 


2000  -i 


N2  HE  MIX 

©  1  SPARK  PIN 
a  2  SPARK  PINS 


1600  - 


1200  - 


800  - 


o 

0 


400  - 


0 

0 


— r 


o 

o 


T 


1000 

FAN 


2000 

SPEED  RPM 


o 

© 


a 

A 


- 1 - 

3000 


A 

A 


CR=1 


4000 


Fig.  6.6  Effect  of  Upstream  and  Downstream  Spark  Pins  on  Arc  Free 

Repetition  Rate  for  a  Mixture  of  40  Torr  and  120  Torr  He. 
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Fig.  6.7  Effect  of  Upstream  and  Downstream  Spark  Pins  on  Arc  Free 
Repetition  Rate  for  80  Torr  CO^- 
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Fig.  6.8  Effect  of  Upstream  and  Downstream  Spark  Pins  on  Arc  Free 

Repetition  Rate  for  a  Mixture  of  40  Torr  C02  and  200  Torr  He. 
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Besides  affecting  the  maximum  repetition  rate  in  CO^  and  ,  the 
discharge  uniformity  was  noticeably  affected  whereas  the  uniformity  of 
the  He  discharge  was  not.  In  the  discharge  contracted  into  a  narrow 
band  downstream  from  the  spark  pin  rather  than  filling  the  electrode 
(Fig.  6.9).  In  CO^  the  effect  was  even  more  interesting.  The  discharge 
formed  a  constricted  glow  column  approximately  1  cm  in  diameter  and 
located  near  the  upstream  edge  of  the  electrodes.  For  a  fixed  gas  flow 
increasing  the  repetition  rate  would  lead  to  two,  then  three,  four  and 
five  columns  spaced  one  behind  the  other  (Fig.  6.10). 

By  observing  the  number  and  spacing  of  the  CO^  discharge  columns  for 
two  fan  speeds  the  gas  flow  velocity  could  be  estimated  and  compared  with 
the  calibration  values.  The  flow  velocity  estimates  are  shown  in  Table  6.1 
and  are  based  on  the  results  displayed  in  Table  6.2. 


fan  speed 
rpm 

calibration  flow 
mm/msec 

observed  flow 
mm/msec 

3529 

i 

25±4 

2000 

! 

I  13  ! 

12±2 

] 

i 

!  ! 

!  \ 

f  i 

L 

Table  6.1  Estimates  of  flow  velocity  based  on  the 
observation  of  column  discharges  in 


The  discharge  spacing  was  measured  with  a  ruler  placed  against  the 
discharge  chamber  window,  and  the  results  are  subject  to  a  parallax  error 
of  about  ±  1/16".  In  addition,  the  spacing  of  the  columns  was  not 
uniform  but  rather  tended  to  decrease  toward  the  downstream  edge  due  to 
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Fig. 


6.9 


Cross-sectional  View  of  Constricted  Discharge 
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Fig.  6.10  Column  Discharges  in  CO^ . 


(a)  Shows  a  single  column  and 


(b)  shows  multiple  columns. 
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Fan  speed  3529  rpm 


#  of 

approx. 

discharge 

A 

estimated 

glow  columns 

spacing  inches 

period  msec 

flow  velocity 
mm/msec 

1 

2.75 

2 

2 

2.25 

23 

3 

1  1% 

1.55 

25 

4 

h  i  ik 

1.10 

29 

5 

!  3  3  3 
^444 

0.80 

24 

Fan  speed  2000 

rpm 

#  of 

approx. 

discharge 

estimated 

glow  columns 

spacing  inches 

period  msec 

flow  velocity 
mm/msec 

1 

5.3 

] 

2 

2 

4.5 

11 

3 

2.9 

11 

4 

XI* 

2.0 

13  ! 

j 

5 

3  i  7  5 

8^88 

1.6 

14  | 

Table  6.2  Discharge  column  spacing 

Since  the  column  spacing  decreases  downstream  the  flow 
velocity  estimate  uses  the  largest  column  spacing.  It  was 
obtained  by  taking  the  largest  spacing  and  dividing  it  by 
the  discharge  period. 
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electrode  curvature.  But  the  data  is  adequate  to  show  that  the  multiple 
discharge  spacing  corresponds  to  an  ionized  plug  of  gas  being  blown 
downstream.  The  initial  plug  comes  from  the  upstream  spark  pin  and 
results  in  a  highly  localized  discharge.  Subsequent  discharges  then  form 
in  this  plug  as  it  blows  downstream  because  of  its  localized  high  pre¬ 
ionization  density,  resulting  in  the  multiple  glows.  Since  CO^  forms 
multiple  glow  columns  in  the  presence  of  an  upstream  spark  pin  instead  of 
bright  arcs  the  maximum  repetition  rate  for  two  spark  pins  shown  in  Fig.  6.7 
is  the  maximum  repetition  rate  for  a  single  glow  column  rather  than  for 
a  bright  arc. 

This  experiment  shows  that  products  from  an  upstream  spark  pin  can 
be  detrimental  to  the  operation  of  a  glow  discharge.  In  particular  CO^j 
N^,  and  mixtures  containing  them  are  affected  while  He  is  not.  This 
could  be  because  the  detrimental  products  are  not  generated  in  pure  He 
or  because  He  is  immune  to  their  effects. 

6.4  REPETITION  RATE  VERSUS  GAS  FLOW  VELOCITY 

Maximum  repetition  rate,  for  a  uniform  discharge  with  no  arcing, 
increases  linearly  with  gas  flow  velocity.  This  was  verified  for  each 
of  the  three  laser  gases  and  for  discharge  energy  densities  ranging 
from  10  to  60  J/&.  Typical  results  are  shown  in  Figures  6.11  to  6.19. 

The  dashed  lines  represent  the  ideal  maximum  repetition  rate,  while  the 
circles  are  observed  data  points.  The  solid  line  is  a  least  squares  fit 
to  the  data  points  constrained  to  pass  through  the  origin. 

Although  repetition  rate  increases  linearly  with  gas  flow,  it  is 
typically  lower  than  the  ideal  repetition  rate.  This  is  in  agreement 
with  observations  by  other  researchers  (IB,  19).  Also,  as  was  pointed 
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Fig.  6.11  Repetition  Rate  Versus  Gas  Flow  Velocity. 

Data  for  60  Torr  N2  and  8  3/1  Energy  Density. 
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Fig.  6.12  Repetition  Rate  Versus  Gas  Flow  Velocity. 

Data  for  60  torr  N2  and  29  J /l  energy  density. 


■ 


REP  RATE  HZ  *10 1 

10.00  20.00  30.00  40.00  50.00 


130 


Fig.  6.13  Repetition  Rate  Versus  Gas  Flow  Velocity. 

Data  for  60  torr  N2  and  58  J/£  energy  density. 
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Fig.  6.14  Repetition  Rate  Versus  Gas  Flow  Velocity. 

Data  for  80  torr  C0?  and  15  J/£  energy  density. 
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Fig.  6.15  Repetition  rate  versus  gas  flow  velocity. 

Data  for  80  torr  and  23  J/£  energy  density. 
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Fig.  6.16  Repetition  Rate  Versus  Gas  Flow  Velocity. 

Data  for  80  torr  C02  and  35  J/£  energy  density. 
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Fig.  6.17  Repetition  Rate  Versus  Gas  Flow  Velocity. 

Data  for  700  torr  He  and  11  J/i  energy  density. 
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Fig.  6.18  Repetition  Rate  Versus  Gas  Flow  Velocity. 


00 


Data  for  700  torr  He  and  37  J / i  energy  density. 
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Fig.  6.19  Repetition  Rate  Versus  Gas  Flow  Velocity. 

Data  for  700  torr  He  and  64  J f%  energy  density. 


137 


out  in  section  6.3,  if  the  energy  density  is  low  enough  the  ideal 
repetition  rate  can  be  exceeded. 

6.5  REPETITION  RATE  VERSUS  ENERGY 

In  this  section  the  variation  of  maximum  obtainable  repetition 
rate  for  a  uniform  glow  is  studied  as  a  function  of  energy.  The 
repetition  rates  are  expressed  as  fractional  frequencies  for  ease  of 
comparison  with  theory. 

The  simplest  estimate  of  maximum  uniform  glow  repetition  rate  is 

given  by  f which  can  be  calculated  from  the  discharge  diameter  and 

the  flow  velocity.  As  was  shown  in  section  6.4  the  observed  repetition 

rates  are  usually  lower  than  f.n  . .  Since  f . ,  .  is  a  linear  function 

ideal  ideal 

of  flow  velocity  it  can  be  expressed  as: 


f 


ideal 


av 


Similarly,  if  the  observed  repetition  rate  is  also  a  linear  function  of 
flow  velocity,  as  appears  to  be  the  case  from  the  data  in  section  6.4, 
it  can  also  be  expressed  as: 


f  ,  =  bv 

obs 


Then  the  fractional  frequency  (equation  6.1.2)  would  be: 


bv  b 
fract  av 


a 


6.5.1 


■ 


. 
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Figure  6.20  Calculating  an  average  fractional  frequency  from 
observed  data. 

for  any  value  of  v.  That  is,  the  fractional  frequency  is  equal  to  the 
ratio  of  the  slope  of  the  line  for  the  observed  repetition  rate  and  the 
slope  of  the  line  for  the  ideal  repetition  rate. 

Individual  observations,  however,  are  subject  to  statistical  vari¬ 
ation  which  introduces  a  random  error  in  the  value  of  f_  calculated. 

f  ract 

But  an  average  value  of  f,.  can  be  obtained.  Given  an  observation 

°  fract 

at  a  particular  velocity  v^  the  ideal  repetition  rate  would  be: 

f  =  av 

ideal  1 


. 
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and  the  observed  repetition  rate  would  be: 


f  =  bv.  +  e. 
obs  1  1 


where  e_^  is  the  error  term  representing  the  departure  of  the  observed 
value  from  the  linear  equation  due  to  random  errors.  The  fractional 
frequency  is: 


f 


f ract . 
i 


f 


F 


obs . 

l 


ideal . 

l 


bv .  +  e . 
i  l 

av . 

l 


and,  for  n  observations,  the  average  value  of  fractional  frequency  is: 


_  n 

f  -  =  ~  l  fr 

fract  n  .  ,  fract. 

i=l  l 


..  n  b  ,+e. 
=  -  S  v1  1 

n  .  1  av . 
i=l  i 


6.5.2 


a  n 


n 


I 

i=i 


e . 

l 


av. 

l 


Equation  6.5.2  describes  the  average  fractional  frequency  calculated 

from  the  observed  data.  Although  it  differs  from  the  fractional  frequency 

1  n  e . 

-L  v*  i 

for  the  linear  model  (equation  6.5.1)  by  the  error  term  —  i  -  they 

n  ,  ..  av . 
i=l  i 

are  approximately  equal  since  the  random  errors  should  tend  to  cancel  and 
the  error  term  become  small.  So  the  average  fractional  frequency  can 
still  be  thought  of  as  approximately  equivalent  to  the  slope  of  a  line 
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fitted  through  the  data  points  and  constrained  to  pass  through  the  origin 
divided  by  the  slope  of  the  line  representing  the  ideal  repetition  rate. 

For  a  given  energy  input  density  this  gives  an  average  fractional 
frequency,  averaged  over  a  range  of  gas  flow  velocities.  The  ideal 
repetition  rate  has  a  fractional  frequency  of  1  and  the  average  fractional 
frequency  shows  how  much  the  repetition  rate  has  dropped  below  the  ideal 
case  for  a  given  energy  density.  This  provides  a  convenient  means  for 
comparing  the  theory  of  Dzakowic  and  Wutzke  with  the  observed  repetition 
rates . 

Figures  6.21  to  6.23  show  the  observed  fractional  frequencies  for 
a  uniform  glow  with  no  arcing  and  the  predicted  fractional  frequencies 
using  Dzakowic  and  Wutzke' s  theory  (equation  6.1.8).  The  observed 
fractional  frequency  falls  off  faster  than  predicted  as  energy  density 
increases . 

6.6  REPETITION  RATE  VERSUS  INDUCTANCE 

In  measuring  the  maximum  uniform  glow  repetition  rate  at  various 
energy  densities  it  is  clear  that  the  repetition  rate  decreases  as  energy 
density  increases.  Since  the  energy  density  was  varied  by  increasing  the 
size  of  the  storage  capacitor  the  pulse  length  also  increased  because 
pulse  length  is  proportional  to  /LC  .  The  purpose  of  this  experiment  was 
to  see  if  the  maximum  repetition  rate  that  could  be  achieved  for  a 
uniform  discharge  was  a  function  of  pulse  length  as  well  as  energy  density. 

The  experiment  was  run  in  60  torr  N2 ,  80  torr  C02 ,  and  700  torr  He. 

For  each  gas  the  charging  voltage  was  kept  constant  at  10  KV  and  the 
storage  capacitor,  C,  was  kept  constant  at  either  40.8  nf  or  81.2  nf .  The 
pulse  length  was  changed  by  varying  L,  the  circuit  inductance  (Fig.  6.24). 
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6  0  T  o  r  r  N  2 


Figure  6.22  Observed  Fractional  Frequencies  for  a  Uniform  Glow  in  60  Torr 
N« .  The  solid  line  is  the  predicted  fractional  freauencies 
using  Dzakowic  and  Wutzke's  theory. 
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80  To  r  r  CO2 


•  arc  between  electrodes 


Figure  6.23  Observed  Fractional  Frequencies  for  a  Uniform  Glow  in 

80  Torr  CO^.  The  solid  line  is  the  predicted  fractional 

frequency  using  Dzakowic  and  Wutzke’s  theory.  The  two 

sets  of  data  points  correspond  to  two  different  arcing  criteria. 


L 


Fig.  6.24  Discharge  Circuit  For  Pulse  Length  Experiments. 

For  a  given  value  of  C  the  pulse  length  was  varied 


by  changing  L. 
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Xn  this  way  the  pulse  length  varied  but  the  energy  per  pulse  remained 
constant.  L  was  varied  from  2  pH,  the  intrinsic  circuit  inductance, 
to  101  pH  producing  pulse  lengths  that  varied  from  2  to  6  psec.  For 
each  value  of  L  maximum  repetition  rates  for  a  uniform  discharge  were 
measured  for  seven  flow  velocities  and  an  average  fractional  frequency 
calculated  as  outlined  in  section  6.5. 

The  results  of  this  experiment  are  summarized  in  the  plots  of 
Figures  6.25  to  6.27.  They  show  that  the  fractional  frequency  attainable 
is  essentially  independent  of  inductance,  or  pulse  length.  The  two  sets 
of  fractional  frequencies  shown  for  40.8  nf  in  80  torr  CO^  require  some 
explaining.  At  an  energy  density  of  about  35  J/£  two  types  of  repetition 
rate  arcing  could  be  seen.  The  lower  fractional  frequency  of  about  .57 
marked  the  transition  region  at  which  an  arc  between  the  cathode  and 
spark  pin  appeared.  Increasing  the  repetition  rate  further  shifts  the 
arc  to  between  the  downstream  edges  of  the  two  electrodes.  This  transition 
occurred  at  fractional  frequencies  of  roughly  .77.  Figure  6.23  also 
notes  this  phenomena.  The  various  arcing  modes  are  discussed  in  more 
detail  in  section  6.10.  Figure  6.26  shows  that  the  onset  of  both  arcing 
modes  is  independent  of  pulse  length. 

6.7  AFTERGLOW  CLEARING  TIME  MEASUREMENTS 

Luminous  afterglows  were  observed  in  N^ ,  CC^,  and  He.  This  made  it 
possible  to  use  a  photomultiplier  tube  as  a  diagnostic  for  studying  the 
convection  of  discharge  products  away  from  the  discharge  region. 

The  apparatus  for  this  study  is  shown  in  Figure  6.28.  A  mirror, 
pinhold,  telescope  and  photomultiplier  were  mounted  on  an  optical  bench. 

The  pinhole  and  telescope  acted  as  a  collimator  to  limit  the  photo¬ 
multiplier’s  field  of  view  to  a  very  narrow  beam  in  space.  By  translating 
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60  Torr  N2 


Fig. 


6.25  Repetition  Rate  as  a 
N^-  The  repetition 


Function  of  Circuit  Inductance  for  60  Torr 
rate  has  been  normalized  by  expressing 


it  as  a  fractional  frequency. 
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Fig.  6.26  Repetition  Rate  as  a  Function  of  Circuit  Inductance  for  80  Torr 
CO^.  The  sets  of  data  at  about  .77  and  . 
represent  two  different  arcing  criteria. 
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Fig.  6.27  Repetition  Rate  as  a  Function  of  Circuit  Inductance  for  700 


Torr  He. 
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the  mirror  along  the  optical  bench  the  afterglow  could  be  examined  at 
various  distances  downstream  from  the  discharge.  The  signal  from  the 
photomultiplier  went  to  the  scope  through  a  short  cable  (0.5  m)  and  a 
100  K£2  pad.  The  short  cable  and  pad  insured  that  the  rc  time  constant 
was  short  enough  to  resolve  the  waveform  and  the  100  Kfi  pad  was  the 
minimum  value  that  produced  a  signal  voltage  large  enough  to  see.  A 
bright  orange  afterglow  could  be  produced  in  60  torr  using  an  80  nf 
capacitor  charged  to  10  Kv.  Similar  but  fainter  afterglows  were  produced 
in  CO^  and  He. 

A  typical  afterglow  pulse  is  shown  in  Figure  6.29.  The  initial 
spike  is  due  to  the  light  flash  when  the  discharge  occurs.  This  discharge 
leaves  a  plug  of  gas  in  an  excited  metastable  state.  The  plug  of  gas 
blowing  downstream  past  the  position  probed  by  the  photomultiplier 
produces  the  square  pulse.  The  dip  between  the  spike  and  the  front  edge 
of  the  pulse  corresponds  to  the  time  needed  for  the  downstream  edge  of 
the  excited  gas  to  reach  this  position.  The  back  edge  of  the  pulse 
corresponds  to  the  time  needed  for  the  upstream  edge  of  the  excited  gas 
to  pass  this  position.  So  the  pulse  length  is  proportional  to  the 
diameter  of  the  discharge,  the  two  being  related  by  the  gas  flow  velocity. 
Ideally  the  pulse  should  be  flat  topped  but  the  droop  corresponds  to  a 
reduction  of  intensity  with  time  as  the  metastable  population  decreases. 

In  principle  the  recombination  coefficient  for  the  metastable  could  be 
estimated  from  the  pulse  decay  using: 

-  =  —  +  6t 

n  n 

o 

but  in  order  to  do  this  the  initial  metastable  density,  n^ ,  must  be  known. 

The  CC>2  afterglow  pulse  is  similar  to  the  afterglow  pulse  except 


’ 


0.5  ms  e  c/d  iv 


Fig.  6.29  Typical  Afterglow  Waveform  in  60  Torr  •  The  discharge 
capacitor  was  81.2  nf  and  the  charging  voltage  was  10  KV . 
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that  it  is  lower  in  amplitude  and  has  a  flatter  top.  The  intensity 
of  the  CO^  afterglow  increases  after  the  discharge  has  run  for  a  while. 
Since  the  discharge  dissociates  the  CO^  into  CO  and  0^  this  suggests  that 
the  afterglow  in  CO^  is  associated  with  the  presence  of  one  of  these 
products . 

The  He  afterglow  proved  to  be  difficult  to  reproduce.  Very  pure  He 
(99.999%)  was  used  in  an  effort  to  try  and  improve  the  afterglow.  But 
the  afterglow  actually  got  fainter  with  repeated  fillings  of  the  very 
pure  gas.  Adding  a  fraction  of  a  torr  of  to  the  He  made  the  afterglow 
come  back  strongly.  So  the  afterglow  observed  in  He  was  due  to  the 
presence  of  impurities. 

Using  an  81.2  nf  capacitor  charged  to  10  Kv , afterglows  were  produced 
in  the  three  gases.  The  afterglows  were  observed  at  a  point  V  upstream 
from  the  back  edge  of  the  electrodes.  The  time  from  the  beginning  of 
the  discharge  spike  to  the  end  of  the  afterglow  pulse  was  observed  as  a 
function  of  gas  flow  velocity.  This  time  could  be  taken  as  a  reasonable 
estimate  of  the  time  needed  to  clear  the  discharge  products  from  the 
discharge  region.  The  inverse  of  the  afterglow  pulse  lengths  are  plotted 
as  a  function  of  gas  flow  velocity  in  Figures  6.30  to  6.32.  For  compar¬ 
ison  the  ideal  repetition  rate  is  shown  as  a  dotted  line. 

In  both  and  CO^  the  inverse  of  the  time  from  the  beginning  of  the 
spike  to  the  end  of  the  pulse  was  proportional  to  the  gas  flow  velocity 
as  expected.  In  He  however  the  afterglow  pulse  length  tended  to  stay 
quite  constant  indicating  that  recombination  occurred  before  all  the  gas 
was  convected  out  from  between  the  electrodes.  As  was  already  mentioned 
He  tends  to  form  a  nonuniform  discharge  rather  than  arcs  if  fideal  is 
exceeded.  If  the  repetition  rate  is  increased  further  arcs  will  form  in 
the  nonuniform  glow.  These  afterglow  measurements  suggest  that  the  high 
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Fie.  6.30  1/x  „  ,  as  a  Function  of  Flow  Velocity. 
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Data  for  60  torr  N2  and  58  J/£  energy  density. 
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FLOW  VEL  M/SEC 


Fig.  6.31  1/t  ,  as  a  Function  of  Flow  Velocity. 

6  afterglow 

Data  for  80  torr  C02  and  70  Jfl  energy  density. 
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Data  for  700  torr  He  and  37  J / £  energy  density. 
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repetition  rates  that  can  be  achieved  in  pure  He  before  arcing  occurs  may 
be  due  to  fast  recombination  of  left  over  ionization  from  previous 
discharge  pulses.  As  a  result  it  would  not  be  necessary  for  the  electrodes 
to  be  cleared  by  convection  before  the  next  discharge  pulse. 

So  in  and  CO^  the  afterglow  lifetime  is  longer  than  the  convect¬ 
ion  time  and  observing  the  afterglow  can  give  some  idea  of  the  time 
needed  to  clear  discharge  products  from  between  the  electrodes.  In  He 
however  the  afterglow  lifetime  is  shorter  than  the  convection  time  so 
the  afterglow  measurements  cannot  provide  any  useful  measure  of  clearing 
time.  Even  in  CO^  and  N  although  qualitatively  the  inverse  of  the 
afterglow  clearing  time  looks  like  a  reasonable  estimate  of  obtainable 
repetition  rate  the  quantitative  agreement  is  poor,  as  can  be  seen  from 
Table  6.3. 


gas  and  energy  density 

— 

fractional 

frequency  based 

on 

1 

Taf terglow 

— 

fractional 

frequency 

observed 

80  torr  CO^  70  J/l 

60  torr  N^  58  J / i 

.67 

| 

.98 

.46 

| 

j 

Table  6 . 3 

6.8  AFTERGLOW  PULSE  LENGTH  VERSUS  ENERGY  DENSITY 

The  same  apparatus  described  in  section  6.7  was  used  for  this 
experiment.  Since  maximum  repetition  rate  decreases  with  increasing 
energy  (section  6.5),  and  Dzakowic  and  Wutzke  suggested  that  this  was  due 
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to  discharge  expansion,  this  experiment  attempted  to  observe  the  expansion 

by  monitoring  the  afterglow  pulse  length.  The  afterglow  pulse  length  is 

t  -t  in  Fig.  6.33  and  should  be  proportional  to  the  diameter  of  the  discharge 
D  A 

products  blowing  downstream,  the  two  being  related  by: 


diameter  =  v  ( t  -t „ 
g  B  A 


) 


where  v  is  the  gas  flow  velocity. 

§ 

The  gas  flow  velocity  was  kept  constant  at  19.7  m/sec  in  N9  and 
22.5  m/sec  in  CO^ .  The  discharge  diameter  was  measured  with  a  travelling 
telescope.  It  was  58  mm  in  and  varied  from  46  to  58  mm  in  CO^.  The 
discharge  energy  density  was  varied  by  changing  the  value  of  storage 
capacitor  and  the  values  t  and  t  were  measured  as  a  function  of  increasing 
energy  density. 


Figure  6.33  Afterglow  pulse 
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The  afterglow  pulse  lengths,  t  -t  ,  are  plotted  as  a  function  of 

Jj  A 

energy  density  in  Figures  6.34  and  6.35.  The  pulse  length  increases  with 
increasing  energy  density.  A  prediction  for  the  pulse  lengths  can  be 
obtained  from  equation  6.1.6.  It  says  that: 


r2  =/N/(l+e)1/y  r± 


and  so  the  expanded  diameter  is  related  to  the  initial  diameter  by 


D2  =^(l+e)1/y  D 


where  D  is  the  initial  diameter  and  D2  is  the  final  diameter.  The  pulse 
length  is  then  given  by: 


At  =  -  CA  = 


V(l+£)1/ 


V 


D 


v 


6.8.1 


Some  interesting  observations  arise  from  applying  equation  6.8.1 
to  the  data.  First  the  afterglow  pulse  is  due  to  a  glowing  plug  of  gas 
that  has  a  different  diameter  than  the  discharge  diameter.  The  minimum 
pulse  length  is  given  by: 


At 


D_ 

v 


6.8.2 


since 


V(l+£)1/Y 


1  as  e  0.  Using  the  observed  discharge  diameter  of 


58  mm  and  velocity  of  19.7  mm/msec  for  N2  gives  a  minimum  pulse  length, 
according  to  equation  6.8.2,  of  2.94  msec.  But  all  the  observed  pulse 
lengths  in  N0  are  less  than  this  value.  Similarly  in  CO  the  discharge 
diameters  of  46  to  58  mm  and  flow  velocity  of  22.5  mm/msec  give  a  minimum 


.  -  w  ' 


E  i  ||  V  + :  ^  e:i,Tti 


msec 
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Fig.  6.34  Afterglow  Pulse  Length  as  a  Function  of  Energy  Density  for  60 

/1+  n 1/y 

Torr  .  The  predicted  curve  is  given  by  At=  - - - 


where  W=31.7  mm  and  v=19.7  mm/msec. 
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80  Tor  r  C02 


Fig.  6.35  Afterglow  Pulse  Length  as  a  Function  of  Energy  Density  for 


80  Torr  C02.  The  predicted  curve  is  given  by  At= 


vfl+e)1/Y 


where  W=41.2  mm  and  v=22.5  mm/msec. 


v 
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pulse  length  of  between  2.04  and  2.58  msec.  Values  both  lower  and 
higher  than  this  range  were  observed.  The  fact  that  higher  values  were 
observed  means  that  some  discharge  expansion  must  have  occurred,  and 
because  lower  values  were  observed,  the  initial  diameter  of  the  afterglow 
plug  must  have  been  smaller  than  the  diameter  observed  with  the  telescope. 


These  results  suggest  that  the  afterglow  pulse  does  not  represent  the 

entire  discharge  volume  but  rather  a  more  strongly  heated  inner  core. 

From  the  data  a  diameter,  W,  can  be  estimated  for  this  core  and  then  its 

Wv{l+e)1/Y 


expansion  predicted  using  the  formula  At 


v 


to  give  the  pulse 


length  due  to  the  expanded  gas.  A  comparison  of  this  model  with  the 
observed  pulse  lengths  (Figures  6.34  and  6.35)  shows  that  it  fits  the 
discharge  well  but  does  not  seem  to  fit  the  CO^  discharge  very  well.  The 
predicted  slope  in  CO^  is  less  steep  than  that  observed  for  the  data. 

There  is,  however,  a  complication  with  the  data  for  CO^  in  that 
there  was  a  large  variation  in  the  measured  initial  discharge  diameters, 
and  in  general  the  initial  discharge  diameter  tends  to  increase  with 
increasing  energy  density.  The  increase  in  initial  diameter  would 
probably  be  accompanied  by  a  corresponding  increase  in  core  diameter  which 
could  account  for  the  discrepancy.  This  can  be  compensated  for  by  calculat¬ 
ing  W,  the  afterglow  diameter,  from  the  observed  pulse  length,  At,  and 

W 

then  plotting  —  where  D  is  the  discharge  diameter.  If  the  discrepancy  is 

W 

due  to  the  increase  of  D  with  energy  then  —  should  stay  constant.  Table 

6.4  shows  the  result  of  this  calculation  for  and  Table  6.5  for  CO^. 

W 

The  results  are  plotted  in  Figure  6.36.  shows  a  constant  —with  any 

variation  being  due  to  random  error.  Although  the  discrepancy  between 


the  data  and  the  model  in  CO seems  reduced  there  is  still  an  increase 
in  core  diameter  with  increasing  energy  density  that  the  model  does  not 
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60  torr 

C  E/vol 

nf  J/£ 

gas  velocity=  19 

v(i+£)1/y  At 

msec 

.  7  m/ sec 

W=  vA  t 

Va+e)1/Y 

D=  58.42  mm 

W/D 

20.4 

16 

1.22 

2.0 

32.30 

.55 

30.3 

24 

1.30 

2.2 

33.34 

.57 

40.8 

32 

1.38 

2.2 

31.41 

.54 

61.2 

48 

1.51 

2.4 

31.31 

.54 

81.2 

64 

1.63 

2.5 

30.21 

.52 

20.4 

16 

1.22 

1.95 

31.49 

.54 

30.7 

24 

1.30 

2.0 

30.31 

.52 

40.8 

32 

1.38 

2.25 

32.12 

.55 

51.1 

40 

1.45 

2.35 

31.93 

.55 

61.2 

48 

1.51 

2.4 

31.31 

.54 

71.1 

56 

1.57 

2.45 

30.74 

.53 

81.2 

63 

1.62 

2.45 

29.79 

.51 

81.2 

63 

1.62 

2.7 

32.83 

.  56 

81.2 

63 

1.62 

2.65 

32.23 

.55 

71.2 

56 

1.57 

2.55 

32.00 

.55 

61.2 

48 

1.51 

2.5 

32.62 

.56 

51.1 

40 

1.45 

2.45 

33.29 

.57 

40.8 

32 

1.38 

2.35 

33.55 

.57 

30.3 

24 

1.30 

2.2 

33.34 

.57 

20.4 

16 

1.22 

2.0 

32.30 

.55 

W/D  average  =  .55 

Table  6.4  Further  analysis  of  afterglow  data  for  60  torr  N2 . 
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80  torr  CO^  gas  velocity  =  22.5  m/sec 


c 

nf 

D 

mm 

E/ vol 
j/  a 

vfl+e)1/Y 

At 

msec 

W=vAt 

W/D 

a/i+c)1/y 

40.8 

50.80 

42 

1.32 

2.15 

36.65 

.72 

61.2 

50.80 

63 

1.44 

2.6 

'  40.63 

.80 

81.2 

53.34 

76 

1.51 

2.9 

43.21 

.81 

81.2 

58.42 

64 

1.45 

2.9 

45.00 

.77 

81.2 

58.42 

64 

1.45 

2.95 

45.78 

.78 

71.1 

53.34 

67 

1.46 

2.7 

41.61 

.78 

61.2 

53.34 

58 

1.41 

2.55 

40.69 

.76 

51.1 

53.34 

48 

1.36 

2.45 

40.53 

.76 

40.8 

50.80 

42 

1.32 

2.3 

39.20 

.77 

30.3 

50.80 

31 

1.25 

2.0 

36.00 

.71 

20.4 

45.72 

26 

1.21 

1.75 

32.54 

.71 

81.2 

55.88 

70 

1.48 

2.95 

44.85 

.80 

W/D  average  =  .76 


Table  6.5  Further  analysis  of  afterglow  data  for  80  torr  CO^. 
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Energy  Density  J/L 


Fig.  6.36  Agreement  Between  Observed  Afterglow  Diameters  and  Gas 
Expansion  Model  for  N.-,  and  CO^ . 
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account  for. 

This  experiment  demonstrated  the  existence  of  long  lived  ionization 
products  which  blow  downstream  as  a  plug.  The  increase  in  plug  diameter 
with  increasing  energy  density  agrees  fairly  well  with  the  expansion 
predicted  by  Dzakowic  and  Wutzke  and  so  supports  their  claim  that  gas 
expansion  is  an  effect  which  lowers  the  achievable  repetition  rate. 
However,  another  effect  that  is  probably  more  important  becomes  evident. 
It  is  a  problem  with  system  geometry.  With  the  spark  pins  downstream 
from  the  discharge  the  plug  of  discharge  products  can  bridge  the  gap 
between  the  electrodes  and  the  spark  pins.  If  its  ionization  density 
is  high  enough  an  arc  will  form  from  the  cathode  to  the  spark  pin.  This 
is  the  kind  of  arc  that  usually  formed  in  and  formed  in  the  other 
gases  as  well  at  high  enough  energy  densities. 

An  additional  observation  in  the  flowing  gas  discharge  system 
supports  this  ionizing  bridge  concept.  Arcing  occurred  when  the  spacing 
between  afterglow  pulses  was  approximately  2  msec  or  less  in  60  torr  N^. 


Figure  6.37  Arcing  occurs  when  the  spacing  between  afterglow 
pulses  was  2  msec  or  less  in  60  torr  N^. 
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A  fractional  frequency  can  be  calculated  that  is  based  on  the 
afterglow  clearing  time,  t^,  plus  2  msec.  In  60  torr  the  discharge 
diameter  is  58  mm  and  the  gas  flow  velocity  is  19.7  mm/msec  which  gives 
an  ideal  repetition  rate  of  339  Hz.  Table  6.6  shows  the  fractional 
frequencies  calculated  from  the  afterglows. 


r 

E 

CB 

msec 

1 

^ f ract 

nf 

Vol 

J/£ 

t  +2msec 

D 

Hz 

20.4 

16 

2.8 

208 

.61 

j  30.3 
* 

24 

2.9 

204 

.60 

oo 

• 

o 

32 

2.9 

204 

.60 

j  61.2 

48 

3.1 

196 

.58 

f  81.2 

l 

64 

3.2 

192 

.57 

| 

Table  6.6  Fractional  frequencies  calculated  from  afterglows. 

The  fractional  frequencies  calculated  from  the  afterglow  clearing 
times  are  plotted  in  Figure  6.38  along  with  the  fractional  frequencies 
obtained  by  observing  the  discharge.  The  two  agree  well.  Two  msec 
corresponds  to  39  mm  at  a  flow  velocity  of  19.7  mm/msec  and  the  spark 
pins  were  38  mm  downstream  from  the  position  probed  by  the  photomultiplier. 
So  in  60  torr  the  expanded  discharge  products  must  clear  not  only 
the  electrodes  but  the  spark  pins  as  well. 

The  bridging  effect  of  the  discharge  products  also  explains  an 
arcing  problem  that  occurred  in  the  laser  system.  High  repetition  rate 
arcing  would  often  occur  from  the  cathode  to  the  Marx  generator  spark 
pins  and  would  burn  a  carbon  track  on  the  perspex  insulation  between 
the  two.  The  carbon  track  would  then  have  to  be  cleaned  off  before  the 
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Figure  6.38 


Comparison  Between  Observed  Fractional  Frequencies  for  60 
Torr  N9  and  Fractional  Frequencies  Calculated  From  the 


Afterglow  Clearing  Times. 
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system  could  operate  again.  In  the  flowing  gas  discharge  system  the 
tracking  problem  was  greatly  reduced  by  replacing  the  perspex  with 
ceramic  sheet  which  is  much  more  heat  resistant. 

6.9  DISCHARGE  INTENSITY  PROFILES 

In  the  repetition  rate  and  afterglow  diameter  versus  energy  density 
measurements  that  were  discussed  in  sections  6.5  and  6.8,  the  energy 
density  was  estimated  using  a  discharge  volume  obtained  by  observing 
the  edges  of  the  discharge  with  a  travelling  telescope.  To  verify  that 
this  gave  a  reasonable  estimate  of  discharge  diameter,  discharge  diameters 
were  measured  with  the  travelling  telescope  and  intensity  profiles  of 
the  discharges  were  measured  with  the  photomultiplier  for  comparison. 

The  apparatus  used  was  the  same  as  for  the  afterglow  measurements  (Figure 
6.31)  except  that  the  photomultiplier  was  terminated  with  a  50^  pad  to 
give  time  resolution  of  the  discharge  pulse.  By  moving  the  plane  mirror 
the  discharge  was  scanned  and  peak  intensities  recorded  at  5  mm  intervals 
to  develop  intensity  profiles  of  the  discharge.  Typical  intensity 
profiles  for  uniform  discharges  are  plotted  in  Figures  6.39  to  6.44. 

The  plots  are  normalized  to  show  a  peak  value  of  1.  The  diameters 
measured  with  the  travelling  telescope  and  the  diameters  representing 
10%  of  peak  intensity  on  the  intensity  profiles  agree  to  within  4  mm  or 
less  for  discharges  of  typically  55  to  70  mm  in  diameter.  So  the  diameter 
observed  with  the  travelling  telescope  is  a  good  indicator  of  discharge 
volume . 

The  intensity  profiles  could  also  be  expanded  into  a  diagnostic 
to  study  discharge  uniformity.  For  example  intensity  profiles  for  dis¬ 
charges  which  became  nonuniform  due  to  excessive  repetition  rates  are 
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700  TORR  HE 
81 .2  NF 
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Fig.  6.39  Intensity  Profile  for  700  Torr  He  Discharge  with 


81.2nf  Capacitor. 
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700  TORR  HE 
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Fig.  6.40  Intensity  Profile  for  700  Torr  He  Discharge  with 


20 . 4nf  Capacitor. 
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60  TORR  112 
81 .2  NF 


Fig.  6.41  Intensity  Profile  for  60  Torr  Discharge  with 


81.2nf  Capacitor. 
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60  TORR  N2 
20.4  NF 


POSITION  CM 

Fig.  6.42  Intensity  Profile  for  60  Torr  Discharge  with 


20.4nf  Capacitor. 
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80  TORR  C02 
81  .2  NF 


POSITION  CM 


Fig.  6.43  Intensity  Profile  for  80  Torr  C02  Discharge  with  • 
81.2nf  Capacitor. 
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80  TORR  C02 
20.4  NF 


POSITION  CM 


Fig.  6.44  Intensity  Profile  for  80  Torr  CO?  Discharge 


with  20.4nf  Capacitor. 


shown  in  Figures  6.45  and  6.46.  They  agree  with  the  qualitative 
variations  visible  to  the  eye.  The  He  discharge  was  brighter  on  the 


175 


downstream  edge  and  the  N?  discharge  on  the  upstream  edge. 


Although  it  is  difficult  to  determine  the  light  intensity  distribution 
from  an  intensity  profile  it  is  easy  to  calculate  the  intensity  profile 
from  a  given  light  intensity  distribution.  For  comparison  with  the 
observed  intensity  profiles  a  profile  is  calculated  for  a  completely  uniform 
cylindrical  discharge.  Assume  the  discharge  has  a  circular  cross  section 
of  radius  1.  If  it  is  completely  uniform,  the  light  intensity  observed 
will  be  proportional  to  the  area  probed  by  the  photomultiplier.  Take  the 
width  probed  by  the  photomultiplier  to  be  w.  The  geometry  of  the  problem 
is  shown  in  Figure  6.47.  Due  to  symmetry  the  area  probed  in  one  quadrant 
only  needs  to  be  found.  The  total  area  will  be  twice  this  amount. 


so 


In  the  top  right  hand  quadrant  an  element  of  area  is  given  by: 


dA  = 


ydx  =  v  1 


-  x  dx 


so  the  area  is: 


2  ^ 

(1  -  x  )  2  dx 


A  =  / 


w 


w 


2 


2 
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700  TORR  HE 
20.4  NF 


POSITION  CM 


Fig.  6.45  Intensity  Profile  for  a  Repetition  Rate  Induced  Nonuniform 
Discharge . 

The  discharge  is  in  700  torr  He  with  a  20.4nf  capacitor. 
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60  TORR  N2 
20.4  NF 


POSITION  CM 

Fig.  6.46  Intensity  Profile  for  a  Repetition  Rate  Induced  Nonuniform 
Discharge . 

The  discharge  is  in  60  torr  with  a  20.4nf  capacitor. 
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Fig .  6.47 


Geometry  of  Uniform  Discharge  Model  used  for 


Calculating  a  Model  Intensity  Profile. 
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From  symmetry  the  area  both  above  and  below  the  x  axis  is  double  this 
so  the  total  area  is  given  by: 

Total  Area  =  (x  +  y^-Cx  +  |)  2  -  (x  -  y^-Cx  -  |)2 

+  sin  1(x  +  y)  -  sin  1(x  -  ^)  6.9.1 

A  modification  of  this  formula  is  useful  near  the  boundary  x  =  1.  In 
this  case  the  integration  limits  are  changed  to  x-w  and  x.  Then  the 
formula  becomes: 

x  o  u 

A  =  /  ( 1  —  x  )  2  dx 

x-w 


-  [ 


Jx  -  2  '  -  x 


xVl  -  x  ,1  .  -1  -i 

— -  +  —  sm  xj 


x-w 


Again,  by  symmetry  the  total  area  is  double  this  so: 


Total  Area 


=  ■  -  -1 


-  -  /  2  -1 
x  +  sin  x  -  (x-w)yl-(x  -w)  -  sin  (x-w) 


6.9.2 


The  total  intensity  is  just  some  constant  times  the  area  so  6.9.1 
and  6.9.2  can  be  used  to  plot  a  representitive  intensity  profile  for  a 
uniform  discharge.  The  discharge  radius  was  typically  35  mm,  the  scan 
spacing  was  5  mm,  and  the  resolution  (w)  was  7  mm.  Normalizing  this  the 
radius  becomes  1,  the  scan  spacing  becomes  .143  and  w  becomes  .20.  The 
intensity  profile  is  shown  in  Figure  6.48. 
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Fig.  6.48  Model  Intensity  Profile  for  a  Uniform  Discharge. 
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The  procedure  can  be  generalized  to  predict  the  intensity  profile 
for  nonuniform  intensity  distributions.  In  that  case  the  intensity 
profile  is  given  by: 


1  (x)  -  /  w  /  / - j  dy  dx 

x  -  —  -  x 

where  f(x,y)  is  the  intensity  distribution  function.  In  the  case  of  a 
uniform  discharge  f(x,y)  =  1  and  the  results  simplify  to  equation  6.9.1. 

6.10  ARCING  MODES 

Two  distinct  types  of  discharge  instability  were  observed  in  the 
flowing  gas  discharge  system.  The  first  is  related  to  energy  density. 

When  the  energy  density  becomes  too  high  a  glow  to  arc  transition  occurs. 
This  transition  occurs  for  an  increasing  fraction  of  the  discharge 
pulses  as  energy  density  rises,  until  at  high  enough  energy  densities  it 
occurs  all  the  time.  At  the  other  end  of  the  scale,  for  very  low  energy 
densities  it  does  not  occur  at  all.  This  transition  occurs  independently 
of  repetition  rate  and  between  the  plane  portion  of  the  electrodes. 

The  other  type  of  discharge  instability  is  repetition  rate  induced 
and  is  related  to  gas  flow  velocity.  It  is  due  to  not  sufficiently 
clearing  the  discharge  products  from  the  previous  pulse  out  of  the  discharge 
region.  This  instability  can  appear  in  one  of  three  typical  modes  which 
are  illustrated  in  Figure  6.49.  They  are:  pulling,  arcing  between  the 
electrodes  on  the  downstream  edge,  and  arcing  from  the  cathode  to  the 
spark  pins.  Although  all  three  repetition  rate  instabilities  can  be  made 
to  occur  in  each  of  the  gases  certain  gases  tend  to  favour  certain 

Pulling  occurs  in  He  for  most  energy  densities  and  in  N* 


instabilities . 
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Figure  6.49  Repetition  Rate  Induced  Arcing  Modes,  (a)  Is  a  uniform 
discharge,  (b)  Shows  discharge  pulling,  (c)  Shows  an 
auxiliary  discharge  (observed  in  C0?),  (d)  Shows  an  arc 
between  the  electrodes,  and  (e)  Shows  an  arc  from  the  cathode 


to  a  spark  pin. 
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at  low  energy  densities.  A  variation  of  pulling  occurs  in  CC^  at  low 
energy  densities.  In  this  variation  a  second  smaller  auxiliary  discharge 
forms  on  the  downstream  portion  of  the  electrodes.  At  higher  energy 
densities  arcing  occurs  and  can  be  either  between  the  electrodes  or 
from  the  cathode  to  the  downstream  spark  pin.  Arcs  in  He  and  CO^  tend 
to  occur  between  the  electrodes  whereas  arcs  in  tend  to  occur  from  the 
cathode  to  the  spark  pin.  The  key  difference  between  the  repetition  rate 
induced  instabilities  and  the  energy  induced  instability  is  that  the 
repetition  rate  effects  occur  on  the  downstream  edge  of  the  discharge  while 
the  energy  instability  occurs  toward  the  center  of  the  discharge  where 
the  electrodes  are  flat. 

Figures  6.50  and  6.51  illustrate  the  typical  repetition  rate 
instabilities  encountered,  and  the  repetition  rates  at  which  they  occur. 

The  instability  repetition  rates  were  measured  by  increasing  the  repetition 
rate  until  a  particular  instability  appeared  to  occur  continuously.  The 
lines  disch,  elect,  and  pins  are  estimated  repetition  rates  based  on  the 
convection  times  needed  for  gas  to  clear  the  discharge  region,  electrodes, 
and  spark  pins  respectively.  These  estimates  ignore  possible  effects 
due  to  discharge  expansion.  The  dominant  instability  is  the  instability 
that  occurs  at  the  lowest  repetition  rate. 

In  He  the  dominant  instability  is  pulling  and  it  occurs  if  the 
discharge  gas  has  not  had  time  to  clear  the  electrodes  between  pulses. 

If  the  repetition  rate  is  pushed  higher  until  the  gas  has  not  yet  had  time 
to  clear  the  discharge  region  arcing  occurs  between  the  electrodes. 

Finally  if  the  repetition  rate  is  pushed  higher  still  arcing  between  the 
cathode  and  spark  pin  will  occur.  In  60  torr  the  dominant  instability 
is  pulling  for  energy  densities  less  than  approximately  45  J /£.  If  the 
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Fig.  6.50  Repetition  Rate  Induced  Instabilities  as  a  Function  of  Energy 

700  Torr  He  for  a  Gas  Flow  Velocity  of  24  m/sec. 
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Fig.  6.51  Repetition  Rate  Induced  Instabilities  as  a  Function  of  Energy 
Density  in  60  Torr  N2  for  a  Gas  Flow  Velocity  of  19.5  m/sec. 


repetition  rate  is  increased  arcing  between  the  cathode  and  the  spark  pin 
sets  in.  For  still  higher  repetition  rates  arcing  between  the  electrodes 
also  occurs.  Note  that  for  energy  densities  above  45  J/£  pulling  stops 
and  arcing  between  the  cathode  and  spark  pin  becomes  the  dominant  mode. 

Figures  6.52  to  6.56  show  the  basic  arcing  modes  and  extend  the 
measurements  to  higher  energy  densities.  The  vertical  bars  map  the 
regions  in  which  energy  induced  instabilities  occur  ranging  from  no  energy 
induced  arcing,  through  sporadic  arcing,  to  a  completely  unstable  discharge. 
Note  that  in  80  torr  CO^  there  is  a  transition  in  arcing  modes  from  arcing 
between  the  electrodes  to  arcing  from  the  cathode  to  the  spark  pin  in  the 
energy  density  range  of  30  to  40  J/i.  The  same  transition  occurs  in  the 
400  torr  laser  mix  at  about  65  J/£. 

The  repetition  rate  instabilities  observed  and  the  repetition  rates 
for  which  they  occur  are  consistent  with  gas  convection  not  having  had 
enough  time  to  sweep  the  discharge  products  out  of  the  regions  between  the 
electrodes  or  the  spark  pins.  In  particular  the  transition  in  arcing  modes 
in  CO^  and  in  the  400  torr  laser  mixture  seems  to  support  this  view 
because  it  occurs  at  a  repetition  rate  for  which  the  gas  has  cleared  the 
electrodes  but  not  the  space  between  the  electrodes  and  the  spark  pins. 

Hence  arcing  jumps  from  being  between  the  electrodes  to  being  between  the 
cathode  and  a  spark  pin.  In  CO^  at  higher  energy  densities  the  repetition 
rate  drops  below  that  predicted  from  the  convection  time  for  clearing  the 
spark  pins.  But  this  can  be  explained  by  expansion  of  the  discharge 
gases  which  will  increase  the  time  needed  to  clear  the  spark  pins  to  somewhat 
longer  than  predicted  from  the  convection  time. 
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700  Torr  He  for  a  Gas  Flow  Velocity  of  24.7  m/sec. 
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Fig.  6.53  Onset  of  Arcing  Modes  as  a  Function  of  Energy  Density  in 
60  Torr  for  a  Gas  Flow  Velocity  of  19.5  m/sec. 
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Fig.  6.55  Onset  of  Arcing  as  a  Function  of  Energy  Density  in  400  Torr 
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Gas  flow  velocity  is  24.8  m/sec. 
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Fig.  6.56  Onset  of  Arcing  as  a  Function  of  Energy  Density  in 
500  Torr  of  1:1:8  C02:N2:He  Laser  Mix. 

Gas  flow  velocity  is  24.5  m/sec. 
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6.11  GAS  POISONING 

Besides  repetition  rate  effects  an  additional  problem  was  observed 
in  the  laser  gas  mixture  that  was  not  apparent  in  the  individual  gases. 

The  problem  is  gas  poisoning. 

A  discharge  was  run  in  a  1:1:8  mixture  of  N^^CO^^He  at  400  torr 
pressure  and  100  Hz  repetition  rate.  The  energy  per  pulse  was  4.6  J. 
Although  the  discharge  was  initially  uniform,  sporadic  arcing  occurred 
after  about  four  minutes.  A  similar  result  was  observed  for  a  500  torr 
mixture.  The  poisoning  of  the  laser  mixture  is  almost  certainly  due  to 
the  dissociation  of  CO^  into  CO  and  0^  as  demonstrated  by  Pace  and  Lacombe 
(22) .  This  poisoning  was  also  probably  an  important  but  unrecognized 
problem  in  the  high  repetition  rate  laser  system. 

Although  dissociation  of  CO^  was  observed  for  discharges  in  80  torr 
CO^  in  the  flowing  gas  discharge  system,  the  effect  on  discharge  uniformity 
seemed  negligible.  This  was  probably  because  the  CO^  discharge  was 
running  at  a  fairly  low  pressure  and  so  was  less  likely  to  collapse  into 
an  arc. 

6.12  CONCLUSIONS 

The  results  of  this  study  show  that  the  maximum  repetition  rate 
limit  for  a  uniform  discharge  in  a  flowing  gas  high  repetition  rate  system 
is  due  essentially  to  the  time  needed  to  clear  the  discharge  region  by 
gas  convection.  The  clearing  time  can  be  complicated  by  several  things. 
Expansion  of  the  discharge  products  reduces  the  attainable  repetition  rate. 
Also  system  geometry  can  reduce  the  repetition  rate.  For  example  in  the 
flowing  gas  discharge  system  (and  the  high  repetition  rate  laser) ,  when 
the  energy  density  gets  large  enough  the  discharge  products  must  clear 
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not  only  the  electrodes  but  the  spark  pins  as  well  resulting  in  a  large 
decrease  in  attainable  repetition  rate. 

In  this  experiment  the  effect  of  convective  clearing  of  discharge 
products  on  repetition  rate  is  somewhat  obscured  by  the  transitions  in 
arcing  mode  with  a  change  in  energy  density  and  by  the  presence  of  energy 
density  induced  instability. 

Besides  convective  limits  on  repetition  rate  another  important  limit 
that  was  observed  was  the  occurrence  of  gas  poisoning  . 


CHAPTER  7 


COMPARISON  OF  SPARK  PREIONIZED  SYSTEM  AND  DIELECTRIC 
DISCHARGE  PREIONIZED  SYSTEM 

One  of  the  problems  encountered  in  the  experiments  with  the  high 
repetition  rate  discharge  was  that  of  arcing  from  the  cathode  to  the 
preionizing  spark  pins.  This  problem  caused  a  significant  decrease  in 
repetition  rate  capability  since  it  meant  that  not  only  the  electrodes 
but  the  spark  pins  as  well  had  to  be  cleared  of  discharge  products 
resulting  in  a  much  longer  clearing  time.  It  seemed  that  this  problem 
could  be  easily  eliminated  by  removing  the  spark  pin  and  going  to  another 
preionization  scheme.  The  other  popular  preionization  technique  is  to 
preionize  from  behind  a  screen  electrode  either  using  spark  sources  or  a 
dielectric  discharge.  The  dielectric  discharge  technique  was  chosen  for 
a  comparison  experiment  because  it  was  easy  to  construct  and  looked  like 
it  would  produce  very  uniform  preionization.  In  addition  Pace  and  Lacombe 
(22)  indicated  that  it  should  result  in  reduced  gas  poisoning  compared  to 
spark  preionization. 

7 . 1  CONSTRUCTION 

The  flowing  gas  discharge  system  was  modified  by  replacing  the  pin 
preionized  electrode  configuration  (Figure  4.7)  with  a  double  discharge 
configuration  (Figure  7.1).  In  the  double  discharge  configuration  one 
of  the  copper  Rogowski  electrodes  was  left  as  it  was  before.  The  other 
electrode  was  replaced  with  a  flat  copper  screen  electrode  that  had  a 
dielectric  coated  electrode  placed  behing  it.  This  resulted  in  a  slight 
increase  in  electrode  separation  from  24  mm  to  27  mm. 

The  dielectric  electrode  (Figure  7.2)  was  formed  by  encasing  a  copper 
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Figure  7.1  Dielectric  Discharge  Configuration.  The  Discharge 

is  Formed  Between  a  Copper  Rogowski  Profiled  Electrode 

and  a  Copper  Screen.  Preionization  is  provided  by  a  smaller 

discharge  between  a  dielectric  coated  electrode  and  the  screen. 
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Fig.  7.2  Dielectric  Electrode. 

A  copper  disc  1/8"  thick  and  2"  in  diameter 
is  encased  in  PVC.  The  PVC  on  the  front  surface 
is  1/8"  thick.  The  post  in  the  back  is  an 
electrical  connection  and  mates  with  a  banana  plug. 


197 


plate  in  PVC.  The  edges  of  the  copper  plate  were  rounded  to  reduce 
electric  field  concentration  and  the  plate  was  bonded  to  the  PVC  using 
epoxy.  The  PVC  on  the  front  of  the  plate  was  1/8"  thick  and  the  gap 
between  the  back  of  the  screen  and  the  dielectric  was  1  ram. 

Except  for  replacing  the  spark  pin  with  the  dielectric  electrode  the 
circuit  remains  unchanged  (Figure  7.3).  Note  that  the  voltage  to  the 
dielectric  discharge  is  reduced  by  using  the  peaking  capacitor  as  a 
capacitive  voltage  divider.  The  peaking  capacitor  corresponds  to  the  spark 
pin  coupling  capacitor  in  the  spark  pin  preionized  system  and  has  the 
same  value  of  capacitance.  Typical  discharge  voltages  and  current  waveforms 
are  shown  in  Figure  7.4.  They  are  similar  to  the  waveforms  for  the  pin 
preionized  system,  the  main  difference  being  that  there  is  a  smaller  initial 
current  spike. 

7.2  GAS  FLOW  CALIBRATION 

The  new  electrode  configuration  altered  the  flow  duct  spacing  so  it 
was  necessary  to  recalibrate  the  gas  flow  velocity  as  a  function  of  fan 
speed.  As  a  result  of  the  new  configuration  flow  velocity  increased 
slightly.  For  example,  in  60  torr  N^,  at  3600  rpm  the  flow  velocity  rose 
from  22  to  25  m/sec.  The  new  coefficients  for  the  calibration  curves 


are  given  in  Table  7.1. 


198 


H' 


i> 


i' 


•H 

a 

o 

n 

•H 

U 

a; 

60 

n 

a3 

a 

CD 

•H 

Q 

C 

O 

•H 

4-J 

03 

N 

•H 

c 

O 

•H 

03 

U 

Ph 

03 

60 

aj 

,£3 

O 

CD 

•H 

Q 

O 

•H 

U 

+J 

o 

03 

i— I 

03 

•H 

o 

co 


60 

•H 


199 


5  K  V/d  i  v 


10  0A  /d  i  v 


0.2  |is  ec/d  iv 

Figure  7.4  Typical  Voltage  and  Current  Waveforms  for  a  60  Torr  Discharge 
in  the  Dielectric  Discharge  Preionized  System.  The  discharge 
capacitor  was  10  nf  and  the  charging  voltage  was  10  KV. 
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gas 

60  torr  N» 

80  torr  CO^ 

700  torr  He 

400  torr  Mix 
699  torr  Air 


y  =  A  +  Bx  +  Cx^+  Dx^ 


A 

0.138677x10 

-0.492731x10 

-0.668171x10 

-0.761284x10 

-0.347010x10 


B 

0.275098x10 

0.549936x10 

0.386930x10 

0.510157x10 

0.597285x10 


C 

0.198219x10 

0.145971x10 

0.164515x10 

0.134211x10 

0.236838x10 


D 

-0.220352x10 

-0.179170x10 

-0.171502x10 

-0.127824x10 

-0.453158x10 


Table  7.1  Calibration  curve  coefficients  for  dielectric  discharge  preionized 
system . 


An  interesting  result  was  also  obtained  by  recording  a  velocity  profile 
between  the  cathode  and  anode  (screen) .  It  shows  an  asymetric  flow  with 
the  screen  having  a  much  thicker  boundary  layer  than  the  smooth  electrode 
(Figure  7.5).  This  is  not  surprising  because  the  screen  is  a  much  rougher 
surface  than  the  solid  electrode. 


7 . 3  RESULTS 

Two  parameters  were  measured  in  the  dielectric  discharge  preionized 
system  for  comparison  with  the  spark  pin  preionized  system.  They  were 
maximum  repetition  rate  as  a  function  of  energy  density  and  maximum 
repetition  rate  as  a  function  of  gas  flow  velocity.  Maximum  repetition 
rate  in  this  case  is  the  repetition  rate  at  which  a  discharge  instability 
appears  to  become  continuous,  which  is  how  it  was  defined  in  section  6.10. 

In  addition  discharge  intensity  profiles  were  monitored  with  the  photo¬ 
multiplier  for  the  purpose  of  establishing  representative  diameters  to 
use  for  calculating  energy  density,  and  also  to  monitor  discharge  uniformity. 
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Fig.  7.5  Velocity  Profile  Between  the  Electrodes. 
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As  before  the  discharge  diameter  is  taken  to  be  at  10%  of  the  peak  of 
the  intensity  profile.  Typical  intensity  profiles  are  shown  in  Figures 
7.6  to  7.8  and  the  discharge  diameter  was  typically  70  mm. 

Plots  of  maximum  repetition  rate  versus  energy  density  are  shown  in 
Figures  7.9  to  7.13  and  can  be  compared  with  the  same  results  for  the 
pin  preionized  system  summarized  in  Figures  6.52  to  6.56.  The  maximum 
repetition  rate  drops  off  with  increasing  energy  density  much  more  rapidly 
for  the  dielectric  discharge  preionized  system  than  does  the  repetition 
rate  for  the  pin  preionized  system.  There  are  two  likely  causes  for  this 
reduction  in  performance  of  the  screen  electrode  system.  One  is  that  the 
screen,  being  less  uniform  than  the  solid  electrode,  probably  causes  arcing 
to  set  in  at  a  lower  energy  density.  The  other  is  that  discharge  products 
probably  get  trapped  in  the  stagnant  gas  region  behind  the  screen  and 
contribute  to  discharge  instability.  An  observation  supporting  this 
hypothesis  is  the  fact  that  for  the  screen  electrode  system  repetition 
rate  induced  arcs  tended  to  take  6  or  7  seconds  to  build  up  in  60  torr 
and  a  similarly  long  time  to  die  out  when  the  repetition  rate  was 
reduced.  Although  this  long  time  for  arcs  to  build  up  and  die  out  was 
observed  in  only> it  does  suggest  that  trapping  of  discharge  products 
in  the  stagnant  gas  behind  the  screen  is  a  mechanism  affecting  arcing. 

In  the  dielectric  discharge  preionized  system,  as  was  the  case  in 
the  spark  pin  preionized  system,  maximum  repetition  rate  generally 
increases  linearly  with  flow  velocity  (Figures  7.14  to  7.17).  There  was 
however  one  notable  exception.  For  45  J/£  in  60  torr  N2  the  repetition 
rate  stayed  almost  constant  at  a  low  value  as  flow  velocity  increased 
(Figure  7.15).  This  also  suggests  that  the  repetition  rate  in  this  case 
is  being  limited  by  the  time  needed  to  clear  discharge  products  out  from 
behind  the  screen  rather  than  by  the  convection  of  gas  between  the 
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.6  Intensity  Profile  for  Screen  Electrode  System. 
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Fig.  7.7  Intensity  Profile  for  Screen  Electrode  System. 

Discharge  in  80  torr  CO2  sn  81.2  nf  capacitor. 
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Fig.  7.8  Intensity  Profile  for  Screen  Electrode  System. 

Discharge  in  700  torr  He  with  an  81.2  nf  capacitor. 
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Fig.  7.10  Maximum  Repetition  Rate  as  a  Function  of  Energy  Density  for  the 

Screen  Electrode  System.  Data  for  60  torr  N~  and  a  flow  velocity 
of  25  m/sec.  CR=1  is  the  repetition  rate  corresponding  to  a 
clearing  ratio  of  1. 
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Fig.  7.11  Maximum  Repetition  Rate  as  a  Function  of  Energy  Density  for 

the  Screen  Electrode  System.  Data  for  80  torr  CO^  and  a  Flow 
velocity  of  29  m/sec.  CR=1  is  the  repetition  rate  corresponding 
to  a  clearing  ratio  of  1. 
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Fig.  7.12  Maximum  Repetition  Rate  as  a  Function  of  Energy  Density  for  the 
Screen  Electrode  System.  Data  for  a  400  torr  1:1:8  laser  mix  and 
a  flow  velocity  of  28  m/sec.  CR=1  is  the  repetition  rate 
corresponding  to  a  clearing  ratio  of  1. 
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Fig.  7.13  Maximum  Repetition  Rate  as  a  Function  of  Energy  Density  for  the 
Screen  Electrode  System.  Data  for  a  500  torr  1:1:8  laser  mix 
and  a  flow  velocity  of  28  m/sec. 
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Fig.  7.14  Repetition  Rate  Versus  Gas  Flow  Velocity  for  Screen  Electrode 
System.  Data  for  80  torr  C02  and  48  3/1  energy  density. 
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Fig.  7.15  Repetition  Rate  Versus  Gas  Flow  Velocity  for  Screen  Electrode 
System.  Data  for  60  torr  N2  and  45  J/l  energy  density. 
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Fig.  7.16  Repetition  Rate  Versus  Gas  Flow  Velocity  for  Screen  Electrode 
System.  Data  for  700  torr  He  and  48  J/fc  energy  density. 
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Fig.  7.17  Repetition  Rate  Versus  Gas  Flow  Velocity  for  Screen  Electrode 
System.  Data  for  400  torr  of  1:1:8  CO^N^He  laser  mix  and 
22  J /  £  energy  density. 


215 


electrodes . 

For  the  dielectric  discharge  preionized  system  gas  poisoning  was 
also  observed  to  occur  in  the  laser  mixture.  Typically  poisoning  took 
about  four  minutes  in  the  400  torr  laser  mixture  using  a  40  nf  storage 
capacitor.  This  system  seems  to  provide  no  improvement  over  the  spark  pin 
preionized  system  with  regards  to  poisoning  time. 

7.4  CONCLUSIONS 

The  dielectric  discharge  preionized  system  did  not  operate  as  well 
as  the  spark  pin  preionized  discharge  system.  The  maximum  repetition 
rate  dropped  off  much  more  rapidly  as  energy  density  increased,  and  the 
maximum  energy  density  that  could  be  obtained  was  lower.  No  significant 
change  was  observed  in  the  time  that  the  system  took  to  poison  due  to 
CO^  dissociation.  The  spark  pin  preionized  system,  therefore,  appears 
to  be  the  superior  of  the  two.  This  is  somewhat  surprising  because  Pace 
and  Lacombe  (22)  claim  that  their  work  showed  the  dielectric  discharge 
preionized  system  to  be  better  than  spark  preionized  systems,  particularly 
with  respect  to  gas  poisoning. 

A  possible  explanation  for  the  poor  performance  of  the  dielectric 
discharge  preionized  system  is  that  discharge  products  can  get  trapped 
in  the  stagnant  gas  behind  the  screen  electrode.  This  could  be  tested  by 
arranging  for  gas  flow  behind  the  screen  and  rerunning  the  experiment  to 
see  if  it  results  in  an  improvement  in  maximum  repetition  rate. 
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CHAPTER  8 


FLOWING  GAS  DISCHARGE  SYSTEM  STREAK  CAMERA  OBSERVATIONS 

INTRODUCTION 

This  chapter  describes  some  observations  that  were  made  using  a 
streak  camera  to  study  discharges  in  60  torr  H  .  The  experiment  had 
two  main  goals.  The  first  was  to  measure  the  gas  expansion  that  had 
been  predicted.  Since  N^  had  a  luminous  afterglow  it  was  hoped  that  the 
expansion  of  the  afterglow  would  be  visible.  The  other  goal  was  to 
produce  some  streak  photographs  which  showed  typical  examples  of  the 
discharge  instabilities  that  occurred  in  this  system. 

8.1  APPARATUS 

These  observations  were  made  using  a  TRW  model  ID  image  converter 
camera  with  a  model  5B  wide  range  plug-in.  Writing  rates  for  the  plug-in 
are  100,  50,  25,  10  and  5  mm/ysec.  A  model  46A  trigger  delay  generator 
was  used  to  trigger  the  camera.  These  units  were  borrowed  from  the 
Plasma  Laboratory. 

A  block  diagram  of  the  circuit  arrangement  is  shown  in  Figure  8.1. 

A  master  trigger  generator  synchronized  the  experiment  by  triggering 
the  discharge  and  the  scope.  The  scope  triggered  the  trigger  delay 
generator  which  in  turn  triggered  the  camera.  Both  the  scope  and  the 
trigger  delay  generator  could  be  used  for  adjusting  the  delay  before  the 
camera  was  triggered  but  only  the  trigger  delay  generator  was  used  for 
this  function.  Electrical  noise  was  a  problem  with  this  experiment  and 
initially  caused  the  camera  to  trigger  spuriously,  and  to  terminate 
streaks  prematurely.  Both  radiated  noise  and  noise  conducted  along  the 
cables  contributed  to  the  problem.  By  keeping  the  cables  short  and 
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Fig.  8.1  Block  Diagram  of  Streak  Camera  Experiment  Showing  Triggering 
Chain  of  Command  and  Signal  Paths.  The  current  probe  and 
camera  monitor  paths  are  signal  paths.  All  others  are 


trigger  paths. 
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shielding  them  with  copper  pipe  the  effect  of  radiated  noise  was  reduced 
to  the  point  where  reasonably  reliable  results  could  be  obtained.  More 
work  was  needed  however.  Better  reliability  could  be  achieved  by  shielding 
the  camera  with  a  copper  mesh  cage  or  a  metal  box,  and  some  of  the  cables 
should  have  been  isolated  using  op to-isolators .  However  since  the  camera 
was  available  on  a  short  term  loan  there  was  not  time  to  take  these 
additional  steps. 

Measurements  were  made  in  60  torr  N^.  Charging  voltage  for  the 
discharge  was  typically  10  KV  and  the  storage  capacitors  were  40  nf  and 
80  nf .  Because  the  spark  pin  preionized  discharge  system  proved  superior 
to  the  dielectric  discharge  preionized  system  it  was  used  for  these 
experiments.  Figure  8.2  shows  the  orientation  of  the  streak  camera  slit 
with  respect  to  the  discharge.  This  slit  orientation  allows  monitoring 
the  discharge  diameter  as  a  function  of  time. 

8 . 2  RESULTS 

The  attempt  to  record  gas  expansion  by  monitoring  the  afterglow 
diameter  was  unsuccessful.  The  main  reason  was  that  the  camera  had 
insufficient  sensitivity  to  record  the  afterglow  for  longer  than  about 
20  ysec  after  the  discharge.  Also  the  image  was  reduced  from  the  original 
discharge  by  a  factor  of  four  which  decreased  the  sensitivity  of  the 
measurement . 

An  estimate  of  expected  expansion  was  calculated  for  an  energy 
density  of  60  J/£,  the  energy  density  produced  by  an  80  nf  capacitor 
charged  to  10  KV.  The  afterglow  diameter  should  expand  from  60  mm  to 
96  mm,  or  a  total  expansion  of  36  mm.  Using  shock  theory  the  particle 
flow  velocity  was  estimated  to  be  .34  mm/ysec,  so  the  diameter  would  expand 
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Fig.  8.2  Orientation  of  Streak  Camera  Slit  with  Respect  to  the  Discharg 
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at  twice  that  rate  or  .68  mm/ysec.  About  53  ysec  would  be  required 
for  the  entire  expansion.  The  longest  streak  available  was  5  ysec  streak, 
since  the  expansion  is  reduced  by  a  factor  of  four  on  the  image,  the 
image  diameter  would  change  by  only  .86  mm  which  is  too  small  to  judge 
accurately  due  to  the  fuzzy  edges  of  the  image.  By  using  delay  the 
afterglow  could  be  examined  for  up  to  20  ysec  after  the  discharge.  Over 
20  ysec  the  estimated  change  in  image  diameter  is  roughly  3  mm  which 
should  be  visible.  However  the  results  are  ambiguous  because  a  variation 
in  light  intensity  due  to  a  change  in  exposure  can  result  in  a  diameter 
variation  of  this  much.  For  example,  the  image  diameter  at  20  ysec  is 

14  mm  and  required  multiple  exposures  to  be  visible.  Multiple  exposures 
may  result  in  a  widening  of  the  streak  due  to  image  jitter  from  shot 

to  shot.  However  the  initial  discharge  image  diameter  was  between  11  and 

15  mm  depending  on  which  exposure  is  used.  If  it  was  11  mm  then  roughly 
the  correct  amount  of  expansion  seems  to  have  occurred  in  20  ysec,  but 
if  it  was  15  mm  then  no  expansion  occurred.  As  a  result  no  definite 
conclusion  about  gas  expansion  can  be  drawn  from  these  pictures. 

The  streak  camera  measurements  produced  pictures  of  typical  discharge 
instabilities.  Figure  8.3  shows  the  discharge  current  waveform  for  a 
uniform  discharge,  with  the  streak  camera  monitor  pulse,  and  the 
synchronization  between  the  two.  A  streak  of  the  uniform  discharge  is 
shown  in  Figure  8.4a.  The  afterglow  appears,  faintly,  following  the 
discharge.  An  energy  induced  glow  to  arc  transition  is  shown  in 
Figure  8.4b  and  a  repetition  rate  induced  glow  to  arc  transition  in 
Figure  8.4c.  It  can  be  seen  that  the  glow  to  arc  transitions  are  character¬ 
ized  by  a  collapse  in  the  discharge  diameter.  Also  that  the  two  types 
of  collapse  occur  at  different  positions,  with  the  repetition  rate  induced 
glow  to  arc  occurring  on  the  downstream  edge  of  the  discharge.  The 
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Fig.  8.3  Synchronization  Between  Discharge  Current  Pulse  (Upper  Trace)  and 
Streak  Camera  Monitor  Pulse  (Lower  Trace) .  The  streak  duration 
is  2.6  usee. 
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Fig.  8. A  Streak  Camera  Pictures  of  (a)  A  Uniform  Discharge,  (b)  An  Energy 
Induced  Glow  to  Arc  Transition  and  (c)  A  Repetition  Rate  Induced 


Glow  to  Arc  Transition.  The  streak  duration  is  2.6  ysec. 
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results  in  Figure  8.4  are  for  a  40  nf  storage  capacitor  charged  to  10  KV. 

In  Figure  8.5  a  higher  charging  voltage,  18  KV,  is  used  to  produce  some 
examples  of  severe  energy  induced  glow  to  arc  transitions.  The  discharge 
current  lasts  about  1  usee  and  the  streak  lasts  approximately  5  usee., 
so  the  arc  channel  stays  luminous  well  into  the  afterglow  region.  If 
the  repetition  rate  is  very  high  an  arc  may  result  instead  of  a  glow  to 
arc.  Figure  8.6a  shows  a  uniform  discharge  and  Figure  8.6b  shows  a 
repetition  rate  induced  arc. 

Similar  results  were  obtained  for  discharges  in  60  torr  using  an 
80  nf  storage  capacitor. 

8.3  CONCLUSIONS 

The  streak  camera  measurements  failed  to  show  whether  or  not  gas 
expansion  occurred.  With  respect  to  discharge  arcing  they  showed  the 
kinds  of  arcing  that  could  be  encountered.  The  energy  induced  instability 
is  a  glow  to  arc  transition.  Repetition  rate  induced  instabilities  can 
occur  either  as  glow  to  arc  transitions  or  as  arcs.  The  two  types  of 
repetition  rate  instability  would  be  indistinguishable  visually  without 
the  streak  camera  because  the  eye  would  superimpose  the  arc  on  following 
glow  pulses  so  both  effects  would  visually  appear  as  a  glow  to  arc. 

These  pictures  do  show  that  repetition  rate  induced  arcs  and  energy  induced 
arcs  can  be  readily  distinguished  from  their  positions  in  the  discharge. 
Energy  induced  arcing  occurs  on  the  upstream  edge  or,  more  often,  toward 
the  center  of  the  electrodes.  Repetition  rate  induced  arcing  always 
occurs  on  the  downstream  edge  of  the  electrodes.  The  fact  that  repetition 
rate  induced  arcing  occurs  on  the  downstream  side  of  the  discharge  supports 
the  conclusion  that  it  is  a  convective  effect  and  is  due  to  gas  convection 
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Fig.  8.5  Examples  of  Energy  Induced  Glow  to  Arc  Transitions.  The 
streak  duration  is  4.7  ysec. 
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Fig.  8.6  (a)  Shows  a  Uniform  Glow  and  (b)  Shows  a  Repetition  Rate  Induced 


Arc.  The  streak  duration  is  4.7  ysec. 
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not  having  enough  time  to  clear  discharge  products  from  the  discharge 
region  between  pulses. 


CHAPTER  9 


CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  STUDY 

9.1  LASER  SYSTEM  RESULTS 

These  experiments  showed  that  a  high  repetition  rate  closed  cycle 
CO^  TEA  laser  is  feasible  to  construct  and  operate.  This  laser  was 
pumped  by  a  10  stage  Marx  generator.  Gas  flow  through  the  discharge 
region  was  approximately  100  m/sec,  giving  an  expected  repetition  rate 
of  about  2  KHz.  A  typical  discharge  pulse  length  was  1  psec.  Although 
pulse  stretching  was  attempted  it  was  not  very  successful,  with  maximum 
practical  pulse  lengths  being  in  the  6  to  7  ysec  range.  Maximum 
repetition  rates  for  reliable  operation  were  in  the  200  to  300  Hz  range. 
The  best  optical  efficiency  was  6%  with  16  J  per  pulse  energy  input  and 
1  J  per  pulse  optical  output.  This  was  achieved  by  using  an  unstable 
resonator  with  a  square  skimmer.  Other  resonators  were  tested  but  found 
to  be  less  effecient.  The  maximum  average  power  output  was  about  200  W 
(1  J/pulse  at  200  pulses  per  second) .  The  optical  system  design  still 
needed  improvement  as  vibration  caused  a  large  pulse  to  pulse  variation 
in  the  optical  energy  output. 

Two  major  problems  existed  with  this  laser  however.  The  first  was 
the  reduction  in  achieved  repetition  rate  from  the  original  estimate  by 
an  order  of  magnitude  -  from  2  KHz  to  200  Hz.  The  second  was  the 
suspected  poisoning  of  the  laser  system  after  a  short  time  due  to  dis¬ 
sociation  of  C02  into  CO  and  0^.  Tbe  discharge  experiments  which 
followed  were  to  identify  the  causes  of  discharge  arcing  and  to  confirm 
the  presence  of  gas  poisoning. 

9.2  SHOCK  WAVE  MEASUREMENTS 

One  explanation  considered  for  the  reduction  in  repetition  rate  was 
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that  possibly  discharge  generated  shock  waves  were  producing  a  nonuniform 
gas  density  and  so  causing  discharge  arcing.  By  monitoring  the  deflection 
of  a  probe  laser  beam  that  traversed  the  discharge  region  shock  waves 
could  be  studied.  They  were  observed  both  in  static  and  flowing  gas. 

In  static  gas  a  large  thermal  gradient  was  observed  after  the 
discharge  but  it  was  absent  in  flowing  gas.  This  is  because  the  hot 
gases  had  to  dissipate  by  diffusion  in  the  static  gas  whereas  in  the 
flowing  gas  they  were  swept  away  very  quickly  by  convection. 

In  both  static  and  flowing  gas  shock  waves  were  observed  that 
originated  from  gas  heating  in  the  cathode  fall  region  of  the  discharge. 
These  shock  waves  were  followed  by  an  oscillation  which  seemed  to 
originate  from  the  Marx  spark  pins.  The  shock  waves  from  the  cathode 
and  from  the  spark  pins  died  out  in  approximately  0.5  msec  so  they  would 
not  contribute  to  discharge  arcing  for  repetition  rates  less  than  2  KHz. 

As  a  result  it  can  be  concluded  that  shock  waves  were  not  a  source  of 
discharge  arcing  in  this  laser  system. 

9.3  FLOWING  GAS  DISCHARGE  SYSTEM  MEASUREMENTS 

The  tests  in  the  flowing  gas  discharge  system  showed  that  several 
factors  affected  the  repetition  rate  induced  arcing  limit  in  the  high 
repetition  rate  laser.  The  repetition  rate  was  linearly  proportional 
to  flow  velocity  which  indicates  that  repetition  rate  induced  arcing  is 
a  convective  effect.  The  types  of  instabilities  which  occurred,  and 
which  are  shown  in  Fig.  6.49,  also  suggest  that  it  is  a  convective 
effect  because  the  repetition  rate  induced  instabilities  always  occur 
on  the  downstream  side  of  the  discharge. 

The  repetition  rate  at  which  repetition  rate  induced  arcing  occurred 
was  independent  of  pulse  length  but  decreased  with  increasing  energy 
density.  Dzakowic  and  Wutzke  have  proposed  that  this  decrease  is  due  to 
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expansion  of  the  discharge  products  so  that  more  time  is  required  to 
convect  them  out  of  the  discharge  region.  However  repetition  rate 
decreased  more  rapidly  with  increasing  energy  density  than  predicted  by 
their  theory.  The  afterglow  measurements  supported  the  concept  of 
discharge  product  expansion  because  the  afterglow  plug  diameter  increased 
with  increasing  energy  density.  The  data  fit  the  expansion  theory  well 
for  but  not  as  well  for  CO^  (Figs.  6.34  and  6.35).  However  there 
was  some  ambiguity  because  the  afterglow  plug  diameter  seemed  to  be 
smaller  than  the  discharge  diameter  as  defined  by  the  discharge  intensity 
profile  measurements.  The  afterglow  measurements  also  supported  the 
conclusion  that  the  repetition  rate  is  essentially  limited  by  the  time 
needed  for  convective  clearing,  since  fractional  frequencies  based  on 
the  time  needed  to  clear  the  afterglow  products  past  the  spark  pins  in 
agreed  well  with  observed  fractional  frequencies  (Fig.  6.38). 

The  streak  camera  measurements  were  intended  to  observe  the 
afterglow  expansion  predicted  by  Dzakowic  and  Wutzke  but  were  inconclusive 
on  this  point.  They  did  however  give  useable  pictures  of  the  various 
arcing  modes. 

The  experiments  in  the  flowing  gas  discharge  system  showed  that 
geometry  is  an  important  factor  in  limiting  the  maximum  repetition  rate. 
Upstream  spark  pins  were  shown  to  be  detrimental  to  system  performance. 

In  the  case  of  a  downstream  spark  pin  the  repetition  rate  at  high  energy 
densities  is  limited  by  having  to  clear  the  discharge  products  not  only 
from  the  electrodes  but  past  the  spark  pins  as  well.  Otherwise  arcing 
would  occur  from  the  electrodes  to  the  spark  pins.  Removing  the 
downstream  spark  pins  and  preionizing  from  behind  a  screen  electrode 
does  not  necessarily  solve  the  problem  either.  The  solid  electrodes 
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with  the  downstream  preionizing  spark  pins  performed  better  than  the 
system  with  preionization  from  a  dielectric  discharge  behind  a  screen 
electrode.  This  was  probably  because  discharge  products  got  trapped 
in  the  stagnant  gas  behind  the  screen. 

Geometry  was  the  most  important  factor  in  reducing  repetition 
rate  as  energy  density  increased  in  the  flowing  gas  discharge  system. 
Consider  first  the  case  of  the  spark  pin  preionized  system.  As  the 
discharge  energy  density  increased  so  did  the  density  of  discharge 
products  left  over  after  the  pulse.  At  low  energy  densities  it  was 
sufficient  to  clear  them  out  from  between  the  electrodes  in  order  to 
produce  a  stable  discharge.  But  as  the  density  of  discharge  products 
increased  a  point  was  reached  where  they  could  form  a  conductive  bridge 
from  the  electrodes  to  the  spark  pins.  When  this  point  was  reached  it 
became  necessary  to  clear  the  discharge  products  past  the  spark  pins 
in  order  to  prevent  arcing  to  the  spark  pins  from  occurring.  This 
results  in  a  greater  decrease  in  repetition  rate  than  expansion  of  the 
discharge  products  does.  In  the  case  of  the  screen  electrode  with  the 
dielectric  discharge  for  preionization  discharge  products  trapped  in 
the  stagnant  gas  behind  the  screen  appeared  to  be  the  main  factor 
limiting  repetition  rate,  so  once  again  the  system  geometry  was  important 

Gas  poisoning  was  also  observed  in  the  flowing  gas  discharge 
system  and  is  believed  to  be  due  to  the  dissociation  of  CO^  into  CO  and 
0^  in  the  discharge.  As  a  result  of  poisoning  the  sealed  off  lifetime 
of  the  system  at  100  Hz  repetition  rate  is  roughly  4  minutes. 


9.4  CONSIDERATIONS  FOR  HIGH  REPETITION  RATE  C02  TEA  LASER  DESIGN 

Several  important  considerations  have  arisen  from  this  study  which 
would  be  useful  in  future  high  repetition  rate  C02  TEA  laser  designs. 
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A  Marx  generator  switched  system  in  which  the  Marx  switches  provide 
preionization  appears  to  be  feasible.  However  the  Marx  pins  must  be 
downstream  from  the  discharge  in  order  to  avoid  adversely  affecting  the 
discharge.  In  addition  the  repetition  rate  will  be  limited  by  the  time 
needed  to  clear  discharge  products  from  between  the  electrodes  and  past 
the  spark  pins  rather  than  just  from  between  the  electrodes.  Using 
a  peaking  capacitor  in  connection  with  the  Marx  generator  improves  the 
discharge  stability. 

Gas  poisoning  was  shown  to  be  a  problem  in  a  sealed  system  and 
dissociation  of  CO^  into  CO  and  0^  is  the  suspected  cause.  Therefore 
future  high  repetition  rate  systems  would  either  have  to  be  run  with  a 
make-up  gas  flow  or  with  some  catalyst  system  to  recombine  the  CO  and 
0^*  Work  is  being  done  by  other  researchers  on  such  catalyst  systems 
(ref .  23) . 

The  conversion  efficiency  of  electrical  energy  to  optical  output 
was  low  for  a  CO^  laser,  being  about  6%.  However  the  gain  was  low 
and  the  gain  path  length  of  60  cm  was  fairly  short  for  a  low  gain  system. 
In  future  lasers  of  this  type  the  optical  efficiency  could  probably  be 
improved  by  using  a  longer  discharge. 

9.5  SUGGESTIONS  FOR  FURTHER  WORK 

Several  possibilities  for  future  work  arise  out  of  this  project. 

One  experiment  would  be  to  arrange  for  gas  flow  behind  the  screen  in 
the  dielectric  discharge  preionized  system  then  compare  its  performance 
with  the  solid  electrode  system.  It  is  possible  that  with  this  change 
the  dielectric  discharge  preionized  system  may  prove  to  be  superior. 

Work  can  also  be  done  on  a  means  of  controlling  gas  poisoning  using  a 
catalyst.  Extension  of  high  repetition  rate  TEA  laser  techniques  to 
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other  gases  is  another  possible  project.  Olson  et  al.  (21)  have  done 
this  for  He-Xe  and  several  other  rare  gas  mixtures. 

The  photomultiplier  monitoring  could  be  expanded  as  well.  The 
intensity  profile  monitoring  could  be  expanded  into  a  tool  for  measuring 
discharge  uniformity.  Also,  the  photomultiplier  could  be  shuttered  to 
protect  it  from  the  initial  discharge  flash  when  doing  afterglow 
monitoring  in  and  CO^.  This  might  permit  observing  the  afterglow 
between  the  electrodes  and  resolve  why  the  afterglow  plug  diameter 
appears  to  be  smaller  than  the  discharge  diameter. 

The  streak  camera  results  could  also  use  more  work.  The  camera 
triggering  reliability  at  high  repetition  rates  should  be  improved.  Also 
the  camera  could  be  used  to  take  high  speed  framing  photographs  which 
would  better  illustrate  the  spatial  development  of  the  discharge  and 
arcs  than  the  streak  photographs. 
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Appendix 

Listing  of  Computer  Routines  used  in  these  experiments. 
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The  APL  functions  were  used  for  the  circuit  model  of  section  3.3. 
A  sample  output  is  included  with  the  listings.  The  function  CUOULI  was 
written  by  Hugh  Woods  and  the  function  PLOT  is  from  the  University's  APL 
Computing  Library. 
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NB 

1  1  60 


VV 

1E"8  1.75E'9  4E'6  20  40000  60000 
NB  DRIVER  VV 
VB  =40000 

TB  =  1.043841304E'7 
VC1TB  =  53000 
ILTB  =  1130. 265456 
7 . 1428571 4 3E20  1 .678571429E14 
22527764.73  3021831.919 

0  '4751361.663 

1776. 61 2791e*< -22527764. 73T) 


AN2e*(-GAxT)  +  AN3<€*< -AL)  )C0S(BExT  -  F'SI) 


28571428.57  1 
3021831.919 
4751361.663 

+  3343 . 061996c  * ( -  302183 1 . 9 1 9T ) COS < 4751 361 . 663x T 


O.OOOOOOOOOEO 
2.609603261E~8 
5 . 21 9206522E "8 
7 • 828809783E ' 8 
1 ,043841304E"7 
1 .043841304E'7 
1 .264240766E'7 
1.4846 40227E"7 
1 ,705039688E~7 
1.925439150E'7 
2, 14583861  IE *7 


O.OOOOOOOOOEO 
2 • 890887418E3 
123622483E4 
409109952E4 
OOOOOOOOOE4 
OOOOOOOOOE4 
3. 230 7 69 78 IE 4 
2.912862630E4 
2.837280698E4 
2 . 879 156434E4 
2. 9649690 71E4 


2 • 366238072E'7 
2.586637534E'7 
2.807036995E"7 
3 . 027436 4 56E~ 7 
3.247835918E'7 
3.468235379E'7 
3.688634840E"7 
3 . 909034302E "7 
4. 129433763E'7 


3 • 052537640E4 
3. 118795814E4 
3. 152308559E4 
3. 148680242E4 
3. 107727825E4 
3 . 031 7334 1 4E4 
2.924358518E4 
2. 789965584E4 
2 . 633191 604E4 


4. 349833 224E'7  2.458679014E4 
4.570232686E"7  2 . 270905856E4 
4 . 7906321 47E' 7  2 . 074079523E4 


5.011031608E'7 
5.231431070E"7 
5.451830531E*7 
5 . 67222 999 3E  ~ 7 
5. 89 2629 45 4E' 7 
6 . 1 13028915E'7 
6.333428377E~7 
6 . 553827838E'7 
6.774227299E'7 
6.994626761E'7 
7.215026222E'7 
7.435425683E'7 
7.655825145E"7 
7 . 876224606E*  7 
8 . 096624067E  ~  7 
8.317023529E'7 
8. 53 7 4 229 90E" 7 
8 . 75782245 1E"7 
8.978221913E‘7 
9. 198621374E“7 
9 • 419020835E"7 
9 . 639420297E*7 
9. 859819758E'7 
1 . 008021922E"6 
1 ,030061868E'6 
1. 05210181 4E "6 
1 . 074141 760E" 6 
1 .096181706E'6 


1 .872072020E4 
1 .668382957E4 
1 .466121558E4 
1 .268002017E4 
1 .076348477E4 
8.931070325E3 
7.1 986292 72E 3 
5.578615506E3 
4 . 080321 732E3 
2 • 710135525E3 
1 .471806999E3 
3.667221093E2 
"6.058234444E2 
'1 .448353995E3 
“2. 164950937E3 
‘ 2.760999982E3 
' 3 . 242951068E3 
'3 . 6 18092827E3 
'3.894343075E3 
“4.080056343E3 
'4. 183849111E3 
“4. 214443 062E3 
'4. 180526361E3 
'4 . 0906327 13E 3 
”3. 953037 730E3 
'3. 775671937E3 
”3. 56604959 4E3 
~3 . 331212395E3 


1.503925361) 
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1 « 118221453E~A 
1 ♦ 1402A1599E“A 
1.142301545E*4 
1 • 18434 14? IE *4 
1.204381437E'A 
1.228421383E*4 
1 • 250461 3 29E "6 
1 .272501274E*4 
1 • 294541222E* A 
1 .31A5811ABE*4 
1 . 338421 1 1 4E * A 
1.340441040E*A 
1 .382701004E*4 
1 .404740952E*4 
1 ,424780898E*4 


*3.077A8A988E3 
*2.81 1455213E3 
*2.537935897E3 
“2 . 261977030E3 
*1 .987857133E3 
" 1 • 719294A40E3 
“1.459464296E3 
*1.21 101904AE3 
*9 . 761 170907E2 
"7.564524119E2 
*5. 53288290 1E2 
*3 . A74928233E2 
*1 .995757092E2 
“4.97254010AE1 
8.215257077E1 
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The  program  REGRESS  was  used  to  generate  the  gas  flow  calibration 
curve  coefficients  shown  in  Tables  5.1  and  7.1. 
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DOUBLE  PRECISION  P(200),T(20) 

INTEGER  ID (4) /*  *,»  V 

DIMENSION  X  (100)  ,  I  ( 100)  ,X  I  (  100)  ,  Y  1  (  100) 

DIMENSION  C  (8)  ,  S  (8)  ,  A(5),  fl(5) 

NAMELIST  /DATA/  ID,  X,  Y,  N 
DO  5  1=  I,  100 
X(I)  =0. 

5  Y  (I)  =0. 

WRITE  (6,  1  00) 

100  FORMAT  ( '  0  ENT  ER  DATA') 

READ (5, DATA) 

WRITE  (6  ,  105)  ID 

105  FORMAT  ('0CALIBR ATION  CURVE  FITTING  FOR'/'  ',4A4) 

WR ITE  (6,  1  10)  X ,  Y,  N 

NO  FORM  AT  ( '  RPM='/20(»  ',SFI0.  I/)/*  VEL  =  '/20('  ',5F10.1/)/'  N  =  »,I4) 
RSQ= 1 0  0 

MD=  1 

15  CONTINUE 
DO  6  1= 1,8 
S  (I)  =0. 

6  C(I)=0. 

DO  7  1=1,5 
B  (I) =0. 

7  A  (I)  =0  . 

DO  20  1=1,  100 
X  1  (I)  =X  (I) 

20  Y  1  (I)  =  Y  (I) 

CALL  RLFOTH  ( X  1  ,  Y  1  ,  N  ,  RSQ  ,M  D  ,  ID  1  ,P  , C,  S  ,  A,  F  ,  I ES) 

C  WRITE (6,  120)  X  I  , Y I 

120  FORMAT  ( *  0  X  1= '  /20  ( •  »,5FI0.6/)/'  YI=*/20(»  ',5FI0.6/)) 

WRITE  (6 ,  12  5)  ID  1 
125  FORMAT ( '0ID 1= '  ,13) 

C  WRITE  (6  ,  130)  C 

130  FORMAT ( ' 0C=' ,2 (/'  » , 5G  1  5. 6) ) 

WRITE  (6,  135)  5 

135  FORMAT  (  '0S=  ' ,  2  (/'  »,5G15.6)) 

A  NUM=N -ID  I  -  I 
V AR  =  S  (ID  1  +  3)  /ANDM 
WRITE (6, 137)  VAR 
137  FORMAT  ('  VAR  IA  NCE  =  '  ,  G  1  5 . 6) 

C  WRITE  (6,  140) A 

C  140  FORMAT  ('0A=' ,2  (/»  »,5EI5.6)) 

C  WRITE  (6,145)  B 

C  I  45  FORMAT  ( ’  0  B=  '  ,  2  (/'  »  ,  5E  !  5.  6)  ) 

WRITE  (6,  150)  IER 
150  FCRMAT('0IER=',I3) 

CALL  FLOP  DC  (X  I  ,N,  A,B,C,  ID  1,  1,?,YI,IER) 

WRITE  (6,  160)  Y  1 

160  FORMAT  (’OY  l='/20('  »  ,  5  F  1 0 . 6 /)  ) 

CALL  PL  DO  PM  (C  ,ID1  ,  A  ,  B  ,  T) 

WRITE  (6  ,  13  0)  C 
MD=  MD ♦  I 

IF  (MD.  LE.  4)  GOTO  15 
STOP 
END 

SAMPLE  DATA  DFCK  FOLLOWS 

4DATA  ID=  *  400  ','TORR',’  HEL','IrJM  ’, 

X=0. ,750. ,  1034. ,  1463. ,2000. ,24 00., 3 000. ,3571.  , 

Y=0. ,2  .4,4.  I  ,6.5,  10.  I,  13. 2,  17. 6,22.2, 

N=  8  A  END 


FILE 


* . 
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The  program  ANALYSE  was  used  to  do  the  fractional  frequency 
analysis  of  repetition  rate  data.  In  addition  to  the  sample  output  shown 
in  the  appendix  this  program  also  generates  a  repetition  rate  versus  flow 
velocity  plot.  The  plot  corresponding  to  the  sample  output  is  Figure  6.16. 
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list  3nalyse2 


> 

1 

C 

PROGRAM  WRITTEN  OCT.  30,  1978 

x. 

2 

C 

LATEST  REVISION  NOVEMBER  10,  1978 

V, 

3 

C 

> 

4 

DIMENSION  CAL ( 4 ) 

> 

5 

DIMENSION  RPM<24),  REPRAT<24),  PREP<24) 

> 

6 

INTEGER  ID  <  4 ) / '  ','  ' , '  '/ 

> 

7 

DIMENSION  XL INE ( 26 ) ,  YLINE<26> 

> 

8 

DIMENSION  FAN (24) , DISCH ( 24 ) , VEL ( 24 ) , FPREP ( 24 ) 

V 

,/ 

9 

NAMELI ST/DATA/ID > CAP > VOLT , DI A > CAL » FAN r DI SCH f N 

X 

,/ 

10 

DO  5  1=1,24 

X. 

/ 

11 

VEL ( I ) =0 ♦ 

x 

12 

FAN ( I ) =0 . 

> 

13 

DISCH  < I ) =0 . 

V 

14 

R'PM  (  I )  =0  . 

> 

15 

REPRAT ( I ) =0 . 

> 

16 

PREP ( I ) =0 . 

'? 

17 

FPREP ( I ) =0 . 

18 

XLINE ( I ) =0 . 

> 

19 

YL I NE ( I ) =0 ♦ 

y' 

20 

5 

CONTINUE 

> 

21 

XLINE ( 25 ) =0 . 

> 

22 

YL INE ( 25 ) =0 . 

> 

23 

XLINE ( 26 ) =0 . 

V 

./ 

24 

YL INE ( 26  >  =0 . 

> 

25 

C 

> 

26 

C 

INPUT  DATA 

V 

,/ 

27 

C 

> 

28 

UR  I TE ( 6 , 90 ) 

\ 

,/ 

29 

90 

FORMAT( 'OENTER  DATA') 

V 

/ 

30 

READ <5, DATA) 

> 

31 

ENERGY=CAP# VOLT* VOLT/2 . 

> 

32 

C 

> 

33 

C 

ELECTRODE  SEPARATION  IS  15/16  INCH=23 . 8125MM 

34 

C 

> 

35 

RAD=DIA/2 . 0 

> 

36 

V0L=3. 141593* RAD*RAD *23.8125*1. E-6 

/ 

37 

EDNSTY=ENERGY/VOL 

> 

38 

C 

> 

39 

C 

CALCULATE  RPM  AND  REPETITION  RATE 

■> 

40 

C 

> 

41 

DO  10  1=1 , N 

42 

RPM( I ) =60, * ( 1 . / <  FAN( I ) * . 001 ) ) 

43 

10 

REPRAT  <I)  =  1./(DISCH(I)#,001 ) 

> 

44 

c 

> 

45 

c 

CALCULATE  FLOW  VELOCITIES  CORRESPONDING  TO  OBSERVED 

46 

c 

X. 

s 

47 

DO  15  1=1, N 

> 

48 

T1=CAL<2)*RPM(I) 

/• 

49 

T2=CAL<3)*RPM(I)*RPM( I) 

.> 

50 

T3=CAL ( 4 ) *RPM( I ) *RPM< I ) *RPM< I ) 

X. 

✓ 

51 

VEL( I )=CAL( 1 )+Tl+T2+T3 

> 

52 

15 

CONTINUE 

> 

53 

c 

> 

54 

c 

CALCULATE  IDEAL  REP  RATE  CORRESPONDING  TO  OBSERVED 

'■> 

55 

c 

REP  RATES 

> 

56 

c 

;> 

57 

DO  25  I = 1 , N 

> 

58 

25 

PREP (I )=1000.#VEL< I)/DIA 

> 

59 

c 

s' 

60 

c 

OUTPUT  DATA  AND  CALCULATED  RESULTS 

/* 

61 

c 

> 

62 

WRITE (6, 130)  ID 

> 

63 

130 

FORMAT( '-REPETITION  RATE  ANALYSIS'/'  ',4A4) 

> 

64 

WRITE (6, 132)  CAP,VOLT,DIA,EDNSTY,N 

' 


. 


65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 
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132  FORMATS  CAPAC I TANCE* ' > E8 . 3 » '  FARADS'/'  VOLTAGE- 'f E8 . 3 » '  VOLTS'/ 
1  '  DISCHARGE  DIAMETER-' fF8.2f 'MM'/'  ENERGY  DENSITY- ' r F8 . 2 > 

1  'J/L'/'  N= ' f I 4 ) 

WRITE  <  6> 134 )  CAL 

134  FORMAT ( ' OCALIBRATION  CURVE  COEFFICIENTS'/'  ' >4015.6) 

WRITE(6> 136) 

136  FORMAT < ' OFAN  PERIOD  AND  INTERPULSE  PERIOD'/'  ' > T7 > ' MSEC ' > T17 > 

1  'MSEC') 

DO  30  1=1 >N 

WRITE<  6.138)  FAN < I ) >  DISCH< I ) 

138  FORMAT < '  ' >2F10.3) 

30  CONTINUE 

WRITE( 6> 140) 

140  F0RMAT( 'OCALCULATED  VALUES'/'  '>T8>'FAN  SPEED '> T24 >' FLOW  VEL'f 
1  T39  > ' REP  RATE' >T48> ' IDEAL  REP  RATE'/ 

1  '  ' >T14> 'RPM' >T27> 'M/SEC' >T45> 'HZ' >T60> 'HZ' ) 

DO  40  1  =  1 >  N 

WR I TE ( 6 > 1 42 )  RPM ( I ) >  VEL(I)>  REPRAT(I)>  PREP(I) 

142  FORMAT < '  ' >4F15. 1) 

40  CONTINUE 

LEAST  SQUARES  FIT  OF  REP  RATE  DATA  TO  COMPARE  WITH 
IDEAL  REP  RATE.  LEAST  SQUARES  FIT  IS  FOR  REP  RATE 
VERSUS  GAS  FLOW  VELOCITY 
PLOTTED  FIT  ASSUMES  A  ZERO  INTERCEPT. 

x=o  • 

Y=0  ♦ 

XY=0  • 

XX=0. 

DO  45  1  =  1 >  N 
XX=XX+VEL< I )*VEL< I ) 

XY=XY+VEL< I )*REPRAT< I ) 

X-X+VEL ( I ) 

45  Y-Y+REPRAT ( I ) 

A  IS  SLOPE  OF  LINE  THROUGH  ORIGIN 
A= Y/X 

AMS  IS  SLOPE  OF  LEAST  SQUARES  FIT 
BMS  IS  INTERCEPT  OF  LEAST  SQUATES  FIT 

DENOM=N*XX-X*X 

ANUM=N*XY-X#Y 

BNUM=Y*XX-X*XY 

AMS=ANUM/DENOM 

BMS=BNUM/DENOM 

CALCULATE  EXPERIMENTAL  FRACTIONAL  FREQ. 

DO  50  1  =  1 >  N 

50  FPREP ( I ) =REPRAT ( I ) /PREP ( I ) 

CALCULATE  MEAN  FRACTIONAL  FREQ  AND  STANDARD  DEV 

SUM=0  * 

DO  55  1  =  1 >  N 
55  SUM-SUM+FPREP ( I ) 

AN  =  N 

AMEAN=SUM/AN 
SUM2=0 . 

DO  60  1  =  1 >  N 
DEV2=FPREP<  D-AMEAN 
60  SUM2=SUM2+DEV2*DEV2 
SIGMA=SQRT ( SUM2/AN) 
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STER«*SIGMA/SGRT <  AN) 

OUTPUT  CALCULATED  DATA 
URITE<6? 143)  AMSfBMS 

143  F0RMAT< '0LEAST  SQUARES  FIT  TO  DATA  IS'/'  Y='fG14.5f'  X+  'fB14.5> 
WRITE <  6  f 1 44 )  A 

144  F0RMAT< '0LEAST  SQUARES  FIT  THROUGH  ORIGIN  IS'/'  Y»'f 
1  G14.5f'X'> 

WRITE (6f 146) AMEANf SIGMA 

146  FORMAT < '  MEAN  FRACTIONAL  FREQ= ' f G1 0 . 3 f '  STANDARD  DEU='fG10.3> 
WRITE  ( 6  f 1 47 )  STER 

147  FORMAT < '  STANDARD  ERR0R= ' f G1 0 . 3 ) 

PLOT  OF  DATA 

INITIALIZE  PLOTTING  SYSTEM 
CALL  PLOTS 
POSITION  NEW  ORIGIN 
CALL  PL0t(3.0f3.0f-3) 

DRAW  X  AXIS 

CALL  AXIS(0. ,0. f 'FLOW  VEL  M/SEC ' t -1 4 f 5 . f 0 . f 0. f 5 . ) 

DRAW  Y  AXIS 

CALL  AXIS(0.0»0. » 'REP  RATE  HZ ' f 1 1 f 5 . f 90 . f 0 . r 1 00 . ) 

DRAW  EXPERIMENTAL  DATA  POINTS 

DO  65  1  =  1 f  N 
XLINE(I)=UEL(I) 

65  YLINE< I )=REPRAT ( I ) 

XL INE ( N+l ) =0 . 

YLINE ( N+l ) =0 ♦ 

XLINE(N+2>  =5 . 0 
YLINE<N+2)=100« 

CALL  LINE(XLINEf  YLINEf  Nf1f-1f1) 

DRAW  LEAST  SQUARES  FIT 
CALL  PLOT <  0 ♦ fO. f3) 

CALL  PL0T<XLINE(N)/XLINE(N+2)f  A* XL INE (N)/YLINE(N  +  2) f2) 

DRAW  IDEAL  REP  RATE 

DO  70  1  =  1 f  N 
70  YLINE< I >=PREP( I ) 

THIS  PART  OF  THE  PROGRAM  THROWS  OUT  REDUNDANT  POINTS  AS  THEY 
CAUSE  PROBLEMS  WITH  DASHLN ♦ 

J=N 

K=N-1 

DO  75  I  =  1 f  K 
LM=N-I 

WRITE  ( 6 f  999 ) 

999  FORMAT  ( '  L00P1') 

72  IF(XLINE<LM) »NE.XLINE(LM+1) )  GO  TO  75 
DO  74  L=LM  f  25 
XLINE(L)=XLINE(L+1 ) 

74  YLINE(L>=YLINE<L+1) 
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>  197 

>  198 

>  199 

>  200 
>  201 
>  202 

>  203 

>  204 

>  205 

>  206 

>  207 

>  208 

>  209 

>  210 
>  211 
>  212 

>  213 

>  214 

>  215 

>  216 

>  217 

>  218 

>  219 

>  220 
>  221 
>  222 
♦END  OF  FILE 
♦ 


J=J-1 

C  WRITE  <  6 » 1 000 ) 

C1000  FORMAT  ('  L00P2') 

75  CONTINUE 

CALL  DASHLN()<LINEfYLINEfJf1) 

C 

C  ADD  IDENTIFICATION 

C 

CALL  SYMBOL (0.5f4.0f,2fIDf0.f16> 

CALL  NUMBER(0.5f3.5f . 2 f EDNST Y f 0 . 0 f -1 ) 

CALL  SYMBOL < 1 ♦ f3 . 5f . 2f ' J/L ' f  0 ♦ f  3 ) 

C 

C  TERMINATE  PLOT 

C 

CALL  PLOT  <  0 . fO, >999) 

STOP 

END 

C 

C  SAMPLE  DATA  DECK  FOLLOWS 

C 

SDATA  ID=/  80  '  t ' TORR '  r '  C02'f'  '» 

CAP=30  »  3E-9 f  U0LT=10000f  DIA=48.26» 

CAL=-0.189129f0»378319E-02f0.135675E-05f-0.183620E-09f 
FAN=51 . f 49, f 41 . ,40. >30. >30. >25. f25. f20. f20. f17.2f 17. f 
DISCH=12 . 4  f 12 . 2  f  9 ♦ 5 f  9 ♦ 1 f  6 ♦ f  5 . 8 f  4 . 65 f  4 . 5  f  3 ♦ 7  f  3 ♦ 6  f3 . 05  f  2 . 95 f 
N=12  SEND 


24$ 
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REPETITION  RATE  ANALYSIS 
80  TORE  C02 

CAPACITAIJCI  =  .  303E-07  FAFADS 
VOLTAGE=. 100E+05  VOLTS 
DISCHARGE  DIAM  ETEP.=  46.26MM 
ENERGY  DENSITY  =  34.78J/L 

N=  12 

CA LIE RATION  CURVE  COEFFICIENTS 

-0.  189129  0 . 37  83 19  E-02  0.  135675F-05  -0 .  183620E-09 

FAN  PERIOD  AND  INTERPULSE  PERIOD 


MSEC 

MSEC 

51.000 

12.  400 

49.000 

12.200 

4  1.  COO 

9.500 

40.  000 

9.  100 

3  0.0  00 

6.  000 

30.  000 

5.800 

25.000 

4.650 

25.000 

4.  5C0 

20.000 

3.  700 

20.  000 

3.600 

17. 200 

3.05  0 

17. 0C0 

2.950 

CALCULATED 

VALUES 

FAN  SPEED 

FLCU  VEL 

FFP  RATE  IDEAL 

RE?  RATE 

R  PM 

M/SEC 

HZ 

HZ 

1176.5 

c  .  8 

SO.  6 

121.0 

1224.5 

6.  1 

82.0 

127.  2 

1463.4 

7.  7 

105.3 

159.  1 

1500. 0 

7.9 

109.  9 

164.  1 

2  0  0C  .  0 

1 1 .3 

166.  7 

234.  9 

2000. 0 

11.3 

172.4 

234.  9 

2400. 0 

14.  2 

21  5.  1 

293.  6 

2400.0 

14.2 

222.  2 

293.6 

3000. 0 

18.4 

270.  3 

381.5 

3000.0 

18.4 

277.8 

381.  5 

3488.4 

21.7 

327.9 

450.  1 

3529.4 

22.0 

339.  0 

455.7 

LEAST  SQUARES  FIT  TO  DATA  IS 
Y=  15.836  X+  -12.569 

LEAST  SQUARES  FIT  n’FF0  UGH  ORIGIN  IS 
Y=  14.88  8  X 

MEAN  FRACTIONAL  FPFC=  0.707  STANDARD  DEV=  0.356E-01 

STANDARD  ERRCF=  0.103E-01 


